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AUS DEN BESPRECHUNGEN 
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Stationare Erwarmung indirekt axial gekiihlter elektrischer Maschinen 


Von 
H. Wen, Rheydt 


Mit 13 Textabbildungen 


(Eingegangen am 9. Mai 1960) 


Ubersicht. Die Warmestromung von den Leitern der Wicklung iiber die Hauptisolation und die Eisenteile 
des Blechpaketes (indirekte Wicklungskiihlung) wird untersucht. Die Rechnung beschrankt sich dabei auf eine 
immer haufiger vorkommende Maschinentype, bei der das Blechpaket in axialer Richtung nicht unterteilt ist. 
Innerhalb des Eisens erfolgt die Kiihlung in axial verlaufenden Kanalen, die in das Blechpaket eingestanzt 
sind. Der Warmestrom verlauft also hauptsachlich in der Blechebene, im Gegensatz zur Kiihlung durch Schlitze 
im Blechpaket. Die Ansatze werden so allgemein gehalten, daB der gréBte Teil der praktisch interessanten 
Lésungen rechnerisch erfaBt werden kann. So wird der Einflu& des Laufers auf die Standererwarmung, die 
veranderte Warmeabgabe im Wickelkopf und die mégliche Warmeabgabe zum Zabnkopf und zum Joch des 
Standers einbezogen. Bei der Beriicksichtigung der zweidimensionalen Warmestr6mung im Nutquerschnitt 
wurde zur Umgehung rechentechnischer Schwierigkeiten eine Naherung verwendet. Durch Vergleich mit einem 
genaueren Ansatz laBt sich ein Anhalt iiber den hierdurch begangenen Fehler gewinnen. Das Temperaturfeld 
wird nach Art einer Randwertsaufgabe berechnet. 


Teil I 


Warmefliisse in axialer Richtung der Maschine bleiben zunachst unberiicksichtigt. Die in 
diesem Teil gewonnenen Ergebnisse sind direkt nur auf Wicklungsanordnungen anwendbar, 
bei denen der im Eisen eingebettete Teil der Wicklung erheblich 
langer ist als der Wickelkopf. Fiir alle anderen FAlle ist eine 
Beriicksichtigung der Warmeabgabe im Wickelkopf und die hier- 
durch verursachten Warmestr6mungen in axialer Richtung not- 
wendig. Die im Teil I abgeleiteten Gleichungen gestatten die 
Angabe der resultierenden Warmeabgabezahl fiir die Wicklung 
innerhalb des Eisenpaketes. Im Teil II werden die Gleichungen 
fiir den Temperaturverlauf in axialer Richtung angegeben und 
diskutiert. 


1. Warmeabgabe in Richtung zum Nutengrund 


Der Rechnung wird ein Modell zugrunde gelegt, dessen Eigen- 
schaften dem Verhalten eines Standers einer Drehstrommaschine 
in thermischer Hinsicht weitgehend entspricht. Mit den Bezeich- 
nungen nach Bild 1 sollen die spezifischen Eisenverluste g, in 
nerhalb des Zahnbereiches 0, + x, und die spezifischen Kupferver- 
luste gy der Nut innerhalb des Nutbereiches 2 x, y, gleichmaBig 
verteilt sein (g, und gy sind raumliche Dichten). Die Warmeleit- 
fahigkeiten in den beiden Hauptrichtungen des Nutquerschnitts 
sind 2, und /,, die Warmeleitfahigkeit des Zahnes in «- Richtung /,,. 
Dank der j mn Vergleich zu seiner Héhe meist geringen Breite des 
Zahnes und der in dieser Richtung sehr kleinen WarmefluBdichten 
k6nnen Temperaturdifferenzen in y-Richtung (des Zahnes) ver- Bild. Lage des Koordinatensystems. 
nachlassigt werden. Der Temperaturabfall im Joch ist meist er- 
heblich kleiner als in den Zahnen, insbesondere dann, wenn — wie in den meisten Fallen bei 
axial gekiihlten Maschinen — das Joch von Kiihlkanalen durchsetzt ist; er soll ftir die fol- 
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[3] 
genden Rechnungen vernachlassigt werden. Bei sehr dicken, nur am AuBenrand gekitihlten 
Jochen kann durch Hinzuaddieren des entsprechenden Betrages zu den hier gewonnenen Er- 
gebnissen die Rechengenauigkeit verbessert werden. Die innerhalb des Nutenraumes definierte 
Temperaturfunktion 0,(x, y) hat der Porssonschen Differentialgleichung zu geniigen, 

Poy Oy | dn 
bP cdi I 


On (1,1) 


wobei mit Hilfe der Transformation y = 2 -7 eine Symmetrierung hinsichtlich der unab- 


hangig Veradnderlichen vorgenommen wurde. Nach Art der vorliegenden Randwertsaufgabe 
empfiehlt sich zur Losung ein Ansatz der Form 


In (%, ¥) = Poln) + Yl) + p(n) - (1,2) 


Beachtet man die im Hinblick auf y = 0 gebotene Symmetrie, so folgt damit aus der Differen- 
tialgleichung die allgemeine Lésung 


y ; ae 
By = Ono 3” oie Dy B, + sinh (v B x) " cos |/729) 


4 3C,- cosh (u 9) -c0s(ep J/72 ). (1,3) 


Die GréBen Py, B,, C, und B bedeuten Integrationskonstanten, die aus den Randbedingungen 
zu ermitteln sind. Da im allgemeinen an der Nut6offnung aus mechanischen und elektrischen 
Griinden eine Ansammlung schlecht warmeleitender Werkstoffe den Warmeiibergang in x-Rich- 
tung praktisch unterbindet, ist es naheliegend, fiir « =o die Warmestromdichte und damit 
den Differentialquotienten 00,,/0% = 0 zu setzen. In der Lésung Gl. (1,3) verschwinden dadurch 


die hinsichtlich « ungeraden Funktionsanteile. Der Versuch einer Anpassung dieser Lésung. 


an die vorhandenen Randbedingungen fiihrt auf kaum iiberwindbare rechentechnische Schwie- 


rigkeiten, solange in der Summe mehr als ein Glied beriicksichtigt wird. Wie die folgende Rech- 


nung zeigt, laBt sich jedoch ein Ansatz mit einem Glied, also der Form 


(yi Ox0— 3h y? + C- cosh (f x) - cos (4 Ve y) : (1,4) 


schon sehr gut an die Randbedingungen der Aufgabe anpassen und fiihrt zu einem recht genauen 
Ergebnis. 
Die Differentialgleichung der Zahntemperatur schreibt sich mit — = (Fs) , den von 
Z VW 
der Nut einstroémenden, auf die Volumeneinheit des Zahnes bezogenen Verlusten 
dhz Vz 2h, 00 x7\ 
ae An bghg ( ay ) (2,5) 
Die Integration dieser Differentialgleichung fiihrt mit der MaBgabe eines an der Stelle x = 0 
verschwindenden Warmestromes und der aus Gl. (1,4) zu berechnenden Ableitung der Nut- 
erwdarmung nach y an der Stelle y, auf eine Losung der Form 


Bees 9qz(1 +s) #2 2hy AEC : : Des 
b= 876 a iE / Tg cosh (B x) - sin (5 Vi » (1,6) 
Als Abkiirzung wurde die Verhiltniszahl s = = nae * als der Quotient aus Nut- und 
8 YA & 


Zahnverlusten eingefiihrt. 

Fir die Bestimmung der Integrationskonstanten der Temperaturgleichungen fiir Nut und 
Zahn spielt die den EinfluB der Nutisolation erfassende Randbedingung eine wichtige Rolle. 
An der Stelle y = y, wird mit /;, der resultierenden Warmeleitfahigkeit und A,, der Dicke der 


Isolation, als Ausdruck der quellenfreien Warmestrémung innerhalb der Isolation, die Bedingung 
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EE Ae IR ad 
Tineet ar G 3 f j qz(t +s) #4 . 
insetzen in Gl. (1,7) und Hinzuaddieren des Ausdruckes —’ eer est mit der noch zu 
2AE 
wahlenden Zahl a, erméglicht eine Trennung in die von x unabhangige Gleichung 
2 2 
qn Vi 2 hy A; qz(1 + Ss) *4 
Oxy — Oz) = "eet pee . 1 
NO Z0 2 Ay ey, I Dae a (1,8) 


und in die von « abhangige, tiber die Nuthdhe gemittelte Beziehung 


1 r 1 Hae 

A-2 { cosh (Bx)-dv= > | (a (=) ax. (1,9) 

Hierbei wurde abkiirzend 
___Cradg aN. ae ctytiy Aes an 

gos g2(t +s) O° (0 /% » f te( |) ee ‘pa am Pm) 4, ae 

eingesetzt. Die Ausfiihrung der Integration in Gl. (1,9) fiihrt auf 


A=(a4 *] iad (1,10) 


3 / sinh (f 4) 


Die Gleichsetzung der beiden Ausdriicke fiir A gestattet die Ermittlung von C zu 


q7its)t Bx, a +4]3 


2 Ag "sinh (B 44) a 
cos alk sd Mey 
hy 


mit der geschweiften Klammer nach Gl. (1,9a). Die Zahla ist keine Integrationskonstante 
und kann deshalb auch nicht aus den Randbedingungen bestimmt werden. Sie ist eine Hilfs- 
groBe, deren Einfitthrung durch die dem Verlauf der Funktion nach nur naherungsweise erfiillte 
Randbedingung (1,7) geboten ist. Obgleich iiber eine Fehlerbetrachtung die Ermittlung von a 
in Abhangigkeit von / x, nahe lage, wird zur Vermeidung weiterer Erschwerungen der Rechnung a 
als von # x, unabhangig und gegeben betrachtet. Bei zahlenmaBigen Rechnungen werden fiir a 
Betrage um 1 eingesetzt. Es liegt nahe, auch die Randbedingung am Nutengrund als im 
Mittel erfiillt vorzuschreiben. Sie lautet: 


C— ; (i501) 


Vi V1 
ae A; * (ad 1 
Ae | ( ie wee al Oy(%) dy — Iz (m4) « (5) 


ie) 


Entsprechend der ausschlieBlich in Richtung Nutengrund verlaufenden Warmestrémung be- 
stimmt man die an der Stelle x = x, herrschende Zahntemperatur, die gleich der Jochtemperatur 
ist, mit dem Produkt aus Warmeabgabezahl und kiithlendem Umfang je Nutteilung a, U, und 
der Summe aus Nut-, Zahn- und Jochverlusten zu 


(1,13) 


Aus Gl. (1,13) folgt durch Einsetzen der entsprechend spezialisierten Gl. (1,6) fiir Jz 


Un + Uz + U7 qz(a + s) x5 ' 2Ay Ae € me ae 
O75 = a, Uy 2Ar " bean VY Ay B cosh (6 x) - sin| 6 iy 92 =) (tea) 


Fiir die Integrationskonstante ?yy errechnet man aus Gl. (1,12) 


Ove = Pucmmt WG As C 7, sin | 6 oe a sinh ( x) aaa COSI Aaa )im 
(1,15) 
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Einsetzen von Gl. (1,14) in Gl. (1,8) ergibt fiir die gleiche GréBe 


ww +%zt%y | wt eG mae NPA, ad 
Oy U, 2dr 2hy AGW 


2Ay Hy ioe cosh (B%) . [, de 
LC by Ip \z Pars sin (i rv Ni]: (1,16) 


Der Vergleich der Ausdriicke GI. (1,15) und Gl. (1,16) erlaubt bei Beriicksichtigung von Gl. (1,11) 
fiir C die Bestimmung von f x, aus der transzendenten Gleichung 


Ay ctgh (8 “(2 hes =| 


ie oN a ur i 
re A i rapes 5) 
arte (0/2) V2 : oy z : ee =O (1,17) 
1+ te (e/% a eae up| \3E 


Hier wurden folgende Abkiirzungen gebraucht: 
Ay AG eee pie s s Of 34), (ieee 
Sn fier fia NE aa ~ t1+s 3(@+£ 1/3) _* a+ 4/3 ~ 


Ono = 


(1,174) 


Um eine angendherte explizite Darstellung fiir 6 x, zu erhalten, fiihrt man fiir Kreis- und 
Hyperbelfunktionen die ersten Glieder der Reihenentwicklung ein: 


il 1 (B #4)? 
Fateh Bm) & pea (t +=). 


Wegen Vee BaaK 1, darf man naherungsweise 


tg ( Vit) ~ on 


setzen, so daB man fiir 6 x, ermittelt: 
Saul 


(Bx)? T- aes rene (1,18) 


wobei zur Abkiirzung eingefiihrt wurde: 


IP A 
1+ oe = u v2 — 
See Oe Sakae Ena y 
[= 1 Me hee ae 1 f Ay 2 Ay v f (1,18 a) 
+ uv t Uv } 
Ay “3,2 3 Ay 3 A? 


Da samtliche Integrationskonstanten nun bestimmt sind, kann der erste Teil der Aufgabe als 
gelést betrachtet werden. 

Bei vielen Anwendungsfallen besteht auBer der Warmeabgabe am Joch der Maschine eine 
zusatzliche Warmesenke im Bereich des Luftspaltes. Die Nuten der Wicklung werden nur 
zum Teil gefiillt und der an der Bohrung zwischen den Zahnképfen entstehende Kanal wird 
zu Kiihlzwecken verwendet. Bei dieser Art der Zahnkiihlung flieBt der Warmestrom in zwei 
Richtungen. Andererseits legen Kiihlungsarten im Elektromaschinenbau vor, bei denen iiber 
den Luftspalt der Maschine noch Verluste der Lauferwicklung mit abgefiihrt werden. 

Von einer umfassenden Erwarmungsberechnung ist zu verlangen, daB sich auch diese Ver- 
allgemeinerungen der Randbedingungen einbauen lassen. Eine Erweiterung der allgemeinen 
Lésung des Abschnittes 1 fiihrt insofern auf erhebliche Schwierigkeiten, als sich durch die 
erhoéhte Zahl von Integrationskonstanten auBerst unhandliche Formeln, insbesondere fiir die 
Bestimmung der f-Werte, ergeben. Es liegt nahe, wenigstens die Temperaturfelder fiir den zu 
Abschnitt 1 entgegengesetzten Sonderfall einer ausschlieBlichen Warmeabgabe am Zahnkopf 
zu berechnen. Da jedoch auch hierfiir mindestens 2 Produktglieder in der Lésung fiir 0, an- 
zusetzen sind, fiihrt auch diese Rechnung auf Ausdriické betrachtlichen Umfanges, so daB auf 
ihre Wiedergabe hier verzichtet werden soll. Zur Erfassung der verallgemeinerten Randbedin- 
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gungen, also der zweiseitigen Warmeabgabe, wird im Abschnitt 2 eine Naherungslésung an- 
gegeben. Ihre Ergebnisse k6nnen mit der genaueren Rechnung des Abschnittes 1 im Sonder- 
fall einseitiger Warmeabgabe in Richtung Nutengrund kontrolliert werden. 


2. Naherungsloésung fiir erweiterte Randbedingungen 


Ahnlich wie im Abschnitt 1 wird wieder eine gleichmaBige Verteilung der Eisenverluste 
innerhalb des Zahnes zwischen +é, und — & angenommen. Es gilt &, + &, = ,; die allgemeine 
Lage des Nullpunktes der -Koordinaten ist im Bild 2 angedeutet. Die Warmeabgabe erfolgt 
an den Enden des Zahnes und wird durch entsprechende Produkte aus Warmeabgabezahl und 
kiihlendem Umfang je Nutteilung («, U, und «a, U,) beriick- 
sichtigt. An der Stelle = &, sind Zahn- und Jochtemperatur 
identisch; das Temperaturgefalle im Joch wird wieder ver- 
nachlassigt. 

Bei Kiihlung in mehreren Kanalreihen miissen deren Stré- 
mungswiderstande ungefahr gleich groB gewahlt werden, wenn 064 Uy 
ein giinstiger Kiihleffekt erzielt werden soll. Unter dieser Vor- 
aussetzung darf man annehmen, daB die Kiihlmittelerwarmung 
in den verschiedenen Kandalen nicht sehr unterschiedlich ist, so 
daB die Umgebungstemperatur an den Kiihlflachen naherungs- 
weise gleich groB (und zunachst gleich Null) gesetzt werden kann. 

Der Einflu8 des Laufers auf die Standerwicklung wird da- 
durch beriicksichtigt, daB an der Stelle & = —&, auch ein Teil 
der Lauferverluste kv, mit k<1 dem Zahn zugefiihrt wird. 
Durch die Anwendung eines vereinfachten Lésungsansatzes 
werden Erleichterungen in der Rechentechnik und Verein- 
fachungen im Ergebnis erzielt. 

Fiir kleine Argumente der Kreis- und Hyperbelfunktionen 
der allgemeinen Losung (1,4) kénnen diese durch die ersten 
Glieder ihrer Potenzreihenentwicklung angenahert werden. Im 
Produkt werden die Glieder mit xy? als kleine GroBen vierter 
Ordnung vernachlassigt, was auf eine Aufspaltung in zwei nur 
jeweils von einer Veranderlichen abhangige Funktionsanteile 
fiihrt. 

Zum selben Ziel fiihrt der Lésungsansatz 


On(%, ¥) = f(*) + Bly) - (2,1) 
Versteht man unter 7, + vy und 7, - vy jene Teilbetrage der Nutverluste, die in x- baw. y-Richtung 
abgefiihrt werden, so darf man mit 


Bild 2, Lage des Koordinatensystems. 


Vow ty oh (232) 
die Differentialgleichung der Nuttemperatur 
Oro aro 
A, as A, at (7, r,) In == 0 (2,3) 
separieren und erhalt als allgemeine Lésung den relativ einfachen Ausdruck 
%%9IN 95 Vy IN 
Oy = Ono ys (ae ee (2,4) 


In Gl. (2,4) sind die ungeraden Funktionsanteile bereits unterdrtickt, wie es die Symmetrie 
zur Mittellinie der Nut und die Randbedingung eines verschwindend kleinen Warmetibergangs 
an der Stelle « = o erfordern. ne 
Die Differentialgleichung der Zahntemperatur (1,5) wird mit der Ableitung (2) aus 
24N 41 YN 7 


Gl. (2,4) und der Abktirzung s = Si aa 


oe -- o. (1 + 7%,s) =0. (2,5) 
E 


dé? 
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Die allgemeine Lésung von Gl. (2,5) lautet: 


gz(1 + %y s) 
9, = Bz) +b E—— 


2Ap 


ip (2,6) 


Beginnt man zur Anpassung an die Randbedingungen mit Gl. (2,6) und fordert fiir § = §,:0, = 0, 


so erhalt man 
qz( (1 + + %y s) 


Oy = Oy = OEE) a ee ae (2,7) 


Die Randtemperatur #, bestimmt man aus der Gleichung der Verlustbilanz an der Stelle &, 


21, In %1 91 — Ag 57 O7(E,) + vp = 0% U,-9,. 
Hieraus ergibt sich: 


Agbz , bzaz(tt+ trys), | 2% INMM | YY 
ar ; 2,8 
Uiheeaee Be Uy 0, Us; a4 pe Olay (2,8) 
Am Zahnkopf, fiir € = —&,, wird #, = 0. Es wird berechnet aus der Verlustbilanz an dieser 


Stelle, 


bg Ag Pz(— &) + hv, = o% Uz De 
und ergibt sich zu 


Anbz | bz 4z(t + YyS) 


Rv, 
Dy a by Oy om “a a, UF Gat Oe Ue © (2,9) 


Die fiir &; und — , angeschriebenen Zahnerwarmungen nach Gl. (2,6) ergeben eine Differenz 


qz(1 + %y s) 
B, — By = b, x, — 2 se Ge (2,10) 
Vergleicht man Gl. (2,10) mit der aus den Gln. (2,8) und (2,9) resultierenden Differenz, so laBt 
U. 
sich daraus eine Gleichung fiir 6, herleiten. Sie lautet mit der Abkiirzung » = =! : 
2 2 


(a + % s) 2bzhr 


Arb 
bya(1 + oe (1 +9)) = |e 4 OA, 9B) + 


14 Un + V7 —kYV dy 
a Uy 


(2314) 


Zu einer identischen Aussage iiber 6, gelangt man durch Heranziehung der Verlustbilanz der 
gesamten Anordnung 


Oy Oy Cy Onis Up = (2,12) 


wobei auf der rechten Seite von Gl. (2,12) die Summe aller Verluste steht, und z. B. die Gln. 
(2,8) und (2,9) in (2,21) eingesetzt werden. 

Es sind also bereits alle den Zahn betreffenden Randbedingungen erfiillt. Da die Zahl der 
noch zu bestimmenden Integrationskonstanten um 1 gr6Ber ist als die Zahl der Randbedin- 
gungen, kann man 0}, = o setzen. Der Nullpunkt der -Koordinaten wird hiermit verschieblich 
zwischen + 0o gegentiber x = 0. Mit 6, = o berechnet man aus Gl. (2,11) bei Einfiihrung der 
bezogenen GréBen 

gy Ss a, Uy % Vp—Vh vy 


1 ay ; 2 wy i By 2Ar bz ? UZ 


den normierten Abstand des Koordinaten-Nullpunktes vom Nutengrund zu 


W, = = Pid th = tas ae 


14+2R+y (1+7ys) (1428, +>)" (2,13) 


Zur Bestimmung der noch unbekannten Integrationskonstanten 7, und vy, stehen zwei Rand- 
bedingungen zu Verfiigung. Langs y = y, ist 


a 1 
=i =A, (cee). dp = { (Oxy(04) — 94) dé. Gas 


2 Sa 
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Bei Beachtung der Transformationsformel zwischen den beiden Koordinatensystemen x = & + &, 
liefert die durch Gl. (2,14) vorgeschriebene Rechnung 


2 2 2 
eR RE Ln Ay Ai "ON *1 | Q2(t + %yS) %y 1 
Oon a, os 2 dy (: ny 2 A; yy = é le T Wy ——]. (2,15) 


Langs des Randes x = x, formuliert man die zu Gl. (2,14) analoge Bedingung 


Vy ) 
1 he OO n 1 fa 
yi | A; 4,( ox ).- dy =; | (Ix (%4) = ;) dy . (2,16) 
Die Ausrechnung fiihrt auf 
ty In Vi, Ye IN Az: 
On — 8, = ae a aes ( | 2 7 ak (2,17) 


Die Gln. (2,15) und (2,17) zusammen mit Gl. (2,13) fiir w, erlauben die Berechnung von 7,. Es 
entsteht ein Ausdruck der Form 


R, 4 an oy; C 
y STP aeny) 
pg a: (2,18) 
Oa Sagara 
Darin bedeuten die Abkiirzungen 

pela 32S), Rae 4337S 
34401 A; 4 4 3Ay M1 Ani” 
pie ee Vr 
wll b= F in 
Regen AS 130494 2h) % 


In Zahlentafel 1 sind die Gleichungen fiir w, und 7, fiir einige Sonderfalle der Anwendung zu- 
sammengestellt. 

Zeile 1 bezieht sich auf die Voraussetzung ausschlieBlicher Warmeabgabe im Standerriicken, 
wie es bei Maschinen kleiner und mittlerer Leistung gegeben ist. Hierfiir wird mit ao, U, — 0, 
also vernachlassigbarer Warmeabgabe am Luftspalt, y->0o, foo und 4, = 2/3. 


Zahlentafel 1. 


W, x 
251 -- Ss Ri) 
0%, U, Ry +> one 
Rv, 3x Ss uN 
1 Ky U, =O 1+ 
vz (1 + Yy S) Pee po ae 
% “/ 
0, U, so ees 
sm A, Up = 0} ave J 8 
A 1 
kv, =0 Re Ry 
= | 
| is t a -b os v7 
0 %,Uy + U7 y Sas UN 
3 << & U,=0| \ 
vz (1 + %yS) Re eye 
a y 
0, U TESS: v 1 
eed owe Sats ee eee 
4 (+) a= a ae Se 
pe nl as (1 + %y S) 2 (By a 1) R ER saprig Seas 
# v6 Sar(Bp 1) 
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w, wird fiir diesen Sonderfall gréBer als 1, d. h. der Nullpunkt der &-Koordinaten riickt aus 
dem Zahnbereich heraus, wenn die auf den Stander tibertragenen Lauferverluste k v, nennens- 
werte Betrige annehmen. Hierdurch ergibt sich andererseits eine VergréBerung des Anteiles 7, 
der Verluste, die iiber den Nutengrund direkt an das Joch abgeleitet werden. Bei gentigender 
Beliiftung des Laufers vereinfachen sich die Gleichungen weiter; es kann dann k v, = 0 gesetze 
werden, was in Zeile 2 aufgefiihrt ist. Da ebenfalls 7, > 0 ist, liegt der heiBeste Punkt der 
Nut in unmittelbarer Nachbarschaft des Keiles. Die Zahntemperatur fallt zum Joch hin ab. 
Zu einer weiteren Vereinfachung der Gleichungen fiihrt der Ubergang zu verschwindend kleinen : 
Eisenverlusten im Zahn, also s > 00; ein Grenzfall welcher der KurzschluBerwarmung elektrischer 
Maschinen weitgehend entspricht. In Zeile 3 wurden fiir den Sonderfall ausschlieBlicher Warme- 
abgabe am Zahnkopf die Gln. fiir 7, und w, angegeben. Voraussetzung fiir diesen Grenzfall ist 
a, U, = 0, »=0, B, = 0, 4: = — 1/3. 7, nimmt dann fiir alle praktisch denkbaren Faille 
negative Werte an und auch der Nullpunkt der Koordinaten kann auBerhalb des Zahnes in 
den Bereich des Joches fallen. Der heiBeste Punkt der Nut liegt nun in der Nachbarschaft 
des Nutengrundes. Zeile 4 gilt dem Sonderfall gleich intensiver Warmeabgabe an Joch und 
Zahnkopf, also y = 1. 

In den Bildern 3 und 4 ist der Verlauf der Nuttemperatur fiir y = o und der Zahntemperatur 
iiber der normierten Nuthdhe zu einigen in der Tabelle enthaltenen Beispielen aufgetragen. 
Bild 3 hat den Sonderfall 2 zur Voraussetzung und wurde fiir folgende Werte berechnet: 


S=4, pHi, v= 0,14, UZ Ciksme 


=—, (Poa We 


(Erklarung der Abkiirzungen in Zahlentafel 2.) 


0 0,2 0,4 0,6 0,8 40 0 0,2 0,4 0,6 0,8 40 
8G) a, 
Dy, Dy 


Bild 3. Verlauf von Nuterwarmung Oy und Zahnerwairmung 07 


Bild 4. Verlauf N a a 
fiir Warmeabgabe am Joch (vy = oo); (Zahlenbeispiel). y Nasties heabe aca Uakelore (/ee Grote 


fiir Warmeabgabe am Zahnkopf (v = o); (Zahlenbeispiel). 


Die eingetragenen @-Werte sind Verhaltniszahlen und durch Division der Temperaturen # 
qZ*4 


mit re gewonnen worden. Der Temperaturverlauf in Nut und Zahn bei einem Kiihlschema 


nach Sonderfall 3 der Tafel1 ist im Bild 4 wiedergegeben. Die ausschlieBliche Kihlung am 

Zahnkopf erweist sich als sehr viel ungiinstiger als die Kithlung am Standerriicken. Es wurden 

wieder die gleichen Zahlenwerte wie bei Bild 3 zugrunde gelegt und f, = 0, fs = 1 gesetzt. 
Bei sich zum Zahnkopf hin verjiingenden Zahnen ergeben sich fiir die Zahnkopfkiihlung 


nears in dieser Rechnung nicht zum Ausdruck kommende Nachteile gegentiber der Joch- 
<tihlung. 


Ne vee eres 
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Entsprechend Sonderfall 5 der Tafel 1 wurde wieder fiir die gleichen Parameter-Werte 
der Temperaturverlauf fiir y= 1 und gleiche Gesamtoberflache fiir die Warmeabgabe, also 
B, = By = 0,5, berechnet und im Bild 5 dargestellt. Durch die annahernd symmetrische Warme- 
abgabe nach beiden Richtungen erlangt man erheblich 
glinstigere Temperaturen als in den vorher betrachteten 40 
Fallen. 

Zahlentafel 2 enthalt die Zusammenstellung der auf 


2 
GZ % : s 
aa bezogenen Temperaturen. Die Zahnerwarmungen 
E 


sind in Abhangigkeit von der Koordinate w = &/x, an- 
gegeben. 

Zur Kontrolle des zahlenmaBigen Einflusses der S 
iiber den Nutengrund abstrémenden Warmemenge © 
werden in Abhangigkeit vom Verhiltnis v = y,/x, fiir 
den Parameter s = 4 die ,,innere‘‘ maximale und die 2 
mittlere Wicklungserwarmung fiir y = 0 mit 


6 aly 0 0,2 04 06 08 4,0 
berechnet. Sie werden bezogen auf jene Wicklungs- Cae 


erwarmung, die sich ohne Beriicksichtigung der Warme- Bild. Verlauf von Nuterwarmung @y und Zahnerwar- 


o3 : Oz fi veiseitige Warmeabgabe (vy = 1); 
abgabe iiber den Nutengrund, also mit 7, =0 zu caesar oar Sia 


(Oxo— O)y = 8 P= (1+ 24) + (1 +9) (2,19) 


ermitteln la8t. Fiir v = 0,14 wurden die Hilfsveranderlichen # und # zu 2,0 bzw. 1,0 gewahlt 
und fiir verschiedene Nutbreiten mit 1/v bzw. v verandert. Dies entspricht einem Vergleich 
verschieden breiter, aber gleich tiefer Nuten mit konstantem Verhaltnis Nut- zu Zahnbreite 
und s = const. AuBerdem wurde vernachlassigbar kleiner Laufereinflu8 und «, - U, = 0 voraus- 
gesetzt. Der Verlauf der Kurven im Bild 6 bestatigt die Erwartung eines relativ groBen Ein- 
flusses des itiber den Nutengrund abflieBenden Warmestromes, insbesondere bei breiten Nuten 
und bei relativ schlechter Leitfahigkeit des Zahnes sowohl hinsichtlich der maximalen (Kurve A) 
als auch der mittleren (Kurve B) Erwirmungen. Die Anwendung der Gl. (2,19) zur angenaherten 
Berechnung der Nuterwarmung fiihrt jedoch auch bei sehr schmalen Nuten auf fehlerhafte 
Ergebnisse. Fiir vy = 0 wird der Mittelwert von dieser Gleichung um 15% zu hoch, der Maximal- 
wert um 10% zu niedrig angegeben. AuBerdem tritt, wie weitere Rechnungen zeigen, bei An- 


Zahlentafel 2. 


A 2 2 
Ox = Ovo ae Bae 


BR Tb, NOR 2a, 


@z = O,+4+ (1+ ty 5) (wi — w?), 
Oz = O,+ (1+ rs) (wy — w%) , 


Ono = O, + rysp—(1 2) (a +79) (o— 


0, = 


i 
ERS ee aes 5 ; 


ne 
Bx 


zal 
v= — 
yy 


Le Glel(2ds))s, 
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ordnungen mit kleinen Verhaltniszahlen s eine stirkere Reduzierung der Erwarmungen auf als 
bei Anordnungen mit grofen s-Werten. Im Bild 6 sind auBer diesen (ausgezogenen) Kurven 
fiir gleiche Voraussetzungen nach den Gleichungen des Abschnitts 1 die entsprechenden bezogenen 
Temperaturen (strichliert) eingetragen. Der Vergleich mit den genaueren Ergebnissen des Ab- 
schnitts 1 laBt erkennen, daB die durch Anwendung der Naherungslosung zu erwartenden 
Fehler, insbesondere bei den praktisch sehr oft vorkommenden relativ schlanken Nutformen nur 
gering sind. Es wird auch deutlich, da® der Fehler der Naherungslésung fast ausschlieBlich 
im Sinne etwas zu hoher Erwarmungen gerichtet ist, so daB bei ihrer Anwendung eine gewisse 
Sicherheit gegeben ist. Wie ein Vergleich der gemittelten Randbedingung Gl. (2,14) mit der 
értlichen Warmeiibergangsbedingung zeigt, sind fiir abnehmende »-Werte langsam zunehmende 
Fehler zu erwarten. 


Bild 6. Einflu8 der Beriicksichtigung der zweidimensionalen 

Warmestrémung in der Nut fiir y= oo. Kurve A, Maximal- 

wert, Kurve B, Mittelwert. Naherung nach Abschn. 2, 
——— Rechnung Abschn. 1; (Zahlenbeispiel). 


Bild 7 (rechts). Lage des Koordinatensystems. 


3. Anwendung des Naherungsverfahrens auf eine Zweischicht-Wicklung 


Drehstrommaschinen werden sehr haufig mit Zweischicht-Wicklungen ausgefiihrt. Sie unter- 
scheiden sich thermisch insofern von den Einschicht-Wicklungen, als sie in der Mitte der Nut- 
hohe eine oder mehrere elektrisch isolierende Trennschichten aufweisen, die insbesondere bei 
Hochspannungswicklungen den Warmeaustausch zwischen den beiden Nuthalften praktisch 
unterbinden. Die Ableitung der Temperaturgleichungen geschah unter der Voraussetzung eines 
vernachlassigbar kleinen Warmeaustausches zwischen den beiden Schichten und verschwindend 
kleiner Warmeabgabe zum Keil. Mit diesen Annahmen ergibt sich bei Anwendung des Nahe- 
rungsverfahrens ein konstanter Temperaturverlauf in x-Richtung innerhalb der Oberschicht. 
Wie Bild 7 zeigt, wurde die Verjiingung des Zahnes in erster Naherung durch die Festlegung 
zweier verschiedener Zahnbreiten im Bereich von Unter- und Oberschicht der Nut beriicksichtigt. 
Die zur Unterschicht gehérigen GréBen wurden durch den Index 1 gekennzeichnet; die GréBen 
der Oberschicht sind mit dem Index 2 versehen. Der Gang der Rechnung verlauft ahnlich wie 
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im Abschnitt 2 fiir die Einschichtwicklung. Die Kiihlmitteltemperatur wurde wie bisher auf 
beiden Seiten der Nut gleich Null gesetzt. Als zusdtzliche Randbedingungen an der Stelle = &, 
also der Mitte des Zahnes, werden ein sprungfreier Verlauf der Zahntemperatur und eine Uber- 
einstimmung der Warmestrome auf beiden Seiten der Trennflache gefordert. Zur Vereinfachung 
wurde der 0-Punkt der x, y-Koordinaten in die Mitte zwischen die beiden Wicklungshalften gelegt 
und die Nutverluste als gleichmaBig tiber den Bereich 2 x, y, verteilt angenommen. Aus der For- 
derung konstanter Gesamtverluste berechnet man die reduzierten Verluste je Volumeneinheit zu 


kr-og?+ X Fr 


2% 


v= 


Darin bedeuten kp den Wirbelstromfaktor, @ den spezifischen Widerstand bei der mittleren Wick- 
lungstemperatur, g die elektrische Stromdichte und 2’ F, die Summe aller blanken Leiterquer- 
schnitte einer Schicht. gy wird auBerdem fiir beide Schichten gleich groB vorausgesetzt, was fiir 
den groBten Teil aller Maschinen naherungsweise zutrifft. Fiir die PEGS aay der raum- 
lichen Verteilung der Zahneisenveluste bestehen zwei Moglichkeiten: 


Man kann annehmen, daB sie in den beiden verschieden breiten Zahnabschnitten, die der 
Unter- bzw. Oberschicht benachbart sind, abschnittsweise gleichmaBig verteilt sind. Die ge- 
samten Eisenverluste des Zahnes sind dann auf diese Abschnitte umzurechnen. Die im Zahnkopf 
erzeugten Eisenverluste werden zu den Verlusten des Abschnittes 2 addiert. 


Bei groBen Verhaltnissen bz ,/bz, und hohen Kiihlkanalen am Zahnkopf kénnen die dort konzen- 
trierten Eisenverluste betrachtliche Anteile der gesamten Zahnverluste bilden. Fiir solche Falle 
laBt der Ansatz die Méglichkeit, die Zahnkopfeisenverluste rechnerisch als vom Laufer her tiber- 
tragene Verluste zu betrachten und sie zu diesen zu addieren. Sie werden dadurch an den Rand 


des Zahnbereiches 2 (§ = —&,) verlegt, wodurch die Genauigkeit der Rechnung erheblich ver- 
bessert wird. 
: pe ay : B 
Die auf 7 bezogenen Nut- und Zahntemperaturen sind in Zahlentafel 3 angegeben. 
E 


Zahlentafel 3. 


6 Ay ia’ v fy? 
= ONio0 — 1x5, Pi hy ae %y —y Sy Pi 8 \y, , 


1+ 2u hy 1+ yy" Sy 
= O1 + Ny S10 a tt er eee Oe: 
w Ss v 
Ws ee en JT 
T 1% T ’ 
By 2B, 20z,/, 


= 6,4 (1+ Ny) (wi — wv) > 


= (A 1S) 


IZ2 y 
= On S- 2 = 
cm qZ1 aP 8 1 


= = @; ae So Pp V 
Z1 


1 
,= O,+ oe ~) 


A 1 boos uv 
roel a+5)(4—2)+04n(2 x4)( 2 22) : 


qz4 1+ %1ySy UZ9 


0z2 = 8, + iBa 2 ee (2 + ry) Wg (Wy + w) + (1 + 59) (a — wt] , 


UZ 2 


2 
49711 
2Ar : 


Bezugstemperatur ist 
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; g 
Die Zahnerwarmungen wurden in Abhangigkeit von der normierten Koordinate w = a an- 


: : ; & 
gegeben. Der Koordinaten-Nullpunkt wird mit w, = ae berechnet aus 
1 


20, = M4 ee aK ; (3,1) 
Darin bedeuten: 
ole) 
ny = = 7 Dali fortotA(r+ 2), as 
3,1a 
(¥z_ + Uno) ty | a + vk v,—vyz 
i 


UZy 


Die in den Gleichungen der Zahlentafel 3 vorkommenden bezogenen Abstande des Zahnmittel- 
punktes und des Zahnkopfes vom Koordinaten-Nullpunkt berechnen sich aus den Beziehungen 


1 
+ W=1, od Neer Seng eat (3,2) 


Der Anteil der Kupferverluste der Unterschicht, welcher iiber den Nutengrund an das Joch ab- _ 
gefiihrt wird, ist zu bestimmen aus 


Bye: 
Vix i, 1 1 (3,3) 
Tee ae In a iy —— ae 
Wie bei der Einschichtwicklung ist 
Mies Kiy — 1. (3,4) 
Es sind noch einige Abkiirzungen aus Gl. (3,3) und den Gln. der Zahlentafel 3 zu erklaren: 
Boeke AA; eis As vA; 
R,=— (2 - 2S), R, =—— 1 = —}, 
Z 34g 1 A; 4 ? 3 Ay V1 +3 A; v1 
UN, _ 9N1 291 UN® UN4 In1 21 
3 = nee nee Soo — — , 
1 = 9z, Oz Uz. YZ Iza bz 
(3,5) 
B Uy B __ Oy Us % p _ Aba p _ An Ze 
_ Ap bz, 2 Ag bz» u Ay ay : 2 Ay % 
SL Ay A; 
Fy, Soll ae Fe 
Oy Uz Mia a4 


Um einen Vergleich mit der Einschichtwicklung durchzufiihren, wird fiir bz, = bz,, also fiir kon- 
stante Zahnbreite und fiir den bereits im Abschnitt 2 zugrunde gelegten Sonderfall ausschlieB- 
licher Warmeabfuhr zum Joch (vy > oo) und 


A 
= =1, RM = 0, Pi =2, Ue O20n 1 — 2 geet 
y 


der Verlauf von Nut- und Zahntemperatur berechnet. Mit den angefiihrten Gleichungen ermittelt 
man 


17 — 0.4005 sy = 0,504, Wy == 1,24", W, = — 0,24, Wy == 0,74. 


Da die Summe der abzufiihrenden Verluste gleich jenen der Einschichtwicklung angesetzt wurden, 
wird auch 0, = 6,33, wie im Abschnitt 2. Die Warmeabfuhr von der Oberschicht findet nur in 
Richtung zum Zahn statt, es entstehen im Vergleich zur Einschichtwicklung héhere Erwar- 
mungen; der Maximalwert ist fiir die gemachten Voraussetzungen um 14% groBer. Bild 8 zeigt 
den Verlauf von Nut- und Zahnerwarmung. Strichliert ist der Verlauf der Nuterwarmung fiir die 
Einschichtwicklung nach Abschnitt 2 eingetragen. AuBerdem ist strichpunktiert auch der Verlauf 
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der Nuterwaérmung fiir die Voraussetzung 7,, = 0, also unter Vernachlassigung des Warme- 
flusses tiber den Nutengrund, wiedergegeben. In vielen Fallen, insbesondere dort, wo es sich um 
schlanke Nuten handelt, wird man (mit etwas mehr Berechtigung als bei der Einschicht- 
wicklung) von dieser Vereinfachung Gebrauch machen diirfen. Ausgehend von den Zahlen- 
werten, die dem Beispiel zugrunde gelegt wurden, ist im Bildg der Zusammenhang zwischen 
hochster Nuterwarmung O,,, fiir 7,,—= 0 und dem Zahnbreitenverhaltnis e = bz ,/bz, wieder- 
gegeben. 


14 


12 


10 


—1,0 0 02 OF 06 08 1,0 1,0 1,2 1,4 1,6 18 2,0 


0 ae oe 
xy € 


Bild 8. Nut- und Zahnerwarmung ftir Zweischichtwicklung, Bild 9. Maximale Nuterwarmung als Funktion der Zahnverbreite- 
—-+— Nuterwarmung fiir ry y= 0, —-— Nuterwarmung fir Ein- rung € fiir 7yz = 0; (Zahlenbeispiel). 
schichtwicklung; (Zahlenbeispiel). 


Teil II 


Im Teil I wurde die radiale Warmestr6mung innerhalb des aktiven Volumens berechnet. Die 
iiber das Eisenblechpaket hinausragende Wicklung gibt einen Teil ihrer Verluste tiber die Wick- 
lungsisolation an das Kiihlgas ab. Im allgemeinen unterscheidet sich also die Warmeabgabe im 
Bereich des Wickelkopfes und des Blechpaketes. In vielen Fallen, insbesondere bei Maschinen 
groBer Leistung und Drehzahl, ist die Warmeabgabe im Bereich des Wickelkopfes besser als im 
Bereich des Eisens. Jedoch auch der umgekehrte Fall ist bei besonderen Umstanden anzutreffen. 
Die exakte mathematische Formulierung der Aufgabe wiirde die Losung eines dreidimensionalen 
Problems verlangen. Statt dessen wird zur angenaherten Berechnung der axialen Temperatur- 
verteilung von der iiber den Nutquerschnitt gemittelten Wicklungstemperatur nach Teil I aus- 
gegangen, wodurch der Rechnungsgang erheblich vereinfacht wird. 


Die Berechnung der in radialer Richtung erfolgenden Warmestrémung ergab (unter Vernach- 
lassigung der axialen Warmestromung) fiir die Nuterwarmung eine Gleichung folgender Form: 


ipa ee oe. (4,1) 


Pana und &#, berechnen sich aus der innerhalb des Nutquerschnitts gemittelten Nuttemperatur 
NES 1 

nach den Gleichungen des Teiles I und entsprechender Ordnung der Glieder. YJ, ist die Tem- 
peratur des Kihlmittels. So erhalt man z. B. fiir Zweischicht-Wicklungen und den hinsichtlich 
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LT 


seiner Genauigkeit oft ausreichenden Sonderfall vernachlassigbar kleiner Warmestrémung tiber 
den Nutengrund (7,, = 0) 


f Bi 
BG Ee ee ; 
aii. 7 Tee 2 (2+ ~)— 2%: bat (1 +2024] 
oy U, ebzorz later + B, (1+ ¢)\2 ae Gy 3 ae 6 ; 
. pb : VR VU, “1 

9! ty ee Cero ane ip Op, (2 oe ite ite v7é 

Seay = = ; 

ee 2(1 +» +B (a +8) a a Bi 12 4 | vzePyJ 

(4,2) 


Die Warmestromung in axialer Richtung (Bild 10) hat zur Folge, daB iiber den Eisenweg der im 
allgemeinen unterschiedliche Betrag o - vy mit o = o(z) abgefiihrt wird, wahrend in axialer Rich- 
tung der Teilbetrag (1— o)-v, der Kupferverluste zu kalteren Teilen der Wicklung abstromt. 
Fiir die Nuttemperatur und die noch unbekannte Hilfsgré8e o konnen zwei Gleichungen aufge- 
stellt werden: 


( ) ig ! 
Iul2) = Ge rge + 8 + cle)» (4,3) 
29 : 
et gee (4,4) 


Ay ist die Warmeleitfahigkeit der Wicklung in z-Richtung, fy ist der Kupferquerschnitt einer ~ 
Wicklungsschicht. Gl. (4,4) kann als Ausdruck des Warmegleichgewichts in Langsrichtung der 
Wicklung betrachtet werden, wahrend Gl. (4,3) sich aus den Untersuchungen tiber die radiale 
Warmestré6mung von der Wicklung iiber das Eisen ergab. 

Um das Gleichungssystem (4,3) und (4,4) nach den Unbekannten #y, o und %, auflésen zu 
kénnen, wird noch eine weitere Gleichung bendtigt. Sie ist gegeben in Form der Differential- 
gleichung der Gaserwaéarmung innerhalb 
der Kiihlkanale. Ist % die Oberflachen- 
temperatur der Kanalwand (#, oder #,) und 

a U, 
a= Cp yur, 
zahl, U, dem Kanalumfang, c, - y der spe- 
zifischen Warme des Kiihlmittels bezogen 
auf die Volumeneinheit bei konstantem 
Druck, v der Kiihlmittelgeschwindigkeit 
und f, dem Kanalquerschnitt, so ist die 
Kiihlgastemperatur zu berechnen aus 


Yor V, 


mit « der Warmeiibergangs- 


Bild 10, Lage der z-Koordinaten und Annahmen tiber } 9 4 
Verlauf der Kithlmitteltemperatur 9g. ae Pi, (Yo ins JG) Ss (4,5) 


Da die geschlossene Lésung der drei gekoppelten Gleichungen auf sehr umstandliche Ergebnisse 
fiihrt, interessiert fiir die Anwendung die den tatsachlichen Verhaltnissen schon sehr nahe kom- 
mende Voraussetzung eines linearen Zusammenhangs zwischen der Gaserwiarmung und der 
axialen Koordinate z (siehe Bild 10). 


Oh eae (4,6) 


Diese Annahme fiihrt auf noch gut iibersehbare Lésungen; sie enthalt auBerdem den fiir sehr 
groBe Kiihlgasmengen annahernd erreichbaren Grenzfall einer von z unabhangigen Kiihlmittel- 
erwarmung, 

Eliminiert man aus den beiden Gln. (4,3) und (4,4) die Hilfsfunktion o und setzt fiir (2) 
die Gl. (4,6) ein, so entsteht fiir Jy die gew6hnliche Differentialgleichung zweiter Ordnung mit 
konstanten Koeffizienten und einem in 2 linearen Stérglied. 

Oy 


dt 


Bey a en 8 2 ov Uy gs On Uy a: 
ay Dy ORR, A ON et a 


; 


ee er ee 
% 


- 


a ae 


">= 4% @ 
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Die aus der partikularen Lésung und dem Integral der homogenen Differentialgleichung (4,5) zu- 
sammengesetzte Lésung lautet: 


bx 
ile = 2 
an 


+ kz-+ A, cosh (ay 2) + A, sinh (a, 2) . (4,8) 


Um die Integrationskonstanten der Nuterwarmung nach Gl. (4,8) zu bestimmen, miissen zuerst 
die Temperaturgleichungen der Wickelkopfanteile hergeleitet werden. 

Betrachtet man den linken Wickelkopf als im Bereich negativer z-Koordinaten am Kiihlgas- 
eintritt legend, so darf man ihm eine Kiihlmitteltemperatur 0,, — 9, zuordnen. Der am Kiihl- 
mittelaustritt liegende Wickelkopf erhalt die Kiithlmitteltemperatur 0,, + 0,. Die gesamte 
Kihlmittelerwarmung 2 3, berechnet man aus den vom Kiihlgas abzufiihrenden Verlusten 2 V, 
der spezifischen Warme des Kihlmittels bei konstantem Druck ec, + y und dessen Volumen je Zeit- 
einheit Q zu 


2 V 
20, = a (4,9) 
Fiir die mit 1 bzw. 2 gekennzeichneten Wickelképfe gelten also die Differentialgleichungen 

a4, 

— a2 [d,, — (8... —9,)] + 0, =0, 

a2), 

dz2 — as [Ps (one o,)| b, = Ox (4,10) 
OF es Xs (Gh +4 UN 
aes eo sri ST yo erie 


Hier wurde bereits vorausgesetzt, daB die Warmeabgabezahl und die Verluste je Volumeneinheit 
beider Wickelkopfe gleich sind und da8 die Warmeleitfahigkeit in z-Richtung und der Wicklungs- 
querschnitt mit den Werten innerhalb der Nut iibereinstimmen. 


Die allgemeinen Lésungen fiir die Wickelkopftemperaturen lauten: 


b 
O,, = 0,,—0, + ee + B, cosh a, (z + 2), | 


‘ (4,11) 


O.. = O19 +O, + pias B, cosh a, (2 — 2,) . 


s 


Die Gln. (4,11) sind der Bedingung verschwindend kleinen Warmestromes in z-Richtung 
(dd, ,/dz = 0, d&,,/dz = 0) an den Enden der Wicklung (z = + z,) angepaBt. Diese Bedingung ist 
allerdings nur dann gerechtfertigt, wenn Unter- und Oberschicht der Wicklung gleich gut gekiihlt 
werden. Da bei der hier vorliegenden Berechnung der axialen Temperaturverteilung von der 
iiber beide Schichten gemittelten Erwarmung ausgegangen wurde, ist die verwendete Form 
(4,11) der Lésung naheliegend. 

Zur Berechnung der Integrationskonstanten A,, A,, 6, und B, stehen an den Stellen - z, die 
Randbedingungen 0, = #, und d#,/dz = dé./dz zur Verfiigung. Die Auflosung der Gleichungen 
fiihrt auf folgende Werte fiir die Konstanten: 


1 UN UN ; Or, —k’ 2, 
es [ a.) , A, a 


i 7 = T Fi 
n+ cosh dy 2, \Oy Uy Oe Ur. m+ sinh ay 2, ’ 


1 Ry _ aN Or, —k’ 2, _ aN tgh ay 2 UN UN ; 
1 — — = | 1 ae =i Sie On d 
sinh a, (2, — 2,) |@, a, m:tghay % as n an Uy OU 
1 k ay O8,—kh’'s% | ay tghay% [on UN | 0.) (4,12) 
a is eT, , 
2 sinh a, (2, — %) a, a, m-tghayy | 4, n oy Uy ge ah he 
an tgh ay 2 an 
i= te = i = te Che ene teh a. (2 — z 
i ' a, tgha, (2,—%) ’ ue a; Ore Naat gh 4, (%, 1)» 


if 


a 


a, 2,: tgh a, 7 : 
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Bild 11 zeigt den typischen Verlauf der axialen Temperaturverteilung an einem Beispiel fiir ver- 
schieden groBe Betrage der Kiihlmittelerwarmung #,. Die Temperaturgleichungen wurden auf 
PL aN ompfiehlt, 


den Wert 6,/a%; bezogen, so daB sich als Parameter der Kurven das Verhaltnis 
das zwischen o und 0,5 variiert wurde. Als weitere Voraussetzung wurde 


g. a 
1 N os 

SO Ne, ana H1, = 0.08 Or = 0,7 
$s 


z as 


angenommen. 


oz, 2; 0 2 2s 


Le 


Bild 11, Axialer Verlauf der Wicklungserwaérmung fiir verschiedene Gaserwarmungen; (Zahlenbeispiel). 


Im Sonderfall vernachlassigbar kleiner Werte fiir #,, wie sie bei Kiihlungen durch groBe 
Kihlmittelmengen naherungsweise vorliegen, vereinfacht sich das System der Konstanten und es 
wird A, = 0, Bb, = B,. Der Verlauf der Temperatur in z-Richtung wird dann ausschlieBlich durch 


gerade Funktionen symmetrisch zu z = o dargestellt. 
Eine besonders wichtige Rolle spielt im Elektromaschinenbau die mittlere Wicklungserwar- 


mung. Sie wird mit den Werten der Integrationskonstanten (4,12) und den Gln. (4,8) und (4,11) als 
Integral-Mittelwert bestimmt aus 


—2; + 4 2s 
Vin | Jf 9,4(2) dz + J Oyl(2) dz + f,2(2) | , (4,13) 


woraus sich ergibt 


UN E ! oy U, e = F(3 X'y 2) 4: i! ee F| + ¢% 
a a Finet Lo» 


‘a ae = 
ES Oe Oral ize oO. Ue \ Be Oo, U; 
an 13a 
| is whe : (4,13) 
tgh ay 2, a, 


1 ae 


an 2; n 
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In Gl. (4,13a) sind die von #, abhangigen Funktionsanteile verschwunden. Sie leisten keinen 
Beitrag zum Mittelwert. Vereinigt man die Einfliisse der Wickelkopfkithlung in den Abkiirzungen 


I 
ay Uy a5 a, Us s a, U, a 4s 


so laBt sich die mittlere Wicklungserwarmung, bezogen auf Kiihlmittel-Eintrittstemperatur, an- 
geben zu 


SL enc (4,15) 


0 = 
E l 
m aun a 


Hierbei ist eine Gliederung nach Art der fiir die einzelnen Temperaturanteile verantwortlichen 
Verluste vorgenommen. Um einen Anhalt tiber die GréBe des Wickelkopfeinflusses zu vermitteln, 
sind in den Bildern 12 und 13 die in den GIn. (4,14) vorkommenden Faktoren M und N fiir einige 
Verhaltnisse z,/z, tiber dem Quotient a,/a, aufgetragen. Es wurde dabei ay z, < 1 vorausgesetzt, 
so da} die Naherungen gelten: 


tgh ay 2, & ay 2, tgh a, (2, — 2) © 4, (% — 41): 


Fir gréBere Argumente der Hyperbelfunktionen verlaufen die M-Kurven im Bereich ay/a,> 1 
steiler, im Bereich a,/a,< 1 flacher als die dargestellten Kurven. Die Abweichungen sind fiir 
dy % = 1 und z,/z, = 2 zwischen den Grenzen 0,5 und 2,0 fiir ay/a, nicht gréBer als 1096. Der 
groBte Wert der Abweichungen entsteht fiir die kleinen Werte von a,/a,. Auch fiir a, z, = 2 kann 
mit den genannten Parametern von z,/z, die Darstellung des Bildes 12 als erste Naherung heran- 
gezogen werden. 


2,0 


M—e 
mas 
=) 
ax 
r=) 


ay 


-__ 
as 


Bild 12. EinfluB der Wickelkopfkiihlung auf den Vy proportionalen Anteil der mittleren Wicklungserwarmung ; 
Erlauterung siehe Text. 


Ahnliches gilt auch fiir den Faktor N, der den Einflu8 des Wickelkopfes auf die von den Eisen- 
verlusten abhangigen Erwarmungsanteile angibt. SchlieBlich sei noch darauf hingewiesen, daB 
fiir eine grobe Naherung zur Erfassung des Wickelkopfeinflusses der Faktor F = 0 gesetzt wer- 
den darf, wenn ay/a, nur wenig von 1 abweicht, was in vielen Anwendungsfallen gegeben ist. 
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Bild 13. Einflu8 der Wickelkopfkiihlung auf den durch die Eisenverluste verursachten Anteil der mittleren Wicklungserwaérmung; 
Erlauterung siehe Text. 


Zusammenfassung 


Es werden Gleichungen zur Berechnung des Temperaturfeldes von indirekt axial gekihlten 
elektrischen Maschinen angegeben. Die naherungsweise Beriicksichtigung der zweidimensionalen 
Warmestromung innerhalb des Kupferquerschnitts der Wicklung erlaubt auch die bei breiten 
Nuten notwendige Erfassung des iiber den Nutengrund str6menden Warmeflusses fiir alle Rand- — 
bedingungen der Warmeabgabe. Eingeschlossen sind dabei die beiden Sonderfalle einer ein- 
seitigen Warmeabgabe. Der EinfluB des Laufers auf die Erwarmung der Standerwicklung 
wurde durch Einfiithrung der iiber den Stander abzufiihrenden Lauferverluste beriicksichtigt. 
Fiir Einschicht- und Zweischichtwicklungen ergeben sich bei Beriicksichtigung der radialen 
Warmestromung in der Nut etwas unterschiedliche Ergebnisse. Sind die Warmeabgabezahlen 
fiir die Wicklung innerhalb des Blechpaketes und im Wickelkopf nicht gleich, so entsteht 
auch in axialer Richtung ein Warmestrom zu den kalteren Teilen der Wicklung hin. Die Beriick- 
sichtigung der Wickelkopfkithlung kann im Ergebnis als Korrekturfaktor fiir die ohne axiale 
Warmestr6mung gewonnenen Ergebnisse angesehen werden; besonders im Hinblick auf die mitt- © 
lere Wicklungserwarmung laBt sich ein gut ttberschaubares Ergebnis ableiten. 

Der EinfluB der Gaserwarmung wird durch die Annahme erfaBt, daB der auf der Gaseintritts- 
seite liegende Wickelkopf ausschlieBlich durch Gas mit Eintrittstemperatur und der auf der 
entgegengesetzten Seite liegende Wickelkopf durch Kiihlgas mit Austrittstemperatur gekiihlt 
wird. Fir die Kiihlung des Eisenteils wird eine lineare Temperaturzunahme in axialer Richtung 
zugrunde gelegt. Die Ergebnisse der Rechnung sind durch Ubergang auf normierte Gleichungen, 
die Durchrechnung einzelner Zahlenbeispiele, sowie durch Hinweise auf Sonderfalle fiir die An- 
wendung vorbereitet. 


Dr.-Ing. H. Wren, Rheydt, ReitbahnstraBe 56 
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Der Erwarmungsvorgang in kreisférmigen dielektrischen Scheiben 
im Hochfrequenzfeld von Hohlleitern mit H,,-Wellen 


Von 
H. BucHHo1z, Darmstadt 


Mit 1 Textabbildung 


(Eingegangen am 28. Mat 1960) 


Ubersicht: In der vorliegenden Arbeit wird der Erwarmungsvorgang im Inneren einer warme- und 
stromleitenden kreisférmigen Scheibe in Form eines geraden Kreiszylinders berechnet, die in einem Hohlleiter 
angeordnet ist. Die sie beheizenden Warmequellen mégen ihrer Grée und 6rtlichen Verteilung nach dem hoch- 
frequenten, axialsymmetrischen, elektrischen Feld einer H,,-Welle entsprechen. 

Die Aufgabe wird mit Hilfe der GREENschen Funktion der Warmeleitung und unter Einsatz der Methoden 
der Laplace-Transformation gelést. Es werden sowohl der Erwarmungsvorgang selbst als auch der stationare 
Endzustand angegeben. Von der die Scheibe radial begrenzenden Zylinderflache wird angenommen, daB sie 
stets auf Raumtemperatur gehalten wird. An den beiden ebenen, kreisf6rmigen Grund- und Deckflachen der 
Scheibe wird ein Warmeausgleich mit dem AuBenraum gemaB der dritten Randwertbedingung angenommen. 

Bei den Angaben fiir die Endlésung wird noch besonders auf den Fall kleiner Erwarmungszeiten eingegangen. 
Dagegen werden alle Materialkonstanten, soweit sie benétigt werden, als temperaturunabhangig angesehen. 

Kreisf6rmige Hohlleiter, die Richtfunkantennen speisen, mitissen gegen das Eindringen von Feuchtigkeit 
in das Innere des Hohlleiters geschiitzt werden, da sie zumeist auf hochgelegenen Dachern von freistehenden 
Gebauden aufgestellt werden und dort den Einfliissen der atmospharischen Luft ausgesetzt sind. Der erforder- 
liche Schutz gegen die Unbilden der Witterung wird in der Regel in der Weise vorgenommen, daB man in der 
Nahe des Endquerschnitts in den Hohlleiter eine durch einen diinnen Rahmen eingefabte Scheibe von einigen 
Zentimetern Dicke einsetzt, die aus einem dielektrischen und zumeist porésen Material besteht. Bei dieser 
Konstruktion ist es in der Praxis haufig vorgekommen, daB schon kurze Zeit nach dem Einschalten die Ein- 
satze in Brand gerieten. Die physikalische Erklarung dafiir hegt sehr nahe. Es wirken dabei die folgenden Ur- 
sachen zusammen: Die dielektrischen Verluste im Scheibenmaterial, die natiirlich eine, wenn auch nur geringe 
Erwarmung hervorrufen, die sehr ungleichmaBige Verteilung der Stromdichte tiber die Scheibenflache und die 
gerade bei porésem Material auBerst schlechte Warmeleitfahigkeit. 

Da der Erwarmungsvorgang im Inneren verlustbehafteter Dielektrika im Felde hochfrequenter elektro 
magnetischer Wellen allgemeines Interesse hat, soll ein solcher Vorgang hier fiir den Fall der Hy,-Welle genaue 
durchgerechnet werden. Es wird dabei die Methode der Greenschen Funktion der Warmeleitung angewendet 


1. Die Laplace-Transformierte der Temperaturverteilung in einer kreisformigen, 
dielektrischen Scheibe unter dem Einflu8B der durch Ho,-Wellen 
verursachten Warmequellen 


Wir betrachten eine aus dielektrischem Stoff gefertigte, kreisf6rmige Scheibe mit dem 
Radius a und der Scheibendicke 2 d. Zur Fixierung der einzelnen Punkte in der Scheibe benutzen 
wir selbstverstandlich ein Zylinderkoordinatensystem (0, gy, 2) und ordnen es in der Weise zu der 
Scheibe an, daB deren Grund- und Deckflache in die Ebenen z = —d und z= + d zu liegen 
kommen und der Scheibenradius 9 =a ist. Zum Zeitpunkt ¢=0 moge iiberall in der 
Scheibe #(x, y, z; 0) gleich Null sein. 

Um den WarmefluB in der Scheibe zu beschreiben, gebrauchen wir- die folgenden Stoff- 
konstanten: 

Die Warmeleitfahigkeit A in W/(cm, grd), 

die spezifische Warme der Volumeneinheit c in Ws/(cm’, grd), 

die Temperaturleitfahigkeit y? = A/c in cm?/s, 

die Warmeiibergangszahl H in W/(cm®, grd), 

die relative Warmeiibergangszahl h = H/A in 1/cm, 

die elektrische Leitfahigkeit x in S/cm; 

die Dielektrizitatskonstante e, und ¢ der Luft und der Scheibe in F/cm mit é/é) © 1; 

die Permeabilitaét uw) der Luft und der Scheibe in H/cm. 


Alle diese Konstanten werden als temperaturunabhangig angesehen. 
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Im Zeitpunkt ¢ = ¢’ > o sei nun im Punkt (x’, y’, 2’) wahrend der unendlich kleinen Zeit dé’ 
eine Warmequelle wirksam, die in dieser unendlich kurzen Zeit die Warmemenge Q’ - dé’ in Ws 
freisetzt. Dann breitet sich in der Zeit ¢ > t’ diese Warmemenge in dem zunachst als unbegrenzt 
angesehenen Raum nach allen Richtungen gleichmaBig aus, und solange sich die Material- 
erfiillung des Raumes nicht andert, ist im Punkt (x, y, 2) zur Zeit t > ¢’die Temperaturverteilung 
durch den Ausdruck 

Q’ © dt’ —yr 2 , 
: sage’ ay’ pte Pf ye ; 7] (472+ (t—?#/) 0 AEN) 
CON, Ve hae ose ) oye € (Ord?) == Ws () 
Pee AY tthe ee 


gegeben. Da nun bekanntlich 


+ co 

[Cee dae Nar iayne (jarc a| = a2) (2a) 
ist, so folgt aus Gl. (1), wenn man nach allen drei Koordinaten x, y, z zwischen den Grenzen 
—oo--:-+ oo integriert, daB die gesamte, wahrend des Zeitelementes dt’ freigesetzte Warme- 
menge mit dem Ausdruck 

+00 +00 +00 
é- ff .f dO,@iv.2549,254—#) da dy dz 08 dr GS 7) 2 Ws (2) 

iibereinstimmt. Damit ist die oben angegebene physikalische Bedeutung der GréBe Q’: dé’ in 
Gl. (1) bestatigt. Die Funktion #(x, y, z) von Gl. (1) gentigt im tibrigen, wie man sich leicht tiber- 
zeugen kann, auBerhalb der Stelle (x’, y’, 2’) tiberall der partiellen Differentialgleichung 


yi - AX = Odd, (3) 
worin 
Leck ae 
Bre at os (3a) 
oder besser 
a One 


602 0 6 @ 02” 


ist, falls etwa das Temperaturfeld auBerdem axiale Symmetrie aufweisen sollte. 

Natiirlich hangt die Frage, ob im vorliegenden Falle in der Scheibe das Temperaturfeld 
tatsachlich eine axiale Symmetrie hat, nicht allein von der Geometrie der Scheibe ab. Auch bei 
Kreisscheibenform und iiberall gleich groBer Scheibendicke ist z. B. dazu iiberdies noch erforder- 
lich, daB auch die Warmequellen im Inneren der Scheibe so verteilt sein miissen, daB sie eine 
solche axialsymmetrische Temperaturverteilung auch wirklich erzeugen. 

Das ist nun in der Tat der Fall, wenn die zur Scheibenebene senkrechte Scheibenachse zur 
z-Ebene des Zylinderkoordinatensystems gemacht wird und diese Scheibenachse und also auch die 
z-Achse zugleich die Achse des Hohlleiters mit der Hy,-Welle darstellt, denn fiir das elektro- 
magnetische Feld der H),-Welle bestehen im eingeschwungenen Zustand dieses Feldes die folgen- 
den drei Gleichungen: 


ny ceed PRO ia, elice an) ieee 
E, (@, 2; 2) =—i-fi,,-\/ 2. T5(Lie,| = 4 tT gy 2 ae Fs (24,), (4a) 


€ 


(0, 2;t) =: Hy,- Berl7a)*—1- Ja (2 Hey SiGe) =O (= Aaa ae (4b) 
nites a a Q@ er ry 
H, (0, 2; t) = H,, Jol ic,) [Zo,] = A/cm, (4c) 
Die hierin noch nicht angegebene Abhangigkeit dieser drei Feldkomponenten von z und ¢ ist 
durch den fiir alle drei Gleichungen gleichartigen Faktor 


a wt wth i-2nx/A@ M3 i\ 
Hy, = Ho, - omit A= ° Gey (4) 
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gegeben, wenn es sich wie im Falle des oberen (unteren) Vorzeichens um eine in Richtung der 
positiven (negativen) z-Achse wandernde Welle handelt. In den Gl. (4), (4a, b, c) bedeutet 

H,, die Amplitude des v-ten Wellenmode der H,,-Welle mit der Dimension A/cm, 

Volo © 120% den Feldwellenwiderstand der freien Raumwelle in Ohm, 


a der Radius des Hohlleiters in cm, der gleich dem Radius der dielektrischen Scheibe ange- 
nommen wird, 

Jor =i, Hie r-te Nullstelle der Funktion Jj (x) = — Jy(x) in bezug auf x, 

7, die freie Raumwellenlange in cm mit der dazugehérigen Wellenzahl ky = 2 2/A, 

Ag, die Grenzwellenlange des 7-ten Wellenmodes der H,,-Welle in cm, wenn der Hohlleiter 
verlustlos ist, und schlieBlich ist 


A) die Lange der Hohlleiterwelle selbst. 


Statt der GréBe Hy, mit der Dimension A/cm, die in den Gl. (4a, b, c) benutzt worden ist, um 
die Intensitat des magnetischen Feldes der H,,-Welle zu kennzeichnen, wird in der Praxis die GroBe 
des Energiestroms der betreffenden Hohlleiterwelle verwendet. Um diesem Gebrauch Rechnung 
zu tragen, wird an dieser Stelle, um jederzeit Umrechnungen zu erméglichen, der Zusammen- 
hang zwischen dem Wert von Hy, und dem zu dieser Teilwelle gehérenden Energiestrom S“ 


angegeben. Er lautet wegen! Hoy Ht, 2 1bh 


~ ly —— dee AeeNe a Se 
SP) = = a? Vugleo ' le hos i, : y( on) 1 Jollee) [S.J =W. (40) 


In der wirklichen, eingangs beschriebenen Ausfiihrung eines gegen das Eindringen von 
Feuchtigkeit geschiitzten Hohlleiters trifft nun allerdings die aus dem Hohlleiter kommende 
Welle in der Scheibe selbst auf ein Material, das im allgemeinen eine andere Dielektrizitats- 
konstante ¢ haben wird als die Luft im Hohlleiter. Da man jedoch ohnehin darauf bedacht sein 
muB8, diese Verschiedenheit in den e-Werten so klein wie méglich zu halten, um Leistungsverluste 
bei den Reflexionen zu vermeiden, mége hier von solchen Reflexionen vor und in der dielektri- 
schen Scheibe abgesehen werden. Die Stromdichte in dem verlustbehafteten Dielektrikum selbst 
1aBt sich dann nach Gl. (4a) mit v6llig ausreichender Annaherung unmittelbar durch den Aus- 
druck 


~ 1/2 , 
inl, 231) =i. He, (22) oe (2%, Alem? (5a) 


1) 


wiedergeben, und somit wird dann durch diese Stromdichte im zeitlichen Mittel in der Zeit- 
und Volumeneinheit die Warmemenge 


S y Os 1 , ,, ’ , Da 1 A : I , Of? 
dle) = Fo = SB yle's 250-180 25) = 2H, worden (A) -J59(L ie.) (8d) 


erzeugt. Die Warmequellen haben also in der Tat eine axialsymmetrische Verteilung, die iiberdies 
nur vong und nicht auch noch von z abhangt. Wiirde fiir das Scheibenmaterial ¢ > é sein, 
so beséBe allerdings die Warmequellenverteilung eine Abhangigkeit von z, aber sie ware auch 
dann noch axialsymmetrisch. Anstatt weiterhin von einer punktformigen Warmequelle zu 
sprechen, sind wir durch die Gl. (5b) auch bereits in die Lage versetzt, den Einflu8 einer kon- 
tinuierlichen Quellenverteilung auf die Erwarmung zu beriicksichtigen. Aus methodischen 
Griinden wollen wir jedoch vorerst noch bei der punktformigen Warmequelle bleiben. 

Von Gl. (1) aus gelangt man dann sofort zu dem Ausdruck, der einer zeitlich ununterbrochenen 
Wirkung der punktférmigen Warmequelle am Orte (0’, y’, 2’) vom Zeitpunkt ¢’ = o an bis zum 


1 Das hochgestellte Zeichen * bezeichnet in dieser Arbeit, daB der konjugiert komplexe Wert der betr. GroBe 
zu nehmen ist, die Uberstreichung eines Funktionszeichens, da es sich um die Laplacetransformierte der betr. 
Funktion handelt. : 
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eT 


jeweiligen Zeitpunkt ¢’ = t entspricht. Es bedarf dazu nur der nachstehend angegebenen Inte- 


gration nach der Zeit ¢’ zwischen den Grenzen 0... ?: 
f t 
dp - oo! Panel (7s) aie Pere kita CR dt’ ey 00 WEooteen 
C> WO Qe a Cis ole a= 8 (x y?+ (¢—?#’))2? Pp. S S 
6 
t 
_-@ sey 2 : ds a 
az a : fexe [ Y (4 Tt) | (x y2 7) 3/? (7 ) 


0 
(7 = 0 + 0? — 200' cosg). 


Die linksstehende, mit c multiplizierte Temperatur 9, ist wegen der endlichen Einwirkungsdauer 
der Warmequelle nun auch von derselben GréBenordnung wie Q’. Mittels der Formel 


fen( fw re) di= = for ; | fi 7 
0 


0 


lee} 


) @>0) &) 


ergeben sich dann fiir die LaPLAcE-Transformierte von @, nacheinander die folgenden Ausdriicke: 


lee) 


C. ea O20, P, 23.01, 9, 25 t) + dt = OD(0, @, 250°, 9, 2; p) -€ 
0 


yo ae : ck ; di 
= | eF. exp {— [(z—z’)? + @ +e *—2e@ 08 PNG 2) eae 
_ _ O __— exp{ le? + 0? — 200 cose + (2 —2')1" - Vey 
4p ny [o? + 9° — 29.0" bcos of 4 (ea 
ar x [ cosce—2) Ky (VS? + ply? -o? + 0? — 200" -cosp)- do, (9) 
0 


denn es ist nach einer bekannten Formel aus der Theorie der LAPLACE-Transformation 


— — #7/(4 0 dt —x- Vpla 
fe Pee Bile Ds part = Alera - 2) -¢ x Vpla (10) 
0 
und zwar auf Grund der allgemeinen Beziehung 
1/2 - (2s def Ke) 9 (lance ar nobel (10a) 
0 


Die Gré8e p in Gl. (9) und in den folgenden Gleichungen ist vorlaufig stets als positiv reell 
anzusehen. 


Mit der letzten Zeile in Gl. (9) ist eine Darstellungsform erreicht worden, die mit Leichtigkeit 
die Eigenschaft der axialen Symmetrie der Quellenverteilung zu beriicksichtigen gestattet, denn 
es ist bekanntlich 


7 ; + 0 : K, (Aa D(a , = , 
Roiree—rean~ $f on fhbar tio Goa 


n=— oOo 


Nun ist nach Gl. (5b) in der letzten Zeile von Gl. (9) 


my 


Q’ = 0 dq’ f do’ F dz’ F q,(0’) : 
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Integriert man nunmehr iiber ’ zwischen den Grenzen o...2a und bringt auch noch an der 


Temperatur O» der verstandlicheren Deutung wegen die Orantine des raumlich unendlich Kleinen 
zum Pedearle so erhalt man: 


c-O,(0, 2:0, 7; p) = 2$ 87,0). | cos f(z — 2’) 


. eee + ply?-0)- Ive + ply a pen Cee fea 


K(/G? + ply? 0’) - Tole? + diy? @) (0 <@’) 


In der Gl. (12) haben wir die LApLAcE-Transformierte der Temperaturverteilung in einem 
unendlich ausgedehnten, warmeleitenden K6rper vor uns, der seit dem Zeitpunkt ¢ = o von einer 
Warmequelle der Ergiebigkeit 7,(0’)W/cm® beheizt wird, die langs eines Ringes mit dem Radius 0’ 
in der Ebene z = # in zeitlich immer gleichbleibender Starke wirksam ist. Die partielle Differen- 
tialgleichung, die die Funktion @,(0,z; 0’, 2’; p) von Gl. (12) fiir alle @ und z mit e +0’ und 
z + 2’ erfiillt, hat im Hinblick auf Gl. (3) die Form: 


2Op L 1 2OR 2OR p 


as -On,=0 (O, = #O,(0, 2; 0’, z))- (13) 


Cieelone Oe eye 

Die vollstandige Lésung der vorliegenden Aufgabe hat nun aber nicht nur die Gl. (13) zu 
befriedigen, sondern sie hat auch noch gewisse Vorschriften an den Grenzen des von der Scheibe 
eingenommenen Raums einzuhalten. Als solche Forderungen moégen hier gelten: 

I. Die Bedingung @,,,(a, z) = 0 fir —d 2 z<+dundo=a, 

II. a,b) die Bedingung @,,,(0, 2)/dz7-+ h- 9..(0, £ d) =o fir z=+d undo=oZa. 


Die Bedingung I. 1aBt sich von der Gl. (12) her ohne viel Mihe erfiillen: Es ist ndmlich offenbar 
nur notig, die unter dem Integral stehende geschweifte Klammer sowohl im Falle @ > @’ wie im 
Falle 9 < o’ durch einen und denselben Term 4A (C) - Io ye? + p/y) zu erganzen, denn in der Tat 
wird bei dieser Wahl die Diff.-Gl. (13) so gut erfiillt wie durch den entsprechenden Term mit Ko. 
Die Bedingung I. selbst fiihrt dann sofort zu einer Bestimmungsgleichung fiir die unbekannte 
Funktion A(¢) von €. Man errechnet fiir sie den Ausdruck: 


A(C) = — Ip w Q’) + Kola w)/Lo(a4 @) 
und damit wird jetzt 
LET pos de ; , 
¢- POR (9, 2502; p34) = ote ile’): | costle— 2’) M (wie, e's a)- at 
0 
i po do dz’ = , F Iy(0’ w) ; K I K d 
Se he Ot | COS Cty [Ko(w g) Io(w a) — Iy(w@) Ko(wa)]- do (14) 
LW, Iy(a w) 


(ee = Ce Piya 0-0) 


Fiir 9 < o’ lautet das Integral in der zweiten Zeile von Gl. (14): 


lee) 


yA LE ) 7 / 
*f costte—2 an [K,(w 0’) I)(w a) — I,(w 0’) Ky/w a)}- dl 


(A Cera 0 0 Js 


Fir a = 0 >’ wird in der ersten Form der Gl. (14) in der Tat d?@, = 0, wie es sein mu8. In der 
obigen zweiten Form dieser Gleichung wird aber in ihrem Geltungsbereich g tberhaupt nicht 
gleich a. Setzt man noch in (14) den Faktor 1/2 vor das Integralzeichen, so muB das Integral 
zwischen den Grenzen —oo-::+ + oo genommen werden, wenn es dem friiheren gleichwertig 


0 
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= LL 


sein soll. Der Integrand in Gl. (14) hat in der ¢-Ebene unendlich viele einfache Nullstellen, und 


zwar dort, wo 


aw = al? + p/yr=71-Jo, i: 1) 28 oere) (144) 
ist. Es ist also mit€ = &+17 
al,=ti-[f, tb: Ply] = 1t- 1,4 (n = 1, 2,3..-) (148) 
Alle diese Pole liegen also symmetrisch zur reellen Achse auf der Achse der rein imaginaren 
Zahlen, und selbst im Falle # = 0 ist stets y, > 0 flirm =0,1,2,.... 
U7 
TY ¢+in=C-Ebene 
a3 


Se t fe es en a 
ee ae on ci = 5 ae ne, He #C, +l +¢, iG 


U7. ES 


Bild 1 erlaéutert anschaulich den Verlauf des komplexen Integrationsweges in Gl. 1 (19). 
In der Frage der Konvergenz ist alles Erforderliche im Text angegeben. 


In Gl. (14) 1aBt sich nun offenbar das Integral auch in der Form schreiben: 


+00 > 


, . . 
ets . | [at Ge CM 6.0 * ands 
2 


== 12,9) 


+o00+te +oo—ie 
1 5 / 3 1 Ae —se 
ee dels —#l. M(w; 0,0’; a) - do + —. i et Fl Mw; 0,0 s@)ere 
2 2 
—co+te —oo—ieé 


OL ae ee a 


Der Ubergang zu den absoluten Betragen von z — 2’ erlaubt es dann sogar, die Integrationswege 
in den letzten beiden Integralen etwa symmetrisch zur rein imaginaéren Achse unter dem Winkel B 
mit 0 < $ <2/2 ins Unendliche auslaufen zu lassen. Das erkennt man sofort, wenn man die 
Integranden in Gl. (14) durch ihre asymptotischen Naherungsausdriicke ersetzt. 

Das zweite Integral in der letzten Gleichungszeile laBt sich noch durch den Ersatz von ¢ durch 
—¢ in das erste Integral dieser Zeile iiberfiihren, so daB8 man fiir die Gl. (14) auch die Schreibweise 


os 
Oo+ arc 
2 


c-O2(0,2;0',2 )p) a) = 2 2 5 BG iB p—x| , Zola’) | Tela ®) Sioa 


ba Ko(aw) Kolo w) 


Oe eae Py) pane (14') 
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wahlen kann. Darin zieht der Integrationsweg in der ¢-Ebene unterhalb der kleinsten Nullstelle 
von I,(a w) vorbei. 

Um nun auch die Grenzbedingungen II. a,b) zu erfiillen, erweitern wir den einzigen, von z 
abhangigen Bestandteil des Integranden in der letzten Gleichung zu dem Ausdruck 


eb“ 4. ibe 4 Baits 


und bilden hiermit gem&B der Bedingung IT. b) €0,/éz LhO, fiir z= +d und z = —d. Auf 
diese Weise erhalten wir die beiden folgenden Bestimmungsgleichungen fiir A und B: 


(A eta) ed 1 (Be etia)  otibd — _ pitd, ptibe ta) (15a) 
(A-eti%). et ild 4 (B. eis), g—ihd _ __ git d, gi (Ce —a) | (15b) 

Dabei ist hierin voriibergehend der Hilfswinkel a durch die Gleichungen 
pois, (16a) sin ga, (16b) cos 2a =, (16¢) 

sina = C/V? +2, (16.0) cosa =hfy~h? +? (168) 

definiert. Man errechnet dann aus (15a, b) fiir A und B die Ausdriicke : 

A= —elorey Sina +6 (@a2) (17a) 
sin 2 (a + € d) 

B= — pletta lo ae @ st (17b) 


sin 2 (« + € d) 
und erhalt damit im ganzen: 
21 

sin 2 (a + € d) 

sin (a + ¢(@ + 2’))-sin(a+€(d—2)) (e>2’) 

sin (a + €(d—2z’))-sin(a +€(@+2)) (<2) (18) 
Eine leichte Rechnung ergibt, daB die linke Gleichungsseite in der Tat die Bedingung IIb) erfiillt. 

Eliminiert man wieder aus Gl. (18) die Hilfsgr6Be « mittels der Gl. (16a, b) und (18a, b), um 

danach die Natur der auf der linken Seite von Gl. (18) stehenden Funktion von ¢ beurteilen zu 
k6nnen, und setzt die Glieder zusammen, so ergibt sich als Lésung im ganzen das folgende langs 
des ganzen Integrationsweges absolut konvergente Integzal: 


; 7,(0’) 


eto z—2"| ne Ais. etic: of. B.e ot? = 


4° 9’ do’ dz’ 


C 2?0,,.(0, zy; 0’, gap: a, = d) at 2 
Dies 


co aro (= —A 
1 (€- cos C(d + 2’) + h-sinf(d + 2’)] [€- cos C(d — 2’) + h- sin C(d — z’)] 
(| 2hC-cos2Cd-+ (h®?—C?)-sin2Cd 


I,(aw) Lp(e w) | Ly(we’) de [c- 2O,.:] = Ws?/cm3 (19) 


(== h)y2) 0 0, =. 0,0 S60, a ot = dh ae Se): 


Wird der frithere Hilfswinkel « der Gl. (16) benutzt, so ist im obigen Integral der Nenner 
2hC-cos2¢d-+ (#—C?)-sm2¢d= (W+2*%)-sn2(a+ Cd). (19’) 


Ist in Gl. (19) @ <0’, so ist das Zeichen @ innerhalb der Determinante durch 9’ und das 
Zeichen 09’ auBerhalb der Determinante durch @ zu ersetzen. Ist z > 2’, so sind im Integranden 
die Vorzeichen der Glieder z und z’ in beiden eckigen Klammern des Zahlers miteinander zu ver- 


tauschen. Im iibrigen stellt jetzt die Funktion d?0,,:(0, 2; 0’, 2; p34, af d) die LapLAcE-Trans- 
formierte der Temperaturverteilung in einer dielektrischen Scheibe dar, die von einer ringférmi- 
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gen Warmequelle am Orte (0’, 2’) innerhalb des Bereicheso = @ < a und —d=zS + d her- 
riihrt und die in zeitlich stets gleichbleibender Starke wirksam ist. 

Der Integrand in Gl. (19) ist offenkundig eine in ¢ eindeutige und ungerade Funktion. Geht 
man also in dem Integral von € zu —€ iiber und wahlt als Integrationsweg den in Abb. 1 dar- 
gestellten Weg, der aus einer oberen und unteren Halfte besteht, so braucht man nur auf der 
rechten Seite von Gl. (19) den Faktor 1/2 einzufiigen, und es stellt dann das so veranderte Integral 
wiederum die Funktion @,,, von Gl. (19) dar. 

Die Pole des Integranden von Gl. (19) sind nach wie vor die Stellen€ = + 7, mit Im(n,) = 0 
von Gl. (14). Dahinzu kommen jetzt noch die reellen Nullstellen des Nenners im ersten Bruch 
des Integranden der Gl. (19). Wie man leicht nachweist, kann nun aber geschrieben werden: 
2hC-cos2¢d- (h®?—C*)-sina¢d=—2:[h-sin¢d+C-cos(d]-[h-cosCd—C-sinéd]. (20) 


n 


Es gibt daher zwei unendliche Folgen von reellen Nullstellen. Fiir die Nullstellen € = &2) = h . &) 
des ersten Faktors ist 


tg (6) d) =—€E/h oder tg (EY “h d) = — &) [es 2) == em (20a) 
und fiir die Nullstellen¢ = &@) =h. E°) des zweiten Faktors: 
tg €2) d) = hie - ‘oder “1g (E2) “h@) == 1/E2) 4 (20b) 


Die numerische Berechnung von Et) und Ee) erfolgt mithin zweckmaBig mittels der beiden 
Gleichungen 
bm ete by =—hd <0 (20a), Ey. tee —+hd>o (208) 
(EO?) — £02) Ad ; m= 1;2,3.. e 


Mit der Einfiihrung der Nullstellen Eu? ist erreicht, daB deren Lage alleine von dem dimensions- 
losen Produkt 4d abhangt. Mit + &”) ist auch — &2) fiir jedes m eine Nullstelle. 


m m 
Es ist nun aber nach einer einfachen Rechnung: 


lzhe-cos2td + (@—C)-sin2d}i | = 2/d-[hd- (1 + hd) + 8,3] (21a) 


Cn = oe m=1,2,3...). 


Damit entsteht im ganzen als Resultat, zu dem die Auflésung der Gl. (19) gefiihrt hat, die folgende 
erste Reihendarstellung: 
PO (0, 250,27; p34, d) = 7 8 2 Gy 
c tot(Q 2; 0", 2°; P; a, + a) poe od q,(2’) 
o(@ WF) Io/o 0%)| sim (EP 2) + sin (EW 2’) «sin? (EP a) (HR a2 4 El" gaye 
) 


Nae Ss Iy(Q’ Om) if m m 
har Tot Mm) [Koa vf) Ko(o w) (Ewa? + hd (ahd) eB) ge 
Io(a wh) — Too wh)| cos (€@ 2) + cos (E9 2’) « cos? (E@ d) (2 a® + £2 gay 
K,(aw™) Kolo w®) (Ea)? + hd- (1 + hd) (2) a)? (22) 


™m 


yrs 
| Ty(a w®) 


ry ER)? pe i Bee, Oe) eee en ae een 


In der unendlichen Reihe der Gl. (22) strebt das vorletzte Glied hinter dem Summen- 
zeichen fiir m — co wegen &,, > oo gegen den Grenzwert 1/d?. Es besteht also fiir die in Gl. (22) 
auftretenden Reihen nur so lange Konvergenz, als 0’ <¢ ist. Fir 9’ = 9 und z = 2’ hort jede 
Konvergenz auf. Das liegt ganz im Charakter der GREENschen Funktion, denn als nichts anderes 
ist der Ausdruck (22) aufzufassen. Die physikalische Dimension jeder Gleichungsseite von (22) 
betragt Ws?/cm%. 

Fiir die physikalische Deutung der oben dargestellten Liésung ist nochmals daran zu erinnern, 
daB sie eine Aussage tiber die LapLacg-Transformierte der zeitlichen und értlichen Temperatur- 
verteilung in einem warmeleitenden Kérper macht, dessen innere Warmequellen durchweg auf 
dem Ringe (0’, 2’) liegen und dort nach Ort und Zeit in tibérall gleichmaBiger und gleichbleibender 
Starke ausgebreitet liegen. 
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Was nun die Warmequellen bei der durch Gl. (5b) angegebenen értlichen Verteilung anbe- 
langt, so erlaubt es die fehlende Abhangigkeit von z’ in der Verteilung der Quellen sofort von der 
ringformigen Quellenverteilung zu der anderen, gesammelteren tiberzugehen, die einem zwischen 


den Ebenen z = + d mit Warmequellen besetzten Zylindermantel vom Radius 0’ entspricht. 
Dazu ist nur notig, in der Gl. (22) cos (€@) - 2’) zu ersetzen durch 
+d 
focas (E@- 2') dz’ = 2 h/(E)?. cos (E% d) [E27] =1fem, [hk] = t/om, ~ (22’) 
tr 


denn die nach den Nullstellen &@) fortschreitende Entwicklung in Gl. (22) fallt bei einer Inte- 
gration zwischen den Grenzen — d--- + dheraus. Man hat also jetzt den einfacheren Ausdruck : 


5 eadon re P 
6 d@s. (0 2303p) 4, + Le ee 


x S48 0(0' whe) | Lola wee) Lolo wy))) cos (E52) + cos? (EG d) hed? + Ea) a (22a) 
w= Tlawe |Ky(aw®) Kew) Eparthaathd Eas 
Ce Get Pi), w= 4, 2,3..:7-0= 00 =a, dees 1d) OS op <6 = 4), 


Die Konvergenzeigenschaften der in Gl. (22) auftretenden Reihe sind durch den Ersatz von 
cos (&2) z’) durch den Ausdruck (22 ‘) wesentlich besser geworden. Sie entsprechen neuerdings 
denen der Reihe 3’(1/m?). 

Wollte man jetzt, wie es nahe liegt, noch daran gehen, in der Gl. (22a) auch die in der Scheibe 
auftretende Quellenverteilung in Richtung 9 gem48 Gl. (5b) zu beriicksichtigen, so stieBe man 
dabei auf ein Integral, das nicht nur nicht in endlicher Form integrierbar ist, wie das auch auf 
dem anderen, sogleich em hlaecaden Wege der Fall ist, sondern das auch wegen des innerhalb 
des Integrationsintervalles 0 = e’ <a notwendigen Wechsels von der Funktion I)(0’ - w?)) 
zu der Funktion K,(0’ - w')) nur schwer behandelt werden kann. Wir bleiben aus diesem Grunde 
die zweite Reihendarstellung, die das Gegenstiick zu der ersten in Gl. (22) gebrachten Dar- 
stellung dieser Art ware, schuldig. 

Aus der Integraldarstellung (19) 1aBt sich nun aber noch eine zweite, von den Gln. (22, 23) 
verschiedene Lésung herleiten, und zwar dadurch, da8 von dem doppelten Integrationsweg in 
Bild 1 der erste unbegrenzt nach oben und der zweite unbegrenzt nach unten gezogen wird. 
Das ist ohne Wertanderung moglich, weil jeweils die beiden Integranden in den genannten beiden 
Halbebenen mit ihren Absolutbetragen dabei gegen Null streben. Beide Wege iiberschneiden 
aber bei dieser Behandlungsweise des Integrals die unendlich vielen, einfachen Nullstellen der 
Funktion J)(a w). Sie liegen in der ¢-Ebene dort, wo 


aw, =a-~O+ ply? = + i jo, (p > 0) (23) 


ist. Lost man nach der Variablen ¢ auf, so gilt fiir diese Nullstellen in der ¢-Ebene die Formel: 


at, == Be ty a ele aes ope ay? =+1 jus “ope (“4 ) (n= 1,2,3.. .) . (23a) 
On 
Bei reellen Werten von # liegen also alle diese Nullstellen und damit die ihnen entsprechenden 
Pole auf der rein imaginaren Achse der €-Ebene. Man berechnet dann leicht die folgenden, von 
den Nullstellen abhangenden Werte: 


/ 


Loltty-@') = Jolion e), (23) 


a 


Ty(ey @) = JolJon2) (230) 


Kyo W,) = Ky (+ a: n “jue aa 
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Mithin ist im besonderen 


Ky(a w,) = —= 3 Yo(Fon) a Alon Jo Foa)) ? (23e) 
dI,(a w) hol semen ae ot vinta : F 
(er) om EE Silla) \): +e (2), (238) 


Die sich dann auf Grund des Residuensatzes in diesem Falle aus Gl. (19) ergebende Reihen- 
entwicklung nimmt darnach die Form an: 
oo Iolo 2) Jo (‘rr : 2) 
a a 


Jon TilGon) 


Soe dota © 
pry-a q,(0") 


€- d*O). (0,25 0, © spss Gage) 


n=1 


d a fG Ce spi aff Ch ye 
EE * Coj (< Xn (: +3) +hd- gin( Sn ( +=) 
x 
2 Tq: d- Coj (= in) + oy Gin (= in) 
) a a a 
E In: Co} (<1 . (: -4)} +hd-Gin (<a: -3)] 
a a d a d 
oe eee J (24) 
d+ Sin (= in) SF IE wn) : Coj - tn) } 
a a a 


1/2 


. [ : =en7 = , 
(font —Fox"|t PP) Sard Se Ee SS [h] = 1/cm; 0< 0,0 <a). 


Sie stellt die LAPLACE-Transformierte der Temperaturverteilung in einer warmeleitenden 
kreisf6rmigen Scheibe dar, die in Richtung der z-Achse die Dicke 2d und in Richtung wachsender 
Werte von @ den Radius a hat. In der Flache 9 = a mit —d =< z < + d ist die Temperatur 


konstant gleich der Umgebung. In den beiden Flachen z = + d befolgt die Temperaturabgabe 
zu allen Zeiten ¢ > o und fiir alle o = @ S a das Gesetz: 
COrt(0, 2/02 Eh: Opt(o,2)=0 (:=+d, h=H/A). (II. a, b) 


Die Warmequellen in der Scheibe wirken von ¢ = 0 an in immer derselben Starke langs eines 
Ringes mit dem Radius 9’ in der Ebene z = 2’, wobei 0 < @’ Sa und —d=2’ = +d sein 
kann. Trotz dem Wurzelausdruck in den Gl. (22) und (24) ist 0, eine eindeutige Funktion 


von ?p. 
Die zweite Zeile der Gl. (24) vereinfacht sich bedeutend, wenn man statt der eben angegebenen 
Bedingung auch in den Flachen z = + d ein Verhalten gemaB der Relation 


O.+(0, + 4) = 0 fir OS pi a (I.) 


vorschreibt. Es ist fiir die Umrechnung durchaus nicht notig, die obigen Rechnungen fiir diesen 
Fall erneut durchzufiihren. Vielmehr lasse man einfach in der obigen Bedingungsgleichung 
(II. a, b) tiberall 4 > oo streben. Die letzte Zeile vereinfacht sich dadurch zu dem Ausdruck: 


Gin g pase - Sin a aes 
S ag s 2 Gn Gk 


eu d 
lon NIN | Gen 
a 


In der Gl. (24) ist aber nun noch zu beriicksichtigen, da in dem wirklich vorliegenden Fall 
einmal die Warmequellen in Richtung z zwischen z = —d und z = + d gleichmaBig verteilt 
sein sollten. Wir haben also in der zweiten Zeile von Gl. (24) in der angeschriebenen Form, die 
fir —d<= 2 <z< +d gilt, nach 2 zwischen den Grenzen — d--- z und in der Form, die fiir 
—d=az<2/< +d gilt — sie geht aus der angeschriebenen Form durch bloBes Vertauschen 
von z und 2’ hervor — nach z’ zwischen den Grenzen z-++-++ d zu integrieren. Auf diese Weise 
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errechnet sich fiir den Zahler in dem Bruch der zweiten Zeile von Gl. (24) zunachst die Dar- 
stellung: 


ality |2[a - Co} (<r) + <4 Sin (Sy,)| [hd Sin (“ ye) + Ss Col (Sa 


d z 8 d paces d 
—2hd- Go} (a= : =| -lhd- Sin (4) + ue) - Co} (= %) 
Der vollstandige, neue Ausdruck fiir die Gl. (24) lautet daher jetzt: 
C+ PO (0,2; 0, 2; p; a, + d) = 


t ( 4 
8 @’+ do’ co Jolioa =): Jolion =) 1 | hd- oi(= xs) ' ; 
=—— 4, ae : ; x 25 
poms qe) a Fen acon tT es hd: Goi( (Stn) + zm Sin (Sz) 
a a 
i 2 
(tn = ion to = ion V3 | ae = ) ; ~—€@S2544;020,05+4;020; iha|=0). 
Jon’ ¥ 


Um endlich die bisher noch immer nicht beriicksichtigte, raumliche Verteilung der Warmequellen 
in Richtung zunehmender Werte von o oder 0’ zu erfassen, setzen wir die Gl. (5b) in Gl. (25) ein 
und integrieren tiber 9’ zwischen den Grenzen 0--- a. Dann bekommt man, wenn voriibergehend 
die zweite Zeile in der obigen Gleichung mit /(z; 2) bezeichnet wird, den folgenden Ausdruck: 


x eee a hy Gay I) = UA 
¢: O,(0,2; p34, + d) = MH, -(F) - Si Oa fz, n 
tot(2 p ) €0 oe Ag | a ae Jon: JiVGon) Son Jar) K ) 


und hierin ist RECN 


(Jom Sir) = J Jol Jon’ *) Jildor x) + dx (x= @'/a), (%,7=1,2,3...)+ (26a) 


Machen wir hier von der Formel 
+a2/2 
1 or 
JaGor': 3) = Jie) == +. f{ Jo(2 jor» ¥ + CoS p) + cos 2p - dp (27) 
—an/2 


Gebrauch, so erhalt man die Darstellung durch das Doppelintegral 


j COS 2 Y (ic %) Jo(2 Jor» % COS p) + x +) ‘dy , (26 b) 


—71/(2 


alr 


Bon jr) a 


in dem sich die innere Integration nach x mittels einer bekannten Formel aus der Theorie der 
Zylinderfunktionen sofort ausfiihren laBt. Dadurch entsteht schlieBlich der auch noch fiir 
240, CoSpP = +7, in bezug auf — stetige Ausdruck: 

1 


‘ Tilon +) Jol ide x cos ) «x- dx 
0 
Bs? jor COS * Jolion) * Jo(2 Jor * COS P) —Jon* Jo(Jon) Jo(2 Jor * COS @) 


Jon is 4507 K cos? p 
Jo(2 jor * cos —) 


See Te fe \ c- 8 
ea sion) a ee (28) 
Folglich ist jetzt nach dem Ubergang von y zu a/2 — p wegen jj, = jy 
+n/2 
BOF NN ‘ NY Jo(2J1y COS G) ; 
Eons Ia2) = Jon|% Ji(Jon) | jen —(25r, «cos p)? COS 2 de 
Ye aay ) 
— Jon * Sion) : Jo(2 Fiz * COSY : Sond 
ee eee 20"! 
fei ) 
a. fonds JiGon) / | SO\2 Jig ° COS P cos 2m: do. 2 
20-27%, Giese esa" (29) 


=e 
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Der Nenner des letzten Ausdrucks scheint eine Entwicklung nach TSCHEBYSCHEFF-Polynomen 
zu empfehlen. Wir setzen jedoch ¢ = exp (ig) und werden dadurch zu der gleichwertigen, 
aber fiir die weiteren Rechnungen zweckmaBigeren Darstellung 


(0 +) fon +4) 
Jon * J1Von) ; p : Ae : Bape - dt (30) 


Ridon: Jur) Ss 2 ACR eae 
471 jit At Vega on) t 


gefiihrt. Hierin ist die Integration zunachst tiber den Umfang des Einheitskreises der t-Ebene 
auszufiihren, auf dem ja bei dem hier auftretenden Integranden keine Pole liegen. Eben deswegen 
kann aber dieser Weg ohne weiteres auch in das Innere des Kreises hineingezogen werden, so- 
lange er dabei auf keine Singularitaten stoBt. Nun sind in bezug auf die Nullstellen des Nenners 
im Integranden der Gl. (29): 


y= = Pur = = arc COS (ont2diup) mit oO < Pyar < 7/2, 
Po = — Prnr> 3 = Pur, Ps = %—Pur 


In der t-Ebene liegen also alle vier Nullstellen des ersten Nenners im Integranden der Gl. (30), 
ndmlich die vier Punkte 


+tPny ge ont Fe. pre? —tPn4 


toe 9 i 5 ts =— , = 


unter der oben angegebenen Voraussetzung in der Tat auf dem Umfang des Einheitskreises. Das 
in Gl. (30) vorkommende Integral mége mittels der Formel 


J cE 2v a A v el Mah, 2 
Jo(+( +3) =F Ve PO =Rat+ yore +ey re By) 
ausgewertet werden, eine Beziehung, die die wesentlich singulare Natur der Stelle t = 0 des 
Integranden von Gl. (30) deutlich erkennbar macht. Mit ¢ = e'* oder t = e*'(/?—®) folgt aus 
Gl. (31) sofort: 


Ji2% = Sin 0) Jae ) + 2+ am cos (2a) - J2(x), (31a) 
] (2a 60S a) =f) es Bre COS (2'4 Py 2 (4): (31b) 


Andererseits laBt sich nun aber schreiben: 


1 Stl 1 = tn, 1 1 
TEt+CH, i)  F, Phew —F—alq—a q—al | A) 
Entwickelt man hierin den Inhalt der eckigen Klammer in der Umgebung des Nullpunktes der 
t-Ebene nach Potenzen von #, so erhalt man nach kurzer Rechnung den Ausdruck: 


1 Z Sy SIN 2 Yy,(1 + S) 


oe Spada ve CIP hee) © G88 =O SiN 2 Py, 


(i < 1, COS, = a 1). (33) 


Damit laBt sich dann auch fiir die Gl. (30) schreiben: 


Hioisdiy tin alee SS eee 


25iy eR = 2049 
(0+) 
2(v +s) 2(y-+s-+ 2) 
x ae (t + # J. de/t (34) 


(4,7) = 1, 2,32.., Joall@ gap) = a) s 
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Hierin verschwindet die zweite Zeile nur dann nicht, wenn vy = — s oder vy = — s — 2 ist, und 


da s nur gleich Null oder gleich positiven ganzen Zahlen werden kann, so ergibt sich daraus 
sofort die folgende Darstellung: 


Elion Fae) = — 2 SN Uli) + Phat 
— fon tion) 1 726,) a ee Thalia 
+5 meee Balin 
— fon Hee) | Tei.) — 20g 20er S)sin 2 Parle +3) - ("Ste aline)| 
oder in anderer Schreibweise: 
Ejow Jur) = fal {ete 2 uy? 3 (sin 29 Om, «3d —S ald dh 
= Iildon) Sow dae) (C08 Gar = Fonl(2 hi.) <1). (35) 


Bei groBeren Werten des Parameters v ist natiirlich die zu dieser Gleichung angegebene Be- 
dingung nicht mehr erfiillt. Ist aber in der Tat 7,,,/27,, > 1, so sind die vier Nullstellen des 
Nenners im Integranden von Gl. (30) 


+n —Vnr Tena 
=e , p= ; ts = — € ‘ i= 


mit dem unter diesen Umstanden reellen Wert: 


Vir = he oj Gonl2 I.) = = in( Zon -- eS) = 7 s 


2I14 


Die vier Nullstellen sind also jetzt alle reellwertig, und von ihnen liegen nur noch die beiden 
Nullstellen ¢, und ¢, oder + /, innerhalb des Integrationsweges von Gl. (30), da dieser ja mit dem 
Einheitskreis identisch ist. Man rechnet aber leicht aus, daB 


Re} Jo(jiré + 1/4) (+ #)/E — fej [Jo(j1/(¢ + 1/¢ ( TN, ee st =} 
Lf +P Q— Hi, lraty B—# \@—P” RP) et frase, 
= — Ontg (2 Par): Jo(2I1,* ©} Par) (36) 


ist und diese beiden Residuen einander gleich sind. Es ware also bei der Anwendung des Resi- 
duensatzes auf die Auflésung von Gl. (30) der doppelte Wert dieses Gliedes als Beitrag der beiden 
Pole in + f, einzusetzen. Wegen 27,,-Cof y,, =4o, ist aber dieser Beitrag gleich Null. 

Die Berechnung des Residuums an der Stelle ¢ = 0, dem Ursprung des Einheitskreises, kann 
in diesem zweiten Falle auf dieselbe Weise wie oben erfolgen. Dabei zeigt sich, daB die neuen 
Ausdriicke aus den entsprechenden friitheren ausnahmslos durch den Ersatz von ,, durch 7: y,, 
hervorgehen. Somit wird dann jetzt: 


8 (ion Jue) = EO tg 2 nr EV) Sit 29 Yur) —— Foi) Tg 2nd - (37) 


Die hierin auftretende Reihe ist trotz der Funktion Git sogar absolut konvergent, wenn auch 
nur recht schwach. Nun ist aber nach Gl. (31a) mit x =7,, unda =7y,,: 


Tol2 iy, « ©OF Yue) = J2 ja) + 2° E (A) Col 29 Yn 2 Tia) - (38) 


v=1 
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Die Gl. (37) kann daher im Hinblick auf die Definitionsgleichung (35) fiir die Funktion S(7,,, 7,,) 
im jetzigen Falle auf die folgende und nunmehr endgiiltige Form 


Sosa) = 


= O92 var | 


2 
Jir 


[es Var =Sonl2J1ip > 1; Var = n( 


+ | Toldon =) Tiling) 
0 


i jo no Ji(Jon) 


ioe 


Jon 
25, 


[7a¢j,.). (4 = Sg2y,) 429d en THis} 


v1 


“VE ) 


gebracht werden. Hierin ist die einzige noch auftretende unendliche Reihe im Gegensatz zu der 
in Gl. (37) sehr gut konvergent, denn fiir reelles y,, und fiir reelles Argument x ist 


T(0)| ~ (20)? -e 


—v-In(2 v] (e+ x)) 


(39a) 


Die Reihe (39) ist also fiir die numerische Berechnung in der Tat sehr gut geeignet, da nicht nur 
die e-Funktion, sondern auch die Besselsche Funktion in den entfernteren Gliedern der Reihe 


gemeinsam verkleinernd wirken. 
Die nachstehende Zahlentafel 1 gilt fiir die Werte 7 = 1 (1) 6 und m = 1 (1) 10. 


Zahlentafel 1. Die Werte von gy, und way nach Gl. (35) und (39). 


3 4 5 6 
1 1,251 598 1,398 555 1,451 996 1,480 427 1,497 701 1,509 460 
2 0,766 540 1,166 407 1,296 056 1,362 134 1,402 429 1,429 624 
3 0,503 060 0,906 189 1,130 504 1,240 050 1,304 976 1,348 388 
4 0,996 245 0,572 810 0,952 053 1,112 389 1,204 718 1,265 516 
5 1,286 603 0,356 230 0,746 875 0,976 032 1,100 380 1,180 370 
6 1,502 665 0,741 724 0,477 502 0,825 547 0,990 128 1,092 099 
Fi 1,676 883 0,972 860 0,290 510 0,650 118 0,871 181 0,999 550 
8 1,823 592 1,148 720 0,617 609 0,418 059 0,738 870 0,901 118 
9 1,950 623 15203 227 0,815 359 0,251 328 0,583 351 0,794 282 
10 2,062 795 1,416 863 0,967 426 0,540 443 0,376 390 0,674 748 


In jeder senkrechten Spalte der Zahlentafel1 stehen unterhalb der fettgedruckten Ziffern die 
Werte von p17. 

Nach diesen Rechnungen kann nun die GréBe S(j,,,, 7,,) als bekannt angesehen werden, und 
es ist damit jetzt endgiiltig fiir die Laplace-Transformierte der Temperaturverteilung in der 
warmeleitenden und zugleich dielektrischen Scheibe, wenn die Abkihlung nach auBen gemaB 
den Randbedingungen I., und die Aufheizung im Inneren gemaB der durch II. a, b) beschriebenen 
Verteilung der Warmequellen erfolgt, im ganzen der folgende Ausdruck entstanden: 

Hoe | (*). Hoy — 7 Jolion: ola) 
Ag is <= Jon JiFon) 


} Sonar day) 


aes Pitesti 
O,,:(0, a; Ds a, Se d) p & A 


h d+ So} (2 i») 


(40) 
aad d : d 
hd- oj (= 2s] a gif Xn Sin & 1] 


a a an _ 
\)=8.+2-5: Pied pet hei fy =f). 


Y Jon y 


(a=8.-8=A.-[ +p-( 


Man kann sich leicht davon tiberzeugen, daB diese Losungsgleichung fiir die Temperaturverteilung 
insbesondere auch die beiden Randbedingungen erfiillt, das sind die beiden Forderungen: 


G,,((4,2; 2; 4, +d) = 0 Air d= ee 


(2 Orolo, 23 23 a, + i) oa Plo, d; p34, +4) ae 
g z=+d 
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Dem Grenziibergang / -> co entspricht in Gl. (40) wie vordem die Annahme, da8 auch in den 
beiden Grenzebenen z = + d die Scheibentemperatur fiir alle Zeiten gleich der AuBentemperatur 
ist, d. h., daB in diesen beiden Ebenen @,,, = 0 ist. In der Tat entartet in diesem Grenzfall die 
zweite Zeile in die einfache Formel: . 

rd — aes 


9 * a rs 
ape Coj (dja 7») (x = iba + 2-5. —¢2:s +4) (40a) 


und es wird dann tatsachlich fiir z = + d die transformierte Temperatur 9,,,(0, + 4; f;a, +d) =o. 
Der Grenziibergang d— oo beseitigt in Gl. (40) das komplizierte Bruchglied innerhalb der 
geschweiften Klammer. Er macht aus der Scheibe einen unendlich langen Kreiszylinder mit 


dem Radius a. 


2. Die stationare und die zeitlich veranderliche Temperaturverteilung in der Scheibe 


Nach dem Umkehrsatz berechnet sich nun aus Gl. 1 (40) der eigentliche zeitliche Verlauf 
der Temperatur O,,, selbst in seiner Abhangigkeit von der Zeit ¢ nach der Formel 
y’ +40 


- ef. Og,2, p34, + d)-dp — (y’ > 0). (1) 


OrolQ, 2515 4, a) =~ — 
274 

y’ —14 00 
Diese Temperaturverteilung setzt sich aus einem zeitlich veranderlichen Anteil zusammen, der 
fiir > co verschwindet, und aus einem stationdaren Anteil, der von der Zeit ¢ tiberhaupt nicht 
abhangt. 

Die rein Ortliche Verteilung dieser stationdren Temperatur entspricht bei der Auflésung 
des Integrales (1) nach dem Residuensatz bekanntlich dem Beitrag der Polstelle £ = 0 zu der 
Summe der Residuen. Dann lehrt aber ein Blick auf die Gl. 1(40) wegen t, = 1 fiir ¢=o0, 
da8 im allgemeinen Fall fiir die stationare Temperaturverteilung das Gesetz besteht: 


D5 AOE ("22) . Ho, ; < Joon * @/4) ; 
él ho Cae fm Jon* J1Von) 
te ee ao) oe) eee o=eSa 
| h d- Goi (“ion Tle Es Cin (<ios)| cael Mee aS +d 
a a a 
(2) 


Die Beziehung fiir die zeitlich veranderliche Temperaturverteilung bestimmt sich generell 
stets aus denjenigen Polen des Integranden von Gl. 1(40), die in der linken £-Halbebene auBerhalb 
der imaginaren Achse liegen. Allenfalls koénnen sie noch in diese Achse hineinfallen. Zunachst 
scheint es so, daB zu diesen weiteren Polen auch die Stelle 


a= fie {1 +P 3) 
gehort, d. h. in der #-Ebene die Stelle 
a aa al eae | Poa (3a) 


Bei genauerem Hinsehen erkennt man aber sofort, daf fiir einen verschwindenden Wert von y, 
auch der zweite Summand in der geschweiften Klammer des Integranden von Gl. 1(40) den 
Wert 1 bekommt. Den Nullstellen von y, entsprechen also keine Pole in der -Ebene. 

Dann aber kénnen wahre Pole des Integranden dieser Gleichung nur noch aus den Nullstellen 
der Nenner in der geschweiften Klammer der Gleichung hervorgehen. Sie miiBten sich also aus 
der Relation 


OF 0. ,2.b->00% a, AE d) = 


S(J n? ji >) 


d ees. ff @ po eee Hh 
Gn 2G (= “s) mw ha (4a) 


berechnen lassen, und es ist die Frage, ob es gelingt, solche Werte von y,, zu finden, die die Gl. (4a) 
befriedigen. Nun kénnen sie sicherlich nicht reellwertig sein. Wohl aber wissen wir bereits von 


2 
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Gl. 1(20b) her, daB die Gl. (4a) in der y-Ebene gewi8 unendlich viele rein imaginare Nullstellen 
hat. Setzen wir demgemaB 


— 474-7 (4b) 


x4 a Xn (4c) 
so miissen die dimensionslosen Zahlenwerte x, die Gleichung 
%:tg%,=thd>o (| =o eee) (4d) 


erfiillen, und fiir die an Hand dieser Gleichung berechneten unendlich vielen Nullstellen x, 
besteht dann mit der GroBe y, von Gl. (4a) der Zusammenhang: 


= GPP, (5a) 
An = Et xy ald. (5b) 
Auf die p-Ebene bezogen, legen oo die gesuchten Pole des Integranden von Gl. (1) an den 
Stellen: 
a2 ‘ a2 ile penne 
bin = x2 ee A ene! ; “4 6 
fa, A a2 Jon a2 ee: go) ( ) 


Die Werte #,,, hangen also, wie es zu erwarten steht, sowohl vom Radius a der dielektrischen - 
Scheibe als auch von ihrer Dicke d ab und auBerdem noch iiber den Koeffizienten y? = A/c von 
der spezifischen Warme c, der Warmeleitfahigkeit 1 und der Warmeiibergangszahl H = h- A. 

Um den endgiiltigen Ausdruck fiir den zeitlich veranderlichen Temperaturanteil zu gewinnen, 
schreiben wir zundchst diesen besonderen Teil des een ee explizite an. Es ist danach 


sp ee aa ae Jo(jon* Q/4) _ . Fi 
O(0; B, t; a, = a) car ale 4 a & A [7 “)-H aan Stace lt Ti jon) : wi dor 93) 
y’ +400 


1 i ebt hd Cof (2/a xq) G) 


See: Papa a d ree tea 
een | ha Col (Stn) + tm Gin (rs) 


mit dem Zeichen y, in der in Gl. (3) angegebenen Bedeutung. Es ist nun aber fiir ein p = fy, der 
zugehorige Wert von 


; a 
Kin es (8a) 
Gl, Iny 
1 ; 
ty = bin, (A 12,30). (Sb) 
Andererseits ist 
“laa Cof( 2 y,)+ 27, Sin zee: Sees a ede 
dp a %) y Xn A Xn \e=2, 2 j2- he 2 Ve, [x4 rt (1 ae )] : (9) 


Damit erhalt man dann im ganzen fiir die von Zeit und Ort abhangende Ausgleichstemperatur, 
die zwischen dem tiberall verschwindenden Anfangszustand und der durch Gl. (2) gegebenen 
Endtemperatur vermittelt, die folgende endgiiltige Formel: 


5, / oo J (i » 2) 
O, 0, zntia;t d) = 8. at hd M0”, (Mor), pe ye pecans: DS (Giaeaie 
rans\v , SA bs Orl/1r >a Jon Jon I1y 


€o at Jon) 
5 onion: vlatt Ss cos (z/d + wa)l(xg - sin ¥) pata ae 7: 
«4 +hd(1 + ha) A d 
TS (44 +hd x +hd a(x + S3,] 
= Or.4(0, ee te Crates d) = Oxtar(0, By; {— OO; 4, af a) . (10) 


[h] = 1/cm; [y*] = cm?/s; [A] = W/(cm- grd). 
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Fir ¢—> oo verschwindet @,,an; innerhalb und am Rande der dielektrischen Scheibe fiir jedes 0 
und z aus den Bereichen 0 = @ S a und —d =z < +d. DaB sich fiir t = 0 die Terme Opanc 
und @.,,, tatsachlich zu o erganzen, folgt aus Gl. (1), wennman darin ¢ = 0 setzt und die Gl. (1) 
nach dem Residuensatz auswertet. Das ist statthaft, weil fiir p— oo der Integrand wie 1/f? 
absolut konvergiert. 

Der Ausdruck (10) ist nun zwar fiir alle / aus dem Bereich 0 = ¢ < ov giiltig. Fiir die numeri- 
sche Berechnung ist er jedoch nur fiir einigermafen groBe Werte von ¢t geeignet. Es empfiehlt 
sich daher, fiir kleine Werte von ¢ nach anderen Entwicklungen zu suchen. Dieser Ausweg ge- 
staltet sich besonders einfach, wenn er im Falle der Randbedingung O0,,.,, (a, + d; t) = 0 be- 
gangen wird. Wie wir bereits wissen, geht die dafiir zu wahlende Lésung aus Gl. (7) durch den 
Grenziibergang #— oo hervor. Gerade fiir kleine Erwarmungszeiten ¢ bleibt unter praktischen 
Bedingungen dieser Wechsel in den Randbedingungen auf das numerische Ergebnis ohne Ein- 
fluB, da sich kurze Zeit nach dem Einschalten ein merklicher Warmestrom noch nicht ausgebildet 
haben diirfte. 

Danach bekommen wir es jetzt gemaB Gl. (7) mit dem Ausdruck 


> Jofion 
# Hoy [App\ Gel ee 
Oz 40; z;t;a,+d)=—4 a. Mo% Oe ( 0 ¥ ee on Jomot 
n a 


SyAl Fie Ao J1(Jon) 
y’ +10 Coj ee 
1 ePst “a 
x 1 -\..dp (11) 
20 Le ‘Pp v r d 
y—iwo Co] |— %, 
a 
zu tun. Nun ist 
To) ce ey ad—zZz a 24 +1)4—z 2A+1)d4+2 
oi CARAS ieee! jj ee, eee 
a = 1 -é =1— > (—)* le a. g ; 
oT ie 2 g x Xd fe) : 
Sp} | —4 MH = 
Co} 5 Xn 1 ae e a (12) 
Andererseits erhalt man nach dem Residuensatz sofort das Resultat, das 
y’ +40 " y’ +100 
pt Esp eS pt.d 
1 ; | e = dp Mes L Aalfat e lon’ a t <2 1 . | é Z (13a) 
LES | PiGon i Pa /y-) Jon LET Sy PP dn 
y’—10 y’—10 ) 


ist, wahrend sich im anderen Falle mit # = (s —1)/A und?’ = ¢/A nach bekannten Formeln [2] 


sofort ergibt: 
yticw to) a e.e) 


—aV1l+pA = 
a ewes ae DGD == C.- Es | ee (4 a) nF ds 


2704 p-(1-+ pA) 251 


ae Ea! .J_—_ 2. euA . erfc nee (=e aes S/T eee one a ee 
=— Ze erfe (. Fale Ce Lie (. gm yA) + é erlc (- am + yt/A } (13b) 
(6>1+yA/’). 


Im vorliegenden Fall ist dabei 
A = a*l(y? Fe) ’ (14a) 
a=((2A+1)d+2)/a-j,,, (14b) 


erfc (x) = 2f/n- fe-* - d& = 1 — erf (x) (14¢) 


und hierin gilt fiir die Funktion erfc(x) die folgende asymptotische Entwicklung: 


Vz - erfc(x) wher : — = L = Nota ( a glee) es a =) +0[55 ) (14d) 
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Damit kann der Ausdruck (11) fiir die Ausgleichstemperatur auf die folgende, gerade bei kleinen 
Zeiten ¢ sehr gut brauchbare Form gebracht werden: 


lo i OV 2 Wey. 
Oi Oreo areae d)=—4a2.% (@) . OES 


eA Ao 


1 
i a (:—+)a—z 
ea SO tad ae a ee 


Jon Ju Jon) =o Vty? 
1 
(2A+1)d—z . A —|¢d—- I 
1 ~ pe yl ene mae eh LC -f ( ig a IL 4 ies 
Sets ete | ee Saag pe 
2 fro yir? 
(24 +1)d+2 (i+s)a+: 
10n yt y2 
Ayn) oy erfc ee aelne 2 
- (244+1)d—z (: gf s)a—s yiy? 
roel (=) "le & erfc aoe On 2 
en yi a 
1 
(2A +1)d+2., (i+ sja+s a 
iy 2 : ty? 
aie f ; ed 3 dons He (15) 
yp 


In der Gl. (15) liegt damit auch eine Formel fiir die Berechnung der Ausgleichstemperatur vor, 
die den erheblichen Nachteil der Entwicklung (10) nicht mehr aufweist. Sie lieBe sich im Hin- 
blick auf die Gl. (14d) fiir die allerersten Zeitteilchen sogar noch weiter vereinfachen. 


3. Angaben iiber die Stoffkonstanten der dielektrischen Scheibe 


Der numerischen Berechnung des Erwarmungsvorganges in der dielektrischen Scheibe, etwa 
nach der zuletzt angegebenen Formel, steht die Schwierigkeit entgegen, daB fiir die maBgebenden 
Materialkonstanten c, A und y? = A/c nur unsichere Zahlenangaben gemacht werden kénnen. 
Zwar sind in den letzten beiden Jahrzehnten gerade mit pordsen Materialien zahlreiche Versuche 
[1] und [5] angestellt worden, um dariiber zutreffende Werte zu erfahren, aber von Fall zu 
Fall diirften sie doch weit auseinander gehen. 

In der nachstehenden Tabelle sind Durchschnittswerte der oben genannten drei Konstanten 
den gut bekannten Werten dieser Konstanten fiir Kupfer gegeniibergestellt worden. Es iiber- 
rascht nicht, daB fiir das porése Material die Warmeleitfahigkeit fast um zwei Zehnerpotenzen 
geringer ist als fiir das gut leitende Kupfer. Das zeitliche Ansteigen der Temperatur bei der Er- 
warmung hangt aber, wie z. B. gerade die letzte Gl. 2(15) erkennen laBt, keineswegs alleine von 
A ab, sondern vornehmlich von der Konstanten A/c, und hier macht der Unterschied nur noch 
eine Zehnerpotenz aus, weil natiirlich die Temperaturerhéhung je Grad bei dem Isolierstoff 
ebenfalls eine weit geringere Warmemenge erfordert als bei einem Metall. 


Zahlentafel 2. Die Stoffkonstantenc, 4, Hund y?sowiex und tg 6. 


Stoffkonstanten 


A 
W/(cm,- grd) 


c 
— Ws/(cm' - grd) 


Material 


Poréses Material 
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Eine andere Materialkonstante des Einsatzes, die fiir die vorliegende Rechnung benotigt 
wird, ist die elektrische Leitfahigkeit des Scheibenmaterials. Dariiber sind jedoch geniigend 
zuverlassige Zahlenangaben zu erhalten, vor allem wenn dabei nicht so sehr nach dem Wert 
von x selbst gefragt wird, sondern nach dem Verlustwinkel 6. Da definitionsgemaB 


tg 0 = x/( , &) (1) 


ist und hier wie schon eingangs ¢, ~ 1 gesetzt werden kann, so laBt sich in allen vorstehenden 
Gleichungen sehr rasch von der Konstanten x zur Konstanten 6 tibergehen. Man hat namlich 
nur notig, in den einzelnen Formeln 


°" =tgd- oe 27 te O08 Es in Ohm- grd/W (2) 
zu setzen mit A) als der Vakuum Wellenlange. 

Fir die zahlenmaBige Berechnung der Nullstellen x, oder &) von Gl. 2(4d) oder 1(20£) 
mége noch darauf hingewiesen werden, da sie bereits in dem unter 3) erwahnten Buch fiir eine 
groBe Zahl von Werten des Produktes / d, dort mit C bezeichnet, durchgefiihrt worden ist. Am 
SchluB des Buches wird auch eine Tabelle dariiber gebracht, die jeweils die ersten sechs Null- 
stellen enthalt. Auf Grund der Zahlentabelle II kann im vorliegenden Falle fiir h d etwa der 


Zahlenwert 0,020 angenommen werden. Es sind dann 


H=—=0;14105 Np) 31470 > Lg 6,2804,, 
Ap 042005 4, ==.12,5080., Ng 15) 70021. 
Die auBerdem erforderlichen Zahlenwerte von 7,, und j,, fiir (7, 7) = 1, 2,3... sind aus der 


Hohlleitertheorie hinlanglich bekannt. 
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Zur Bestimmung des Laufschienenfeldes am wandernden Lichtbogen 
Von 
Dipl.-Ing. DirK Hesse, Darmstadt 


(Mitteilung aus dem Institut fiir Hochspannungs- und MeBtechnik der T. H. Darmstadt) 


Mit 30 Textabbildungen 


(Eingegangen am 29. February 1960) 


Ubersicht: Messungen der Lichtbogenwanderungsgeschwindigkeit als Funktion des Stromes sind nur 
dann miteinander vergleichbar, wenn auBer der Induktion der Blasspule die GréBe des Laufschienenfeldes 
beriicksichtigt wird. Hier soll daher ein Naherungsverfahren beschrieben werden, das eine Bestimmung 
des Laufschienenfeldes in der Umgebung des Lichtbogens erméglicht. Die danach ermittelte Induktionsver- 
teilung wird mit der durch einen Hallgenerator ausgemessenen verglichen. Dabei wird der Einflu® der Strom- 
linienverteilung in der Laufschiene bei punktformigem Stromaustritt am LichtbogenfuBpunkt erértert. Ein 
Vergleich der an einem Modellichtbogen aufgenommenen Induktionsverteilung mit der am wirklichen Licht- 
bogen gemessenen lat Schliisse auf die Breite des stromfiihrenden Kanals und die Stromdichteverteilung im 


Bogen zu. 
1. Aufteilung im Lichtbogeneigenfeld und Laufschienenfeld 


Fiir die folgenden Betrachtungen wird, wenn nicht anders angegeben, das im Bild 1 fest- 
gelegte Koordinatensystem benutzt. Der Koordinatenursprung liegt in der Mitte der linken 
Schiene. Die Stromeinspeisung in die Schienen erfolge bei y = oo. 

Das magnetische Feld am Lichtbogen kann man sich aus zwei Anteilen zusammengesetzt 
vorstellen: Der erste besteht aus dem durch die Bewegung der Ladungstrager im Lichtbogen 

at 7 verursachten Feld, dem sogenannten Lichtbogen 

a eigenfeld. Nimmt man die Gestalt des Lichtbogens- 

als Zylinder an, so sind seine magnetischen Feldlinien 

Kreise. Auf den Lichtbogen wirken daher von allen 
Seiten radial gerichtete Krafte ein, die abgesehen 
von einer komprimierenden Wirkung auf das Licht- 
bogenplasma keine Veranderung der raumlichen 
Lage des Lichtbogens zur Folge haben. Der zweite 
Anteil stammt vom Strom in den Laufschienen. 
Die Kraftlinien des Laufschienenfeldes durchsetzen 
den ganzen Lichtbogen in angenahert gleicher 
Richtung. Sie bestimmen zusammen mit dem Licht- 


Sa bogenstrom J und der Lichtbogenlange e die GréBe 


Bild1. Lage des Koordinatensystems zu den Laufschienen. = : 
b Schienenbreite; d Schienenstarke; e Schienenabstand = Licht- der Blaskraft Pim: 
bogenlange; a@ Abstand von Schienenmitte zu Schienenmitte. 
P= [J c Binet] ‘e (1) 


B miter St der mittlere Induktionswert, der sich aus der Integration des Laufschienenfeldes iiber 
der Lichtbogenlange e ergibt. 


2. Das Lichtbogeneigenfeld 


Bei Annahme zylinderférmiger Gestalt des Lichtbogens gilt fiir die Induktionsverteilung 
auBerhalb des Bogens (Bild 2) 


le 1 
b= 2 (iA (2) 


Die Induktionsverteilung innerhalb des Bogens ergibt sich zu 


> Mo L MoI 
B= ee (3) 


207" 2S 


Beide Gleichungen gelten fiir einen oo langen geradlinigen Lichtbogen. 
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Fir einen Bogen endlicher Lange (Bild 3) folgt aus dem Brot-Savartschen Gesetz die In- 
duktionsverteilung auBerhalb des Bogens zu 
PETE MT |e é 


x . —_— *| 
SS Sill Aik — ; : 
Ge | arc tan 5 + sin arc tan 5 i (4) 


wenn die x-Koordinate fiir Punkt P zwischen O und ¢ liegt. 


8 


LILI 


+2 


Bild 2. Induktionsverlauf des Lichtbogeneigenfeldes. Bild 3. Zur Bestimmung der Induktion auBerhalb eines endlich 
langen Lichtbogens. 


Lichtbogenltinge e=ee 


Lichtbogenldnge 
om 


2 
0 0 +1,0 2,0 cm 
; beg —e 
Bild 4. Induktionsverlauf im Abstand y=1cm von der Licht- Bild 5. Induktionsverteilung bei endlicher Lange und Zylinderform 
bogenachse. B = f(x) fir J = 100A u. Lichtbogenlange 1 cm. des Bogens. B= f(x; y) fir x=0 bisx=eu.y=rbisy= +o. 


Im Bild 4 ist der Verlauf der Induktion langs des Bogens, also B = f(x), fiir einen Abstand 
4 = 1cm von der Lichtbogenachse dargestellt. Wahrend die Induktion B bei einem unendlich 
langen geradlinigen Bogen konstant ist, fallt sie bei endlicher Lange von der Bogenmitte nach den 
Bogenenden hin ab. In Bild 5 ist die raumliche Induktionsverteilung B = f(x; y) des Lichtbogen- 
eigenfeldes skizziert. 


3. Das Laufschienenfeld ohne Beriicksichtigung des Laufschienenprofils 


Alle Stromlinien seien zu einem Stromfaden in Leitermitte zusammengefaBt. In der Hohe 
des Lichtbogens soll der Stromfaden rechtwinklig abbiegen. Der Lichtbogen wird dabei als 
unendlich diinner Stromfaden angenommen. Die Induktion des Blasfeldes in Hohe der Licht- 
bogenachse (y = 0) ergibt sich aus dem Brot-Savartschen Gesetz. Fiir die von der linken Lauf- 
schiene im Punkt P, (Bild 6) hervorgerufene Induktion gilt: 


co 


pele ai x dy (5) 


4% 2% 1 y2)3/2° 


fe) 


Als Gesamtinduktion in P,, hervorgerufen durch den Strom in beiden Laufschienen, erhalt man 


B= 1 oa ) (6a) 


4a \w% a= *% 
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Fiir einen Punkt P, mit der Ordinate y = + oo ist die Induktion gerade doppelt so groB, némlich 


s Ho 1 1 
B= M12 +5) (6b) 
Fiir die Kraft auf den ganzen Lichtbogen folgt durch Integration iiber die Lichtbogenlange 
a—b/2 
_ top f (4 & _ He yp 28? 
tienes | G+<5)}- tome! In iD yp (7) 


[2 


Auf einen Lichtbogen von 1000 A mit einer Lange e = 1 cm wirkt danach bei einer Schienenbreite 
von b = 0,5 cm eine Blaskraft von 33 p. 
Bei Beriicksichtigung der endlichen Ausdehnung des Licht- 
*y bogens ist auch die Laufschieneninduktion als Funktion des Ab- 
| standes von der Lichtbogenachse von Interesse. B = f(y) bei kon- 
stantem x-Wert erhalt man durch Aufteilung in zwei Bereiche 
(Bild 7). 
Bereich: y= bissy — cos 
Bereich 27 y = 0sbisyy— 


Die Induktion, die von den Laufschienenabschnitten im ersten 
Bereich im Punkt P erzeugt wird, bestimmt sich analog nach 


Gl. (6a) zu 
a He (2 see ). 


4% x a—* 


Die Induktion, die von den Laufschienenabschnitten im zweiten 
Bereich hervorgerufen wird, erhalt man mit Hilfe der Gl. (5). Nur 
muB die obere Grenze des Integrals jetzt y = w heiBen. Die Ge- 
samtinduktion im Punkt P ergibt sich aus der Uberlagerung der 
Induktionsanteile beider Bereiche zu 


1) 


Peal ) 
1 = marae -- sin are tan =) + é 


: jf 
ait sin arc tan Z 
Bild 6. Zur Bestimmung des Laufschie- 47 |% —* ( ae t a, 
nenfeldes in Hdhe der Lichtbogenachse ( 8) 


(y = 0). 


Bild 7. Zur Bestimmung des Laufschienenfeldes Bild 8. Induktionsverteilung des Laufschienenfeldes bei punktformigem Stromaustritt 
fiir beliebige y-Werte. aus den Laufschienen. B = f(x; y) fir z= 0. 


In Bild 8 ist die rdumliche Induktionsverteilung des Laufschienenfeldes als B = f(x; y) fiir 
z= 0 skizziert. Den Verlauf der Gesamtinduktion in der Umgebung des Lichtbogens erhalt 
man, wenn man Bild 5 und 8 iiberlagert. 
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4. Das Laufschienenfeld unter Beriicksichtigung des Laufschienenprofils 


Bei kleinen Schienenabstanden ist es nicht mehr zulassig, die Induktion zwischen den Lauf- 
schienen durch Konzentration aller Stromlinien zu einem Stromfaden im Leitermittelpunkt zu 
bestimmen. Wenn ein Rechteckleiter von y = + co bis y = —oo gleichmakig vom Strom 


magn.Kraftlinie 


Bild 9. Zur Bestimmung der Induktion bei Rechteckleitern Bild 10. Lage der Stromfaden zur naherungsweisen Bestimmung 
mit gleichmaBiger Stromverteilung. des Induktionsverlaufes zwischen den Laufschienen. 


durchflossen wird, so ist die im Abstand x und z vom Koordinatenursprung (Bild 9) erzeugte In- 


duktion in z-Richtung 
d 2 
(: 2 b 
7a (« ~) 


( b \2 =) b 2 
Pea d vt) 4249 (« =| 
1 = oe : (: | jin Z : (: )In 


4 d d £ d d 

zg+— zg —— Z.— g—— 

are tg — arc tg ; — (« —) arc tg sana te a tale (ga) 
Ss aoe x -- — & ei NTP eR 


( zn y+ d a 
pei esha sae pais 
lM iE || ae 2 4 b 2 b 2 
B= Se - a ) AEG 2(x | 7) are tg ; + 2(¥——) are tg mie (gb) 
sats Baas aes Ges 
5 4 sie - (fares 


Nimmt man wieder an, daB alle Stromlinien in Hohe des Lichtbogens rechtwinklig abbiegen, so 
betragt die Induktion in der Lichtbogenachse gerade die Halfte der nach Gl. (gb) berechneten 
Induktion. Denn die Laufschienen sind jetzt nur von y = + oo bis y = 0 vom Strom durch- 
flossen. Die Gesamtinduktion des Laufschienenfeldes ergibt sich wieder aus Uberlagerung beider 
Laufschienenanteile. Das gleiche Verfahren wurde auch von KUHNERT [1] angewendet, ohne es 
jedoch durch Messung nachzupriifen. 

Wahrend die Bestimmung der Blasfeldinduktion fiir die Ordinaten y = + co und y= 0 
einfach ist, erfordert sie fiir andere y-Werte einen gréBeren Aufwand. Der Induktionsverlauf 
B = f(y) laBt sich aber annahern, wenn man alle Stromlinien jedes Rechteckleiters zu einem 
Stromfaden an einem solchen Ort zusammenfaBt, daB zwischen den Laufschienen eine der wahren 
Induktionsverteilung méglichst ahnliche entsteht. Um den Abstand von diesen beiden Strom- 
faden, in folgendem mit a, bezeichnet, zu ermitteln, geht man von Gl. (6b) aus. Danach gilt fiir 
den Induktionsverlauf zwischen den Laufschienen 


pati (+ : ) (10) 
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wenn der Strom von y = + 06 bis y = —oo flieBt (Bild 10). Berechnet man nach Gl. (gb) die 


wahre GroBe der Induktion in der Mitte zwischen den Laufschienen, so 1aBt sich mit ihr aus 
Gl. (10) a, bestimmen, wenn x mit a,/2 eingesetzt wird. Der Ursprung des Koordinatensystems 
ist dabei mit dem Stromfaden der linken Laufschiene zusammengelegt. 


Bs df (11) 


a, = : 
: a Binitte 


In Bild 11 sind fiir ein bestimmtes Schienenprofil die nach verschiedenen Verfahren ermittelten 
Induktionsverteilungen zwischen den Laufschienen dargestellt. Kurve1 folgt aus Gl. (gb). 
Kurve 2 ist das Ergebnis der oben angegebenen Naherungslésung. In der Mitte zwischen den 


1§ 


1 wahre Induktionsverteilung 
nach Gl. (gb); 


2 Naherung nach Gl. (10) 
4, = 37 mam; 


10 


3 Naherung nach E1IpINGER 
und RIEDER, 4, = 46,5 mm; 


4 Stromfaden in Leitermitte 
konzentriert a = 57 mm. 


pO are are 
Tm an lion e TS — 


e=7mm -¥—— @, ——»| 
Bild 11. Induktionsverteilung zwischen zwei Rechteckleitern Bild 12. Versuchsanordnung zur Messung des Laufschienen- 
bei J=100A. B=f(x) fir z=o0. Schienenbreite b=50mm; feldes in der Umgebung des Lichtbogens. 


Schienenstarke d= 5 mm; Schienenabstand e = 7 mm, 


Schienen stimmt nach Voraussetzung der wahre Induktionswert mit dem der Naherungslésung 
iiberein. Kurve 3 gibt den Induktionsverlauf nach einem von EIDINGER und RIEDER [2] an- 
gegebenen Naherungsverfahren wieder. Kurve 4 entsteht dann, wenn man die Abmessungen 
des Laufschienenprofils vernachlassigt und alle Stromfaden in Leitermitte konzentriert. 

Hat man den Ort gefunden, in dem alle Stromlinien zu einem Stromfaden vereinigt sind, so 
1aBt sich der Induktionsverlauf 6 = f(y) fiir Rechteckleiter nach Gl. (8) ermitteln, wenn man 
statt a jetzt a, einsetzt 


[ee | : [ | sin are tan “ ) + é ( + sin are tan Z I. (12) 
4 | % x a—* ay — % 

Fiir zwei verschieden breite Reckteckprofile wurde der nach Gl. (12) berechnete Induktions- 
verlauf des Laufschienenfeldes B = f(y) fiir « = a,/2 und z=o0 durch Messungen mit dem 
Hallgenerator nachgepriift. Der Schienenabstand e betrug fiir beide Profile 7mm, der Strom 
J = 100 A. Um bei der Messung des Laufschienenfeldes das Lichtbogeneigenfeld nicht mitzu- 
messen, wurde an Stelle des Lichtbogens ein Kupferbiigel so eingesetzt (Bild 12), daB die von 
dem ihn durchflieBenden Strom erzeugten Feldlinien keine senkrechte Komponente zur Ebene 
des Hallgenerators bilden. Somit kann von ihnen auch keine Hallspannung erzeugt werden. 

In Bild 13 stellt Kurve 1 den nach Gl. (12) berechneten Induktionsverlauf fiir das breite 
Profil (6 = 50mm) dar. Bei Ansatz des Kupferbiigels‘an den Laufschienenkanten wurde der 
Induktionsverlauf gemessen, wie ihn Kurve 2 wiedergibt. Das Ansteigen der Induktion in der 
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Nahe von y = 0 laBt sich folgendermaBen erklaren: Das vom Lichtbogen mit der Lange e hervor- 
gerufene Eigenfeld ist durch die Anordnung des Kupferbiigels zwar ausgeschaltet, aber die 
Stromlinien in jeder Laufschiene, die eine axiale Verlangerung des Lichtbogens bilden, rufen eine 
ahnliche Feldverteilung wie der Lichtbogen hervor. Dieses Feld iibt aber keine Kraft auf den 
wandernden Lichtbogen aus. Die wahre GréBe des Blasfeldes ]aBt sich in dieser Versuchsan- 
ordnung daher nur fiir y = 0 bestimmen. 


/ 
| | 
S ld 
25 f| KS ce ts + Aquipotentiallinien 
| 


ted | 
20 k 


10 TF 


-10 boa tt tt 
43 ade = BIS OE A A a St a EEL) 
cm 
Bild 13. Induktionsverteilung des Laufschienenfeldes B = f(y) fiir x= a,/2 
bei punktformigem Stromaustritt an der Stelle y=o. J=100A; Lauf- 
schienenprofil: b= 50mm d= 5 mm; Schienenabstand: e = 7 mm. 


Bild 14 (rechts). Stromlinienvertcilung in einem rechteckigen Leiter bei 
punktformigem Austritt des Stromes an der Kante des Leiters. Aufgenom- 
men im elektrolytischen Trog. 


Um den Induktionsanteil dieser Stromlinien auch bei anderen y-Werten auszuschalten, 
wurde der Drahtbiigel mit seinen Enden im Koordinatenursprung fiir die linke Schiene und im 
Abstand a, in der rechten Schiene angeschlossen (gestrichelter Verlauf des Biigels in Bild 12). Die 
bei dieser Anordnung erhaltenen MeBwerte sind in Kurve 3 (Bild 13) wiedergegeben. Es ergibt 
sich eine verhaltnismaBig gute Ubereinstimmung mit dem berechneten Verlauf. 

Bei der zweiten MeBanordnung liegt aber eine grundsatzlich andere Stromverteilung in den 
Laufschienen in der Nahe von y = 0 vor als es bei der ersten MeBanordnung der Fall ist. Man 
erkennt dies aus den unterschiedlichen Induktionswerten fiir y = 0 (B = 8,5 G bei der ersten 
MeBanordnung, B = 11,5 G nach der zweiten MeBanordnung). Die tatsachlich in der Lichtbogen- 
achse vorhandene Blasfeldinduktion ist kleiner als die halbe Induktion, die sich fiir y = oo ergibt. 

Die Stromlinienverteilung bei punktformigem Stromaustritt aus der Laufschienenkante, 
aufgenommen im elektrolytischen Trog, ist in Bild 14 wiedergegeben. Die an der rechten AuBen- 
kante des Schienenprofils liegenden Stromlinien eilen iiber den Lichtbogenansatzpunkt hinaus 
und kehren nach Durchlaufen einer Schleife zum Stromaustrittspunkt zuriick. Die von den riick- 
laufenden Stromlinien am Lichtbogen erzeugte Induktion schwacht die von den vorlaufenden 
Stromlinien hervorgerufene. Die an der AuBenkante der Laufschiene tiber den Lichtbogenfu- 
punkt hinaus vorlaufenden Stromlinien verringern zwar diese Induktionsabsenkung der riick- 
laufenden Stromlinien. Sie kénnen sie jedoch wegen ihres gr6Beren Abstandes vom Lichtbogen 
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nicht ganz aufheben. Der induktionsmindernde Einflu8 der riicklaufigen Stromlinien wird daher 
um so gréBer werden, je breiter die Laufschienen sind. 

Fiir den Anschlu8 des Kupferbiigels im Koordinatenursprung 7 = 0 und x = a, ist die 
Stromlinienverteilung in Bild 15 skizziert. Das Zuriicklaufen der Stromlinien geschieht in der 
Nahe der y-Achse. Die von ihnen erzeugte Induktionsabsenkung am Lichtbogen wird jetzt durch 
die an der Laufschienenkante vorlaufenden Stromlinien weitgehend kompensiert. Es muB sich 
in diesem Fall fiir y = 0 eine héhere Induktion ergeben. 


Ari 6=10mm 
a= 3mm 


«S. vomHaligenerator 
| angezeigter Mittelwert 


Profil 6=50mm 

4 d= 5mm 

| SvomHallgenerator 
ngezeigter Mittelwert ~Bnitte 


20 


Schienenabstand e 


Bild 15. Skizzierter Stromlinienverlauf bei AnschluB Bild 16. Induktionsverteilung zwischen den Laufschienen B = f(x) 
des Kupferbiigels in + =o u.x =a. fiir y= oo und J = 100A. 


Bei der Messung am breiten Profil (b = 50mm) konnte der inhomogene Feldverlauf langs 
des Hallgenerators vernachlassigt werden (Bild 16). Dagegen liegt der beim schmalen Profil 
(6 = 10 mm) gemessene Wert deutlich iiber dem fiir x = a,/2 errechneten (fiir y = co berechnet: 
Bitte = 53.0 G, gemessen: 56G). In Bild 17 ist Kurve 1 wieder der nach Gl. 12 berechnete 
Induktionsverlauf. Kurve 2 gibt die bei AnschluB des Kupferbiigels an den Laufschienenkanten 
gemessenen Werte wieder. Wegen der geringen Schienenbreite ist jetzt in der Umgebung von 
y = 0 keine Induktionsiiberhéhung mehr festzustellen. Fiir y = 0 wird ein Wert von B = 20G 
gemessen statt 28 G, die sich ergeben wiirden, wenn am Lichtbogen die Induktion gerade halb 
so groB ware wie fiir y = o. 

In Bild 18 sind die Ergebnisse der Laufschienenfeldberechnung fiir die Schienenbreite b = 
10mm mit denen von MULLER [3] verglichen. Es ist B,,i,,. als Funktion des Laufschienen- 
abstandes aufgetragen. Da MULLER das Blasfeld nur fiir eine unendlich diinne Schienenstarke 
angibt, wurde bei den eigenen Berechnungen die Schienenstarke d als Parameter gewahlt. Mit 
abnehmender Schienenstarke nahern sich die Induktionswerte der von MULLER angegebenen 
Kurve. 

Bild 19 zeigt den EinfluB der Schienenstarke d auf die GréBe des Laufschienenfeldes deut- 
licher. Hier ist B = /(d) fiir eine Schienenbreite b = 10 mm und einen Schienenabstand e = 7mm 
aufgetragen. Mit wachsender Laufschienenstarke fallt die Induktion etwa hyperbolisch ab. 

Tragt man die Induktion B,,;,; als Funktion der Schienenbreite } auf (Bild 20), so nehmen 
die Abweichungen zwischen MULLERs und den eigenen Berechnungen mit breiter werdendem 
Profil zu. Bei einer Schienenbreite von 6 = 50mm wurden fiir B,,iy,; 8,5 G gemessen, die 
eigene Berechnung ergab 11 G und nach MULLER folgt, wenn man die von ihm angegebenen Kur- 
ven entsprechend extrapoliert, 16G. Die induktionsvermindernde Wirkung der riicklaufigen 
Stromfaden in Fu8punktnahe, bedingt die mit wachsender Schienenbreite gréBer werdenden 
Abweichungen von der Rechnung. 
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Bild 17. Induktionsverteilung des Laufschienenfeldes bei punkt- 

formigem Stromaustritt an der Stelle y=o. J=100A; 

Laufschienenprofil: 6 = 10mm d= 3mm; Schienenabstand: 
é=7mm. 
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0 10 20 30 _@ 40 Bild 20. Mittlere Induktion des Laufschienenfeldes am Lichtbogen 
™mMm (y = 0, = 0) in Abhangigkeit von der Schienenbreite. By) js¢.)= f(0); 
Bild 19. Induktion des Laufschienenfeldes am Lichtbogen J = 100 A; Schienenabstand e = 7 mm; 
(y=o, =o) in Abhangigkeit von der Schienenstarke d, Berechnung nach MtjttER d=omm _ (1) 
B=f(d); J = 100A; Schienenbreite 6 = 10mm; Schienen- Eigene Berechnung d=o0,6mm (2) 
abstand e = 7 mm. undd=3mm _ (3). 


5. Verlauf der Gesamtinduktion am Modellichtbogen 


Addiert man Laufschienen- und Lichtbogeneigenfeld, so erhalt man die endgiiltige Induk- 
tionsverteilung zwischen den Laufschienen am Lichtbogen. Eine Kontrolle dieser rechnerischen 
Zusammensetzung wurde durch Messung mit dem Hallgenerator an einem beide Laufschienen 
verbindenden Modellichtbogen durchgefiihrt. Als Modellichtbogen diente ein Kupferbolzen von 
5mm Durchmesser, der mit einem Zapfen in eine Bohrung der Laufschiene eingesetzt wurde. 
Der Laufschienenabstand betrug e = 7 mm, die Laufschienenbreite ) = 10 mm und die Lauf- 
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schienenstarke d= 3mm. Der Hallgenerator (Type SBV 511, Abmessungen 24mm), der 
als TangentialfeldmeBsonde ausgebildet ist, konnte bis unmittelbar an den Rundleiter herange- 
fiihrt werden. Er wurde parallel zur y-Achse in der Hohe z = 0 
in der Mitte zwischen den Laufschienen verschoben (Bild 21). 

In Kurve 4 in Bild 22 sind die hierbei erhaltenen MeBwerte 
wiedergegeben. Kurve 3 ist aus der rechnerischen Addition von 
Laufschienenfeld (Kurve 1) und dem Lichtbogeneigenfeld (Kurve 2) 
entstanden. Das Lichtbogeneigenfeld wurde nach Gl. (4) fiir x = a,/2 
berechnet. Fiir groBe Entfernungen vom Lichtbogen, miissen 
die Kurven 3 und 4 mit Kurve 1 zusammenfallen, wenn das Licht- 
bogeneigenfeld auf Null abgeklungen ist. 

Wahrend die Kurven 1—4 fiir z = 0 gelten, wurde im Hinblick 
auf die Messungen mit dem Hallgenerator am wandernden Licht- 
bogen der Verlauf der Gesamtinduktion in z-Richtung fiir x = a,/2 
und z = 2,5 mm aufgenommen (Kurve 5). In Bild 23 sind auBer 
der Gesamtinduktion (Kurve 3) das Lichtbogeneigenfeld (Kurve 2) 
Bild 21. Messung der Gesamtinduktion und das Blasfeld (Kurve 1) (wieder mit dem Hallgenerator aus- 


am Modellichtbogen mit dem Hallgene- i 
rator (HG). gemessen) fiir z = 2,5 mm dargestellt. 
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Bild 22. Verlauf der Induktion in z-Richtung am Lichtbogen in Abh. vom Licht ogenabstand By, = f(y). Profil: d=3mm b)= 10mm; 
e= 7mm 4, = 15,1mm; Lichtbogendurchmesser 5 mm; J = 100 A; z fiir 1 bis 4 =< o. Kurve 1: Schienenblasfeld in der Mitte zwische . 
den Schienen aufgenommen mit dem Hallgenerator (kein Lichtbogen vorhanden) ; --— Kurve 2: Lichtbogeneigenfeld; —+— Kurve 3: Gosamites 
Induktionsverlauf By, = f(y). Addition von 1 u, 2; «rece Kurve 4: Gesamter Induktionsverlauf nachgemessen mit dem Mallconcpaiars 

—++— Kurve 5: Gesamter Induktionsverlauf B, = f(y) fiir z= 2,5 mm gemessen mit dem Hallgenerator. : a. 
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6. Verlauf der Gesamtinduktion am wandernden Lichtbogen 


Wahrend der Hallgenerator bei der Bestimmung des Induktionsverlaufes in der Umgebung 
des Modellichtbogens parallel zur y-Achse bewegt wurde, liegt er nun fest und der Lichtbogen 
lauft tiber ihn weg. Die vom Hallgenerator erzeugte Spannung wird einer Oszillographenschleife 
zugefiihrt, die den zeitlichen Verlauf der Induktion B = f(t) wiedergibt. Bei konstanter Wan- 
derungsgeschwindigkeit des Lichtbogens laBt sich die zeitliche in eine rdumliche Induktionsver- 
teilung umwandeln, so daB diese mit der des Modellichtbogens verglichen werden kann. 
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Bild 23. Verlauf der Induktion am Lichtbogen in Abh. vom Lichtbogenabstand Bz = f(y). Bild 24. Messung der Gesamtinduktion 
Profil: d=3mm; bJ=1iomm; e=7mm; Lichtbogendurchmesser 5mm; J = 100A; am wandernden Lichtbogen. 
furs (2.== 2,5 mm. Kurve 1: Schienenblasfeld gemessen; — — Kurve 2: Lichtbogen- 
eigenfeld berechnet; —-— Kurve 3: Gesamter Induktionsverlauf gemessen. 


Der Hallgenerator wurde bei diesem Versuch in der Mitte zwischen den Laufschienen in der 
Hohe z = 2,5 mm fest angebracht, so daB seine Fassung biindig mit der Laufschienenkante ab- 
schneidet (Bild 24). Ein besonderer Schutz der Hallgenerators gegen die heiBen Lichtbogengase 
erwies sich bei der hohen Wanderungsgeschwindigkeit des Bogens als nicht erforderlich. Es 
wurden die gleichen Schienenabmessungen und derselbe Schienenabstand wie beim Modellversuch 
gewahlt. Aus dem Oszillogramm (Bild 25) wird fiir den Lichtbogenstrom ein Wert von 1230 A 
entnommen. Die Geschwindigkeit des Lichtbogens, die einmal aus dem Anstieg der Weg-Zeit- 
linie und zum anderen mit unterhalb der Laufschienen angebrachten Photozellen bestimmt 
wurde [4], betragt 90 m/s. Wahrend bei Strémen iiber 1000 A keine Anderungen in der Geschwin- 
digkeit beim Uberlaufen des Hallgenerators auftraten, betrug die Geschwindigkeitsabsenkung 
bei Stré6men unter 1000 A bis zu 50% bei J = 300 A. Der Lichtbogen wurde von der Seite und von 
vorn durch eine sich schnelldrehende Schlitzscheibe photographiert (Bild 26). Er zeigt eine stark 
von der Zylinderform abweichende Gestalt, deren Entstehung in [4] naher erlautert ist. Der 
Hallgenerator ist genau iiber der 3. Photozelle angebracht. Die Lage der Photozellen sowie die 
Laufrichtung des Bogens sind durch Leuchtmarken auf der Photographie gekennzeichnet. Wenn 
auch in der Weg-Zeitlinie keine St6rung in der Lichtbogenwanderung zu erkennen ist, so ver- 
gréBert der Bogen doch beim Uberlaufen des Hallgenerators seine Breite und verliert gleichzeitig 
seine sonst so symmetrische Form, wie man der seitlichen Photographie entnehmen kann. Bei 
der Photographie von vorn ist der Lichtbogen durch den Hallgenerator abgedeckt, so daB von 
ihm nur ein Stiick sichtbar ist. Die der Photographie entnommene Breite 0 des Bogens unter 
dem Hallgenerator betragt 4cm, die Dicke d = 1,9 cm, gegeniiber etwa b = 2,5cm und d= 2,2cm 
vor und hinter dem Hallgenerator. 

In Bild 27 ist der dem Oszillogramm entnommene zeitliche Verlauf der Induktion im ver- 
groBerten MaBstab aufgetragen. Gegeniiber dem Modellichtbogen besitzt der wahre Lichtbogen 


+2 


1 LichtHoe Archiv fiir 
D. Hesse: Zur Bestimmung des Laufschienenfeldes am wandernden Lichtbogen _pjektrotechnik 


Lichtbogenstrom 
1=7030A 


Induktionsverlouf 
amLichtbogen 


Bild 27. Induktionsverteilung am wandernden Lichtbogen 
B= f(t) (dem Oszillogramm entnommen). 


von der Seite von vorn 


Bild 26. Photographie des wandernden Lichtbogens. 
b Lichtbogenbreite; d Lichtbogendicke. 


Bild 28. ——— Jnduktionsverteilung am wandernden 
Lichtbogen B= f(y) (umgerechnet aus Bild 28); 
---- Laufschienenfeldverlauf B= f(y) (entnommen 
aus Bild 24, Kurve 1 auf J = 1230 A umgerechnet), 


breite Bild 28. 


wesentlich gréBere Abmessungen. Daher sind die Induktionsspitzen am Bogenrand wesentlich 
kleiner. Rechnet man den zeitlichen Induktionsverlauf in eine raéumliche Induktionsverteilung 
um (Bild 28), so ist der Abstand der beiden Induktionsspitzen gleich der elektrischen Breite des 
Lichtbogens. Unter ihr soll der Raum verstanden werden, in dem der Ladungstragertransport 
zwischen Katode und Anode stattfindet. Aus dem Bild 28 entnimmt man fiir den Lichtbogen 
eine elektrische Breite von etwa 4,5 cm. Sie ist identisch mit dem Abschnitt b des Photozellen- 
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impulses der 3. Photozelle (Bild 25 und Bild 29). Der Abschnitt ¢ rithrt nur noch von einer 
leuchtenden Gaswolke her, die nicht mehr am Leitungsmechanismus beteiligt ist. DaB die 
Photographie eine etwas zu kleine Lichtbogenbreite anzeigt, diirfte auf die zu harte Entwicklung 
von Film und Bild zuriickgefiihrt werden. Aus Photographie und Photozellenimpuls erkennt 
man ferner, da sich der vordere Rand des Lichtbogens scharf zur Umgebung abzeichnet, wahrend 
der hintere Rand des Bogens ungleichmafig und durch Nebel verwaschen erscheint. Ein scharfer 
Ubergang zwischen leitendem und nicht leitendem Plasma ist nicht vorhanden. Dementsprechend 
ist die positive Induktionsspitze nicht so scharf ausgepragt wie die negative. 


Laufschiene 


B=100%o 


B=20% 
T B=12,5% 
Coes eae REY 
b--=—c cou 
Bild 29. Impuls der 3. Photozelleam Hallgenerator. Bild 30. Skizzierter Verlauf des Laufschienenfeldes Bz = f(y) 
(Entnommen aus dem Oszillogramm.) bei mehreren parallel brennenden FuBpunkten. 


In Bild 28 ist das am Modellichtbogen aufgenommene und auf einen Strom von J = 1230 A 
umgerechnete Laufschienenfeld gestrichelt miteingezeichnet. Fiir beide Feldverlaufe fehlt aller- 
dings ein gemeinsamer Bezugspunkt auf der y-Achse. Um jedoch zwischen ihnen einen Vergleich 
zu erméglichen, wurden die Punkte, an denen die jeweiligen Induktionen gerade auf die Halfte 
ihres Wertes von y = oo gesunken sind, zusammengelegt. Das am Modellbogen aufgenommene 
Laufschienenfeld zeigt einen wesentlich steileren Anstieg als das Gesamtfeld am Lichtbogen. 
Wahrend bei der Aufnahme des Laufschienenfeldes 1m Modell der Strom konzentriert an einem 
Punkt aus der Laufschiene austrat, ist dies beim Lichtbogen nicht mehr der Fall. Es existieren, 
wie in [4] gezeigt ist, an der Anode mehrere parallel brennende LichtbogenfuBpunkte. Fiir diesen 
Fall kann man sich naherungsweise eine in Bild 30 wiedergegebene Strom- und Induktionsver- 
teilung vorstellen, die auch die unterschiedliche Steigung in der Verbindungslinie zwischen beiden 
Induktionsspitzen in Bild 28 erklart. Der Fu8punkt am vorderen Lichtbogenrand fiihrt den 
groBten Teilstrom. Von ihm rihrt daher auch der Teil des steilsten Induktionsanstieges her. 
Hinter ihm folgen FuBpunkte, die einen kleineren Strom fiihren. Die Teilstréme werden immer 
geringer, je naher der FuBpunkt am hinteren Rand des Bogens liegt. Fiir jeden Teilstrom kann 
_ man nun eine eigene Induktionsverteilung aufstellen. Der Verlauf der Gesamtinduktion ergibt 
sich dann aus der Uberlagerung der Teilinduktionen. Der zur positiven Induktionsspitze hin 
immer flacher ansteigende Induktionsverlauf hat seine Ursache in den zum Bogenende kleiner 
werdenden Teilstr6men. 

Die Stromdichteverteilung im Lichtbogen ist dem von der Photozelle angezeigten Verlauf der 
Leuchtintensitat ahnlich. Die Stromdichte erreicht nach sehr steilem Anstieg kurz hinter dem 
vorderen Lichtbogenrand den maximalen Wert, um dann langsam bis zum Lichtbogenende ab- 
zunehmen. Da sich der hintere Lichtbogenrand durch keine scharfe Grenze zwischen dem strom- 
fiihrenden Plasma und dem Neutralgas auszeichnet, wird die Stromdichte nicht von einem end- 
lichen Wert auf Null springen. Sie wird vielmehr einen stetigen, wenn auch steilen, Ubergang 
aufweisen. Der wandernde Lichtbogen besitzt also keinen rotationssymmetrischen Aufbau. 
Eine Bestimmung des Lichtbogeneigenfeldes aus der Differenz von Gesamtfeld und Laufschienen- 
feld ist nur dann méglich, wenn die genaue Stromverteilung auf die einzelnen FuBpunkte be- 


kannt ware. 
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Messungen und Berechnungen des Induktionsverlaufes am Lichtbogen nach dem hier ange- 
gebenen Verfahren sind nur dann sinnvoll, wenn der Bogen die ktirzeste Verbindung zwischen 
beiden Laufschienen bildet. Dies ist nicht immer der Fall. Die Form des wandernden Licht- 
bogens wird durch den Abstand und das Profil der Laufschienen sowie durch die Hohe des Stro- 
mes bestimmt. Fiir die obenangegebene Versuchsanordnung stellt cer Lichtbogen erst etwa 
oberhalb 600 A die kiirzeste Verbindung zwischen beiden Laufschienen dar. Bei kleineren Str6- 
men erfahrt der Lichtbogen eine z. T. recht erhebliche seitliche Auslenkung aus der Schienenebene. 
Verbunden ist damit eine starke Absenkung des Blasfeldes, das den Lichtbogen in Laufschienen- 
richtung antreibt. Bei dieser Lichtbogenform tritt ein vollig anderer Bewegungsmechanismus 
des Bogens auf. Eine Bestimmung des Blasfeldes ist dann nach der hier angegebenen Methode 
nicht mehr méglich. 


Zusammenfassung 


Die am wandernden Lichtbogen vorhandene Induktion wird in zwei Anteile aufgespalten, 
von denen der eine vom Lichtbogen selber, der andere von den stromdurchflossenen Lauf- 
schienen erzeugt wird. Nur der zweite Anteil bildet in Verbindung mit dem Strom die den Licht- 
bogen antreibende Blaskraft. Die Bestimmung des Laufschienenfeldes bei kleinen Schienen- 
abstanden kann nur unter Beriicksichtigung des Laufschienenprofils erfolgen. Da die exakte 
Berechnung des Laufschienenfeldes in der Umgebung des Lichtbogens schwierig ist, wird ein - 
Naherungsverfahren angegeben. In der Héhe der Lichtbogenachse betragt die Induktion gerade 
die Halfte von dem Wert, der vorhanden ware, wenn die Laufschienen von + co bis —oo vom 
Strom durchflossen wiirden. Der Verlauf des Schienenfeldes konnte in einer besonderen Versuchs- 
anordnung, in der das Lichtbogeneigenfeld ausgeschaltet wurde mit dem Hallgenerator gemessen 
und mit dem berechneten verglichen werden. Wahrend sich bei schmalen Profilen eine gute 
Ubereinstimmung zeigt, ist dies bei breiten Profilen infolge der in der Berechnung nicht beriick- 
sichtigten Stromlinienverteilung in FuBpunktnahe nicht mehr der Fall. 

Die rechnerische Addition von Laufschienenfeld und Lichtbogeneigenfeld steht im Einklang 
mit der mittels Hallgenerator an einem Modellichtbogen gemessenen Induktionsverteilung. 

Dagegen ergeben sich aus dem Vergleich der Induktionsverlaufe zwischen Modellichtbogen 
und wanderndem Lichtbogen grundsatzliche Unterschiede. Ihre Ursachen sind in der von der 
Zylinderform stark abweichenden Gestalt des laufenden Bogens mit ungleichmaBiger Strom- 
dichteverteilung und den zahlreichen parallel brennenden FuBpunkten zu suchen. Die GréBe des 
Laufschienenfeldes an der Stelle der héchsten Stromdichte ist etwa auf die Halfte des Wertes 
abgesunken, der sich fiir y = + oo ergibt, sofern der Lichtbogen die kiirzeste Verbindung zwi- 
schen den Laufschienen bildet. 

Die vorliegende Arbeit entstand im Zusammenhang mit Untersuchungen iiber den Bewegungs- 
mechanismus vom wandernden Gleichstromlichtbogen im Institut fiir Hochspannungs- und 
MeBtechnik der T. H. Darmstadt auf Anregung von Prof. Dr.-Ing. FrUHAUF. 
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Transformatoren und andere Vierpole in der Netzberechnung 
mit Inzidenzmatrizen 


Von 
Hans EpELMAnNN, Erlangen 


Mit 6 Textabbildungen 


(Eingegangen am 7. Mai 1960) 


Ubersicht: In den vorausgegangenen Arbeiten wurden nur Methoden zur Berechnung von Netzen, welche 
sich ausschlieBlich aus Zweipolen zusammensetzen, behandelt. Transformatoren kénnen zwar, sofern ihr 
Ubersetzungsverhaltnis reell ist und sie nicht gleichzeitig auch Isoliertransformatoren sind, durch eine Zwei- 
pol-Ersatzschaltung wiedergegeben werden. In allen anderen Fallen empfiehlt es sich, zur Aufstellung der Glei- 
chungssysteme fiir Strome und Spannungen entweder Kopplungselemente in die Zweig-Impedanz-Matrix ein- 
zufiihren oder die Strom- bzw. Spannungsiibersetzungsverhiltnisse durch Zwangsbedingungen zu_beriick- 
sichtigen. Im letzteren Fall werden die Inzidenzmatrizen modifiziert. Fiir alle erwahnten Methoden werden 
hinreichende (in der Praxis meist erfiillte oder leicht herstellbare) Bedingungen fiir die Auflésbarkeit der 
Gleichungssysteme angegeben. 


1, Das Reziprozitétsgesetz bei Zweipolnetzen und Transformatoren 


In vorausgegangenen Arbeitent wurden die Grundlagen der Netzberechnung mit Inzidenz- 
matrizen behandelt. Dort wurden Netze betrachtet, welche sich ausschlieBlich aus Zweipolen 
zusammensetzen (RLC-Netze). Die Maschen-Impedanzmatrix, wie auch die Knotenpunkts- 
und Schnittmengen-Admittanzmatrix solcher Netze sind symmetrisch (zur Hauptdiagonale). 
Fiir die Symmetrie dieser Matrizen ist also die Eigenschaft der Netze, ausschlieBlich aus Zwei- 
polen aufgebaut zu sein, hinreichend, jedoch nicht notwendig, wie weiter unten gezeigt wird, 
denn sobald man Vierpole hinzunimmt, welche selbst durch symmetrische Impedanz- oder Ad- 
mittanzmatrizen dargestellt werden, ergeben sich ebenfalls wieder symmetrische Impedanz- bzw. 
Admittanzmatrizen. Netze mit symmetrischen Impedanz- bzw. Admittanzmatrizen nennt man 
auch reziprok, da sie dem Reziprozitatssatz geniigen, welcher besagt, daB man bei Einspeisung 
in ein Netz am 7-ten Klemmenpaar, alle tibrigen Klemmenpaare gedffnet, mit 


4;=% und bei Messung von u, = up 


am k-ten Klemmenpaar bei einer zweiten Einspeisung am k-ten Klemmenpaar, alle tibrigen 
Klemmenpaare geoffnet, mit 
1, = 1% die gleiche Spannung namlich u; = u, 
am i-ten Klemmenpaar miBt; denn es ist im 1. Versuch 
Uy = %y 0; = Xi 0 

und im 2. Versuch 
j U, = 254 ty = 2p to» 
woraus sich aus der Gleichheit von z,; und z;, die Gleichheit von u, und uw, sofort ergibt. Diese 
Eigenschaft gilt in allen Netzen, in welchen die ohmsche Verluste und ebenso auch die magneti- 
sche und elektrische Energie ZustandsgréBen der Einspeisespannungen bzw. -strome sind.? 

Das Reziprozitatsgesetz gilt also immer dann, wenn die Elemente, aus denen das Netz auf- 


gebaut ist, selbst dem Reziprozitatsgesetz gehorchen, und zwar vollig unabhangig von Passivitats- 
eigenschaften; denn es gibt passive Netze® (z.B. gewisse gyratorische Vierpole), welche dem 


1 Allgemeine Grundlagen der Netzberechnung mit Inzidenzmatrizen“ [72], im folgenden kurz AGNI 
genannt, und ,,Die Inzidenzmatrizen der Schnittmengen eines Netzes‘‘ [82], im folgenden kurz ISNE genannt. 

2 Vel. die grundlegenden Gedanken im Empzschen Beweis fiir die Reziprozitat beim Transformator in 
H. G. M6LueEr [74]. 

3 Ein Netz soll passiv heiSen, wenn es bei Erregung mit Spannungen oder Strémen der Form Re c e? bei 
Re p = o im Mittel cine nicht negative Summenwirkleistung aufnimmt. 


33% 


~, 
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Reziprozitatsgesetz nicht gehorchen, umgekehrt gibt es auch aktive Netze (z.B. solche mit nega- 
tiven ohmschen, induktiven oder kapazitiven Elementen), welche dem Reziprozitatsgesetz ge- 
horchen. 

Ein Vierpol, der dem Reziprozitatsgesetz gehorcht, ist der Transformator, genauer gesagt, 
der Transformator mit reellem Ubersetzungsverhaltnis. Netze, welche auBer ohmschen Wider- 
standen, Kapazitaéten und Induktivitaéten noch Transformatoren enthalten, gehorchen somit, 
wie auch weiter unten gezeigt wird, dem Reziprozitatsgesetz. Da das Reziprozitatsgesetz hier 
gilt, wird man versuchen, Transformatoren durch Zweipole darzustellen. 


2. Zuriickfiihrung auf die Inzidenzmatrizenmethode 


2.1 Darstellung des Transformators durch eine Schaltung mit Zweipolen 
2.11 Der ideale Transformator 
Der ideale (d. h. der zugleich streuungs- und magnetisierungsstromfreie) Transformator mit 


reellem Ubersetzungsverhaltnis 1a4Bt sich, falls man auf galvanische Trennung verzichten kann, 
tatsachlich durch eine Kettenschaltung von Zweipolen darstellen, wodurch sich schlieBlich wieder 


L 2 ty 
au ahs (1-r9)r (aoe 4 


Bild 1. Eine aus Zweipolen zusammengesetzte Ersatzschaltung fiir einen idealen Transformator 1:w 
(untere Klemmen durchgeschaltet). 


ein Netz ergibt, das ausschlieBlich Zweipole enthalt. Allerdings sind jeweils zwei der Impedanzen 
der Ersatzschaltung negativ. Wollte man einen idealen Transformator durch positive RLC- 
Elemente nachbilden, so ist dies nur dadurch méglich, daB man die ,,positiven‘‘ Impedanzen 
durch Induktivitaten, die ,,negativen“ Impedanzen durch Kapazitaten darstellt (oder umgekehrt). 
Hierdurch wird der Transformator allerdings frequenzabhangig. In der Rechnung kann man 
natiirlich auch mit negativen ohmschen Widerstanden arbeiten. DaB die in Bild 1 dargestellte 
Ersatzschaltung einem idealen Transformator dquivalent ist, erkennt man am einfachsten 
folgendermaBen: Legt man an das linke Klemmenpaar eine Spannung , an, so ergibt sich an 
den offenen Klemmen rechts eine Spannung 


Wx 


fie Uy = WU. (1a) 


1—w)*¥+we 


Diese Spannung ist konstant und unabhangig vom Belastungsstrom 7,, denn der innere Wider- 
stand der gesamten Anordnung vom rechten Klemmenpaar aus betrachtet ist 


Ry = (w’— w) x +- 0. (2a) 


Es ist also allgemein 
Up = W Uy. (3a) 


Speist man links einen konstanten Strom 7, ein, so kann auch gezeigt werden, daB der iiber das 


rechte Klemmenpaar flieBende Strom unabhangig von dem Abschlu8widerstand ist, denn der 
KurzschluBstrom ist 


194, = —— ay ae = = aay (1b) 
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und der innere Widerstand dieser Stromquelle ist 


R, = (w® — w) x at A (68), (2b) 


wx (1—w)x—-* 


wodurch der Strom 7, unabhangig von Belastungswiderstand und damit von w, wird, d.h. es gilt 
allgemein 
1 


=U (3b) 


Ww 


— 


Die Gl. (3a) und (3b) sind aber die Gleichungen des idealen Transformators, wodurch gezeigt ist, 
daB die in Bild 1 wiedergegebene Schaltung einem idealen Transformator aquivalent ist. 


2.12 Der halb-ideale Transformator mit Streuungsimpedanz 


Jeden realen Transformator kann man nun aus einem idealen Transformator und aus den 
Streu- und Magnetisierungsimpedanzen zusammensetzen. Da man in Netzuntersuchungen der 
Energietibertragungstechnik gewohnlich die Magnetisierungsimpedanz unberiicksichtigt laBt, ist 


fi 4:7 Be 


Bild 2. Eine aus Zweipolen zusammengesetzte Ersatzschaltung fiir einen halb-idealen Transformator 1: w mit Streuung 
(untere Klemmen durchgeschaltet). 


hauptsachlich der Sonderfall ,,Transforrmator mit Streurmpedanz“ wichtig. Hierbei gelingt es, 
durch Kombination einer in der Ersatzschaltung des idealen Transformators enthaltenen nega- 
tiven Impedanz mit der Streuimpedanz eine Vereinfachung zu erzielen. Man ist némlich in der 
Wahl der GréBe x véllig frei. Aus diesem Grunde ist es zweckmaBig (w?— w) x gleich dem 
negativen Wert der auf die rechte Seite des Transformators mit Streuung (nach Bild 2) be- 
zogenen Streuimpedanz zu setzen, so daB diese sich zu null aufheben und der rechte Langszweig 
fortfallt. Damit ergibt sich die Forderung 


(w2 — w) x = — (w?z,+ 2), 
so daB man erhalt 
W 24 + fa 
a. Oe ok; 
< w—t1 =a (4) 


Hierbei wurde zur Abkiirzung die mittlere Streiumpedanz 
%44 = wat = (5) 


eingefiihrt. Somit ergibt sich die in Bild 2 rechts wiedergegebene Ersatzschaltung, die merk- 
wiirdigerweise einfacher ist als diejenige fiir den idealen Transformator. Der Sonderfall des 1:1- 
Ubertragers ist in dieser Schaltung natiirlich auch enthalten. In diesem Fall entfallen die Quet- 
impedanzen und es verbleibt die Streuimpedanz in Gestalt von z. Da8 die Schaltung ftir den 
idealen Transformator komplizierter, d.h. um eine Impedanz reicher sein muB, licgt offenbar 
daran, daB die gesuchte Ersatzschaltung folgende Bedingungen erfiillen muB: 


1. soll die Ersatzschaltung beim Leerlauf auf der einen Seite die Eingangsimpedanz fiir die 
andere Seite oo ergeben, 
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2. soll die Ersatzschaltung bei KurzschluB auf der einen Seite die Eingangsimpedanz fiir die 


andere Seite 0 ergeben, 
3. die Schaltung soll aus endlichen und von null verschiedenen Widerstanden aufgebaut sein. 


Insbesondere die 3. Forderung ist hier auch wesentlich, denn versucht man nach Bild 2 den 
Grenziibergang |z,| +0 zu machen, so entartet die Schaltung. 

Auch hier kommt man nicht ohne negative Impedanz aus. Vom Grenzfall des 1: 1-Transformators 
abgesehen, ist immer entweder die linke oder rechte Querimpedanz negativ. Mit dem Vorhanden- 
sein negativer Impedanzen taucht aber auch die Frage auf, ob in dem Netz Kreisstréme flieBen 
k6énnen, wodurch sich die bereits in AGNI erwahnten Auflésungsschwierigkeiten der Gleichungen 
ergeben kénnen. Hierzu ist aber folgendes zu sagen. Die Ersatzschaltung des Transformators 
fiir sich, d. h. ohne angeschlossene Netzteile, also auf beiden Seiten leerlaufend, 14Bt einen Kreis- 
strom in der Masche — x, (1 — w) x, w x nach Bild 1 zu, da die Summe der Maschenimpedanzen 
null ergibt.4 Sobald aber wenigstens auf die Primar- oder Sekundarseite ein Netzteil mit verlust- 
behafteter Eingangsimpedanz angeschlossen ist, ergeben sich insgesamt positive Verluste. Aus 
energetischen Griinden kann ein solcher Strom nicht dauernd flieBen. 

Damit ist der Fall des galvanisch nicht getrennten Transformators (d.h. untere Klemmen 
durchgeschaltet) auf ein reines Zweipolnetz zuriickgefiihrt, dassich unter denin AGNI angegebenen 
hinreichenden Bedingungen ohne Schwierigkeiten rechnen 1a8t, da die Determinante des ent- 
stehenden Gleichungssystems nicht verschwindet. 

Es seinoch erwahnt, da8 der halb-ideale Transformator mit Magnetisierungsimpedanz 
dual hierzu in Form eines T-Gliedes dargestellt werden kann. 


2.2 Darstellung des Transformators durch Einfiihrung von Kopplungs-Koeffi- 
zienten in der Zweigimpedanz- bzw. Admittanzmatrix 


2.21 Die Maschenmethode 


In Gl. (4) von AGNI wurden die Zweigimpedanzen 2; ; als Diagonalelemente der Zweigimpedanz- 
matrix eingefiihrt, da die Zweigspannungen nur von den Zweigstrémen abhangen. Liegt nun 
zwischen zwei Zweigen ein Transformator, so ergeben sich auf Grund der magnetischen Kopplung 
auch Glieder auBerhalb der Hauptdiagonale. Auf diese Weise lassen sich also nach der Maschen- 
methode auch Transformatoren (oder sonstige Vierpole) einfiithren, sofern deren Impedanz- 
matrix existiert. Natiirlich kann man auf diese Weise auch andere Vierpole (z. B. Leitungsstiicke) 


Bild 3. Transformator 1: w mit primar- und sekundarseitiger Streuimpedanz 2,0 und zgq 


und mittlerer Hauptimpedanz zp (auf die Primarseite bezogen = ee , auf die Sekundarseite bezogen = |w| zp) 


|| 


einfithren, Nun kann es zwar vorkommen, daB die Impedanzmatrix eines solchen Vierpols nicht 
existiert, dies ist z.B. beim idealen Transformator der Fall, aber auch beim Transformator ohne 
Magnetisierungsstrom. Diese Schwierigkeit laBt sich am leichtesten dadurch umgehen, daB man 
eine von null verschiedene Magnetisierungs-Impedanz in der GroBenordnung der in dem betrach- 
teten Transformator wirklich auftretenden Magnetisierungsimpedanz einfiihrt. Diesen prakti- 
schen Weg sollte man tiberhaupt in allen Fallen gehen, in welchen sich infolge eines exakten Null- 


“ 


4 Wegen dieser Eigenschaft wird diese Transformator-Ersatzschaltung auch gern Resonanz-JJ- oder auch 
Resonanz-A-Schaltung (engl.: resonant-4) genannt. Vgl. E. CLARKE [73], S. 174. 
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setzens einzelner GréBen Anderungen im Rechengang und somit lastige Fallunterscheidungen 
ergeben. Aus diesem Grunde soll hier zunachst der allgemeine Transformator mit Streuung und 
Magnetisierungsstrom eingefiihrt werden. 


Bild 3 zeigt das Ersatzschaltbild des allgemeinen (Zweiwicklungs-) Transformators. Dort 
kann man die folgenden Beziehungen ablesen 
WwW 


* : w i ; 
[w| En, Ye > UL S [| 41 + (26 + % ||) 29 - (6) 


Ue (:1 zt ce tae 
Hierbei wurde, um auch phasendrehende Transformatoren behandeln zu kénnen, der allgemeinere 
Transformator mit komplexem Ubersetzungsverhaltnis in Betracht gezogen5. Der Transfor- 
mator mit reellem Ubersetzungsverhaltnis ist als Sonderfall in diesem allgemeineren Fall ent- 
halten (natiirlich auch mit negativem Ubersetzungsverhaltnis, Klemmenvertauschung!). Das 
Gleichungssystem (6) in Matrizen geschrieben lautet 


Ep Ww 
6 i ae —-%, 
Uy ‘ [| Je] aa Cat Se 1 : 
epee w SOE ae 3 ete Le (7) 
? eh Zag + % || ? ‘ ae = 


|| 


Wir haben damit die Z-Matrix des Transformators mit Streuung und Magnetisierungsstrom und 
fiir ein komplexes Ubersetzungsverhaltnis gewonnen. Die Matrix ist fiir nicht reelle Uber- 
setzungsverhaltnisse w nicht symmetrisch. Mit |z,| +0 erhalt man auch den Sonderfall der 
Entkopplung. Als Diagonalelemente verbleiben die Streuimpedanzen z,, und 2,, die jetzt die 
Rollen von Zweigimpedanzen spielen. 


272) 223, 432 233 


Generator Last Transformator Leitung Last 


Bild 4. Berticksichtigung eines Transformators in der einphasigen Darstellung eines Energictibertragungssystems bestehend 

aus Generator, Eigenverbrauchslast, Transformator, Leitung und Hauptlast durch Einftihrung von Kopplungskoeffizienten. 

a Originalschaltbild, b Zugehériger orientierter Graph mit vollstandigem Baum (dick ausgezogen) und System unabhangiger 
Zweige (diinn ausgezogen). 


Wir betrachten hierzu ein Beispiel. Bild 4 zeigt ein einfaches kettenférmiges Netz eines 
Energieiibertragungssystems, das einen Transformator enthalt mit einem Ersatzschaltbild nach 
Bild 3. Der Transformator erzeugt eine Kopplung der Zweige 2 und 3. Die Erdriickleitung soll 


5 Vgl. zur Theorie des idealen Transformators mit komplexem Ubersetzungsverhialtnis [77]. 
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als ein einziger Knotenpunkt aufgefaBt werden. Die zugehdrige Zweig-Zyklen-Matrix hat 


dann folgende Gestalt: 


1 ) fo) fo) O 
) 1 ) fo) O 
Oo fo) 1 fo) ) 
es fo) ) ) 1 ) (8) 
) fo) ) ) 1 
re = O 0 fo) 
fe) OO i Siu eens 
) 0 fo) 1 —1 


Eliminiert man die induktive Kopplung der Zweige 2 u 3, so sind die verbleibenden Teilnetze 
entkoppelt. Aus diesem Grund erscheint bei entsprechender Umordnung und zwar spalten- und 
zeilenweise in der Reihenfolge 1, 2, 6; 3, 4, 5, 7, 8 die Matrix C in der Gestalt 


0 Coo 


Da auch die Maschen entkoppelt sind, ergibt sich bei entsprechender Partitionierung auch im 
unteren Teil der C-Matrix eine Kastchenbildung der gleichen Art, wie es in (8) durch gestrichelte 
Linien angedeutet ist. Die verallgemeinerte Zweigimpedanzmatrix weicht jetzt also 
von der Diagonalform ab und enthalt auBerhalb der Hauptdiagonale die Kopplungselemente 293 
und 239, die sich analog zu Gl. (7) errechnen. Auf diese Weise hat die Zweigimpedanzmatrix 
folgende Gestalt 


LO (0 50. Os LOles OmmEO 
On 253 Zag) (Oe O0e, Om O ae FO 
On ahh cag On a1 Oy AO ROE 
5 = |S t0 Oat OO” Oma (0) 
0) 0% =O) © Oi, 29, BOr e70" meO 
0. VOe OF Oy Lome en) oo 
On" "0. “Os 50'S Oe 50 F ceO 
| Oo -'O1 Os 0) 20 FO On era. 


Somit kénnen die Gleichungen von AGNI iibernommen werden Der einzige Unterschied besteht 
darin, da8 Z nun auch auBerhalb der Hauptdiagonale von null verschiedene Elemente besitzt. 
Mierdurch werden natiirlich die Regeln fiir den Aufbau der Maschenimpedanzmatrix Z° kompli- 
zierter. Im einzelnen gilt aber natiirlich auch hier 


Ua Zia. ec, (10) 
wobei in e nur ¢, = 0, wahrend ¢,=\...== és = 0 Ist. Mut 
i=Ci?, (11) 
u° = C,u (12) 
und 
Ur ==0) (13) 


erhalt man das fiir die Maschenstréme giiltige Gleichungssystem 
U = 292 SS Cre 
oder 
Zoi? = — Cre’, (14) 
wobei ; 


Lo =/COZ GS (15) 
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Setzen wir die Matrizen nach den Gln. (8) und (9) in Gl. (15) ein, so ergibt sich fiir die Maschen- 
impedanzmatrix des Beispiels nach Bild 4 


£11 + 266 266 e) o) ce) 
“66 Zo + 266 493 Le) e) 
Zo = O 230 233 1 2qq 2q9 0 . (16) 
O 0 2qq 24g 2a + See 288 
O (e) ) 2g 255 + 2g 


Es fallt auf, daB oberhalb und unterhalb der Hauptdiagonale nur je eine Reihe besetzt ist, was 
typisch ist fiir Kettenschaltungen. Die Matrix Z° ist nur dann symmetrisch, wenn die auf diese 
Weise eingefiihrten Vierpole ebenfalls symmetrische Impedanzmatrizen haben. Fiir die Elemente 
dieser Matrix gelten in Erweiterung der entsprechenden Regeln in AGNI (S. 422 unten) die folgen- 
den Regeln: 

1. In der Hauptdiagonale ist das Element z?; gleich der Summe aller Impedanzen, welche der 
a-ten Masche angehdren, vermehrt um die algebraische Summe derjenigen Zweigkopplungs- 
_ impedanzen z,,,, welche sowohl mit dem u-ten, als auch mit dem »-ten Zweig in der 7-ten Masche 
enthalten sind; das Vorzeichen dieser Zweigkopplungsimpedanzen ist positiv oder negativ je 
nachdem ob, die Orientierungen der gekoppelten u-ten und v-ten Zweige innerhalb der 7-ten 
Masche gleich oder entgegengesetzt gerichtet sind. 


| 2. AuBerhalb der Hauptdiagonale besteht das Element z?, aus zwei Bestandteilen. Der eine 
_ Bestandteil ist gleich der algebraischen Summe aller Baumimpedanzen, welche gleichzeitig der 
a-ten und k-ten Masche angehéren, wobei das Vorzeichen einer solchen Baumimpedanz in dieser 
Summe positiv oder negativ ist, je nachdem, ob die Maschen gleichsinnig oder gegensinnig tiber 
die betreffende Baumimpedanz laufen. Der andere Bestandteil ist die algebraische Summe der- 
jenigen Zweigkopplungsimpedanzen -+ z,,, welche sowohl mit dem y-ten Zweig einerseits als 
auch mit dem »-ten Zweig andererseits in der i-ten und k-ten Masche enthalten ist, wobei das Vor- 
zeichen der Zweigkopplungsimpedanzen z,,, positiv ist, wenn die Orientierungen des mw-ten und 
y-ten Zweigs mit denjenigen der 7-ten und k-ten Masche beiderseits gleich- oder beiderseits gegen- 
sinnig zueinander sind und anderenfalls negativ ist. 


Der Beweis fiir diese Regeln ergibt sich dadurch, daB man das Element der 7-ten Zeile und 
k-ten Spalte in der Maschenimpedanzmatrix betrachtet. Dieses Element ist gleich 


z &z 
a Ce 0, ee 
M=1v=1 
In dieser Doppelsumme ergibt sich nur dann ein Beitrag von + z,,, wenn bei festgehaltenem 
und k, d. h. in der 7-ten Spalte und k-ten Spalte von in der w-ten und »-ten Zeile gleichzeitig von 
null verschiedene Elemente + 1 vorhanden sind. Aus der Kombination der Vorzeichen von c¢,,, 
und ¢,, ergeben sich die Vorzeichenregeln. 


Auch hier erhalt man die Maschenstréme eindeutig, falls die Determinante des Systems 
nicht verschwindet, aus dem Gleichungssystem mit y Unbekannten nach Gl. (14) 


Ze —= Gren (14) 


Alsdann erhalt man die Zweigstréme aus Gl. (11) und die Zweigspannungen aus Gl. (10), womit 
alle wesentlichen GréBen des Netzes ermittelt sind. 


2.22 Die Knotenpunktsmethode 


Natiirlich 148t sich auch die Knotenpunktsmethode so erweitern, daB Transformatoren 
(oder sonstige Vierpole) eingefiithrt werden kénnen. Hierzu sind nur in der Zweigadmittanz- 
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matrix die entsprechenden Kopplungselemente y,, zu besetzen. Die Admittanzmatrix des Trans- 
formators nach Bild 4 erhalt man durch Inversion der Impedanzmatrix. Es gilt dann mit 


w* 
; aq + %, [72 Toy 5 
vA 1 ' 
t= ta el ty 2 (‘\== —~— ; (8). (17) 
yo ee iy ea 
|20| a “16 ais |x| 


Im Gegensatz zur Maschenmethode kann man nun hier den Grenziibergang |z,|-> 00 ausfiithren, 
ohne daB Elemente der Matrix unendlich werden. Allerdings wird die Matrix dann singular, 
denn ihre Inverse existiert ja nicht. Fiir den halb-idealen Transformator ohne Magnetisierungs- 
strom gilt dann: 


Sl ie 
( “a Ze || 20 ta (18) 
(2) = w 1 Us 

|w or [| 26 


Wenden wir uns wieder unserem Beispiel nach Bild 4 zu. Fiihrt man die Admittanzmatrix 
des Transformators in die Zweig-Admittanzmatrix ein, so ist die Knotenpunkts-Zweig-Matrix 
mit der dort angegebenen Knotenpunktsnummerierung (die gesamte Erdriickleitung wieder als - 
ein Knotenpunkt Nr. 1 aufgefaBt) 


1 1 O 1 at 1 1 
Sih O O Os ail O 10) 
Kees uae (19) 
O O 1 1 O ler al O 
O 10) 0 ib =i 10) OF. 


und die Zweigadmittanzmatrix unter sinngemaBer Anwendung von Gl. (18) bzw. (19) 


{ vaa ) ) ) ) ) 0) ) 
0 BAYS BLY ee 0 0 O O 
0 Yao 7 338012 0 0 0 O 

ae ) ) ) Veq 0 ) ) ) as 

) fe) ) ) Neg sO re) ) 
) ) 0 Suko ) Vee~ 0 ) 
) ) ) ) ) ) eR ae 
) ) ) fe) ) ) 0: ee Vee 


Auch hier kénnen jetzt die Gleichungen von AGNI iitbernommen werden. Fiir die Zweigstr6me 
gilt ebenso 


i= Y (u—e), (21) 
wobei wieder in @ zu setzen ist ¢; += 01nd ¢, —=.4....—= yea O's Mt 
feet GI (22) 
ee ate . (23) 
und 
poss (24) 
erhalt man das Gleichungssystem fiir die Knotenpunktsspannungen 
it = Yn Keyieme0 (25) 
oder 
V2 ere (26) 
wobei 


Ye oS Key Ke (27) 
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Durch Einsetzen von (20) und (21) und Ausmultiplizieren erhalt man 


r Yi + Yoo + Vso 4 A ee ON iat 
+ Yes + Y33 + V55 7 Sa eile ae cae ss = SA33 == SE 
Bega yen toy. Ys2— Yee Y77 
Digi Bop RES ONE, San y A 
) a maar OS aa) 6G + Yee 7 sa (20) 
V3 ane 38 Vap Vee toy 4a ot ey 
4 
AG + Vaq ‘ 
me.5 Siam 88 ce) saa AA gree Co 
+ Yee J 


Wegen der besonderen Eigenschaft der Matrix K, daB jede Zeile die negative Summe aller iibrigen 

Zeilen ist, ist Yt singular, d.h. es ist auch in Y* jede Zeile oder Spalte die negative Summe der 
Uy 

itbrigen. Damit kénnen die Elemente von ut =| : | nicht eindeutig, sondern nur bis auf eine 
Uy 

willkiirliche Konstante bestimmt werden. Jede (k —1)-reihige Hauptuntermatrix ist i. allg. 

Rangmatrix (wenn man Resonanzfille ausschlieBt). Fiihren wir z.B. den Knotenpunkt 1 als 

Bezugsknotenpunkt ein und setzen uj = 0, so erhalt man dadurch ein im allgemeinen auflés- 

bares Gleichungssystem fiir die Knotenpunktsspannungen, da8 man in K die 1. Zeile streicht, 

diese Matrix heiBe K, und in Y* die 1. Zeile und die 1. Spalte streicht, diese Matrix heiBe V3 

und anstelle von (26) einfiihrt® 


Y+t+=KYe. (26a) 
Hieraus erhalt man die Zweispannungen explizit aus Gl. (24) und die Zweigstréme explizit aus 
Gl. (22). 
Die kombinatorischen Regeln zur Aufstellung von Y* lauten folgendermaBen: 


1. In der Hauptdiagonale ist das Element y,; gleich der Summe aller Zweigadmittanzen, 
welche mit dem 7-ten Knotenpunkt verbunden sind, vermehrt um die algebraische Summe der- 


| jenigen Zweigkopplungsadmittanzen + y,,, deren u-ter und y-ter Zweig mit dem 7-ten Knoten- 
| punkt verbunden ist, wobei das Vorzeichen von Vx» positiv oder negativ ist, je nachdem die 


Definitionspfeile des w-ten und »-ten Zweigs gleich- oder gegensinnig mit dem 7-ten Knotenpunkt 


verbunden sind. 


2. Das auBerhalb der Hauptdiagonale stehende Element y,, ist gleich der negativen Summe 
aller Zweigadmittanzen, welche den 7-ten und k-ten Knotenpunkt miteinander verbinden, ver- 
mindert um die algebraische Summe derjenigen Zweigkopplungsadmittanzen -+ y,,, deren 
p-ter Zweig einerseits und y-ter Zweig andererseits mit dem 7-ten und k-ten Knotenpunkt ver- 
bunden sind, wobei das Vorzeichen von y,,, positiv ist, wenn die Definitionspfeile des y-ten und 
y-ten Zweigs beide zu den Knotenpunkten 7 und & hinweisen oder beide wegweisen und negativ 
ist, wenn diese Definitionspfeile beieinem der Knotenpunkte 7 und & hinweisen und beim anderen 
wegweisen. 


Der Beweis fiir diese Regeln kann analog wie in Kap. 2.21 gefiihrt werden, inden man 


Vik aa pa, Ds, Fin Vuv Ray 
p=H1vy=1 


betrachtet. 


6 Vgl. AGNI, Kap. 3. 
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2.23 Verallgemeimerungen 


Die Einfithrung von Transformatoren in die Zweipolnetze durch Besetzung der Kopplungs- 
elemente in den Zweigimpedanz- bzw. Admittanzmatrizen kann natiirlich auch auf andere Vier- 
pole ausgedehnt werden. Es sei in diesem Zusammenhang erwahnt, daB die Einfiihrung von 
Vierpolen in die Zweig-Admittanz-Matrix um die Knotenpunktsmethode zu verallgemeinern, 
identisch ist mit der Methode der Vierpolnetze von U. Hater [75]. Eine Ausdehnung auf Netze 
mit Vierpolen in allgemeiner Lage und ihre digitale Berechnung findet man bei G. NEEs [81]. 

Natiirlich kann man auf diese Weise auch andere 2 1-Pole (m > 2) einfiihren, sofern die Impe- 
danz- bzw. Admittanz-Matrix existiert (insbesondere auch Mehrwicklungstransformatoren). In 
all diesen Fallen wie auch bei Vierpolen mu8 aber ausdriicklich betont werden, daB diese Be- 
ziehungen zundchst nur giiltig sind, wenn den 2 m-Polen (m > 2) an »—1 Klemmenpaaren 
Isoliertransformatoren nachgeschaltet sind. Im Falle einer galvanischen Durchverbindung z. B. 
auch bei Spartransformatoren treten noch zusatzlich Bedingungen auf, z. B. Gleichheit der 
Potentiale an den galvanischen Verbindungen zweier Klemmen, die in 2 n-Pol-Matrix nicht zum 
Ausdruck gebracht sind. Es kann nun aber sein, da8 in der die 2 7 Pole umgebenden Schaltung 
diese Verbindungen sowieso schon enthalten sind; dann braucht man sie natiirlich nicht zu beriick- 
sichtigen. 


2.24 Eindeutige Auflésbarkeit der Maschen- und Knotenpunkts- und Schnittmengengleichungen ber 
Beriicksichtigung von Kopplungen 


Die in Kap.5 von AGNI ausgesprochenen hinreichenden Satze iiber die eindeutige Auflésbar- 
keit der Maschen- und Knotenpunktsgleichungen kénnen leicht auf den Fall einbezogener Kopp- 
lungen verallgemeinert werden. Es hat dann z. B. in dem Beweis in Gl. (63) von AGNI anstelle 
der reinen Quadratsumme mit positiven Koeffizienten im 1. Term eine positivdefinite quadra- 
tische Form zu treten, die man erhalt, wenn man fordert, daB der hermitesche Teil irgendeiner 
im System unabhangiger Zweige zugeordneten Untermatrix der Zweigimpedanzmatrix positiv- 
definit ist, das tibrige Netz aber keine aktiven Teile enthalt. Damit kénnen wir die erweiterten 
Satze folgendermaBen formulieren: 


A. Damit das System der Maschengleichungen 


PAN Ora (14) 


| eindeutige Lésungen hat, ist hinreichend, da8 


1. der hermitesche Teil der zugrundegelegten Zweigimpedanzmatrix positiv-semidefinit ist 


_ (bei passiven Netzwerken automatisch erfiillt), 


2. der hermitesche Teil irgendeiner einem System unabhangiger Zweige zugeordneten 
Untermatrix positiv-definit und damit Rangmatrix ist. 


B. Damit das System der Knotenpunktsgleichungen 


Y+ Ut =K,Ye (26a) 


| eindeutige Lésungen hat, ist hinreichend, daB 


1. der hermitesche Teil der zugrundegelegten Zweig-Admittanz-Matrix positiv-semidefinit 


_ ist (bei passiven Netzwerken automatisch erfiillt), 


2. der hermitesche Teil irgendeiner Baumuntermatrix Rangmatrix und damit positiv-definit 


| ist. 


Es soll hier nur noch erwahnt werden, daB man zum Beweise des Satzes B dual zum Satz A 
vorgehen mu8 und der in AGNI, Kap. 5.2 angegebene Beweis, der die Kombination der Baum- 
zweige verwendet, sich nicht in dieser Form erweitern 1aBt. Das entspricht aber einem Beweis 
fiir die Existenz der Lésungen der Schnittmengen-Gleichungen 


¥? i? = Feaer (29) 


Wegen det Y* = det Y+sind die beiden Existenzsatze aquivalent. 
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3. Beriicksichtigung halbidealer Transformatoren 
durch Einfithrung von Zwangsbedingungen. Modifizierte Inzidenzmatrizen 


Die bisher besprochenen Verfahren zur Einbeziehung von Transformatoren zeichnen sich da- 
durch aus, daB die Methode der Inzidenzmatrizen ohne besondere Modifikation tibernommen 
werden kann. 

Bei der Darstellung durch Zweipolnetze ist es zweckmaBig, die Impedanz x so fest- 
zulegen, daB eine Zweigimpedanz oder -admittanz verschwindet. Hierdurch ergeben sich ge- 
wisse Sondervorschriften fiir die Festlegung der Impedanzen der Ersatzschaltung des Trans- 
formators. Grundsatzlich sind alle Sonderfalle halb-idealer und idealer Transformatoren zulassig. 
Transformatoren mit komplexem Ubersetzungsverhaltnis kénnen nach diescr Methode nicht be- 
handelt werden. 

Bei der Darstelleng durch Einfithrung von Kopplungskoeffizienten bleibt der 
Graph unverdndert, jedoch kann es jetzt vorkommen, da8 Z oder Y nicht existieren oder weder 
Z noch Y existiert. Auf den letzteren Fall ist also streng genommen weder die Maschen- noch die 
Knotenpunktsmethode (bzw. Schnittmengenmethode) anwendbar, doch lassen sich diese theore- 


tischen Hindernisse durch folgenden praktischen Weg beseitigen: Je nachdem, ob oder |2,,| 


Zn 
oder/und |z,,| gleich null zu setzen sind oder auch alle genannten GréBen, fithrt man stattdessen 


hinreichend kleine Admittanzen bzw. Impedanzen ein’. Auf diese Weise enthebt man sich solchen 
Sonderbehandlungen. AuBerdem wird der Fall des idealen Transformators nach dieser Methode 
iiberhaupt erst rechenbar. In Tabelle 1 sind links die Graphen des allgemeinen Transformators 
und im Anschlu8 daran die in der entsprechenden Darstellung des allgemeinen Transformators 
und der Sonderfalle zusatzlich hinzukommende Anzahl von Zweigen Az, von Knotenpunkten 
Ak und von Maschen Ar zusammengestellt : 


Tabelle 1. 1. Darstellung durch Zweipolnetze (Kap. 2.1) 
Graph des allgem. halb-idealer Transformator nclealen 
Transtoumators Allgem. Transformator 2)| > 00 |2o,| %q,| — 00 Transformator 
Az = 5(4) 3 | 3 4 
Alje == 21) 1 fe) | 1 
Ar — 3(3) 2 3 | 3 
Pe ts OE = : is 
Bemerkungen: durch entsprechende | bei entsprechender | bei entsprechender 
Festlegung von * kann | Festlegung vonw | Festlegung von + 
| ein Zweig und ein Kno- | 
| tenpunkt eliminiert wer- 
den, dann gelten die ein- 
geklammerten Werte 
2. Darstellung durch Einfiihrung von Kopplungskoeffizienten (Kap. 2.2). 
F nee : 
Graph des allgem. Allgem. halb-idealer Transformator Fdeales 
Transformators Transformator 74 — 00 | lZo,|, 2g,| > 0 Transformator 
| 2 
Nl 
NO re) Oo == 
Ak =o 2 2B — 
Vie = PD 2 2 — 
Bemerkungen: | — Z existiert nicht. Nur | Y existiert nicht. Nur die | weder Z noch 
die Knotenpunkts-u. | Maschenmethode ist an- | Y existiert 
Schnittmengenmetho- | wendbar. | 
| de ist anwendbar 


7 z. B. 4 = 1 MQ, go, bzw. 2c, = 10772. Die Admittanzen 
gewahlt werden, weil sie sonst bei Rundungen gegentiber den 
einer echten Nullsetzung gleich. 


bzw. Impedanzen dirfen auch nicht zu klein 
iibrigen GréBen verlorengehen, dies kame aber 
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3. Beriicksichtigung durch Einfithrung von Zwangsbedingungen (Kap. 3): 
Graph des allgem Allgem halb-idealer Transformator idealer 
; s allgem. J d 
Transformators Transformator |7a] —> 00 Fo, | . |2o2| — oo Transformator 
S| SS ee EE ees 
AV si ==s 2 _— 
Ane — fo) — 
Ale =P a 2 ie 
Bemerkungen: — — nur die Maschenmethode = 
ist anwendbar 
Graph des halb-idealen | Allgem. halb-idealer Transformator idealer 
Transformators Transformator ie | eo) | 2a | Zg,| > © Transformator 
mit Streuimpedanz 1|> il> | 
Az = — 2 — = 
| | 
Ak = — | fo) | — _ 
| Ar = — | 2 | = oe 
e i ee a ae <8 
Bemerkungen: — | nur die Knotenpunkts- = = 
und Schnittmengen- 
methode ist anwendbar 


In den Bemerkungen sind Hinweise dariiber, welche Regeln zu beachten sind. Hierbei fallt 
auf, da8 durch die Einfithrung von Kopplungskoeffizienten der geringste Aufwand an zusatz- 
lichen Elementen im Graphen entsteht, was sich im Falle einer erforderlichen Elimination oder 
Inversion besonders giinstig auswirkt, weil hierdurch die Ordnung der Matrix nicht so groB wird. 
Doch muB hierzu gesagt werden, daB diese Vereinfachung auf Kosten der Einfachheit der Z- bzw. 
Y-Matrix geht. Diese Matrizen haben jetzt namlich nicht mehr Diagonalgestalt. Die Speicherung 
in einem Digitalrechner ist damit nicht mehr in Form eines Vektors moglich; die Speicherung 
als Matrix ware aber im Hinblick auf die vielen Nullen auBerhalb der Hauptdiagonale unratio- 
nell, so daB man zu Sonderlésungen dieses Speicherproblems greifen mu8 auf Kosten der Eleganz. 


3.1 Zwangsbedingungen in der Maschenmethode 


Will man auf die Behandlung von Transformatoren mit komplexem Ubersetzugsverhaltnis 
nicht verzichten, dabei aber gleichzeitig die Schwierigkeiten in der Speicherung der Zweigimpe- 
danzmatrix mit Kopplungselementen umgehen, so empfiehlt sich die Einfiihrung von Zwangs- 
bedingungen. 


Hierzu betrachten wir wieder das Netz nach Bild 4. Der in diesem Netz vorkommende Trans- 
formator mit dem Ubersetzungsverhaltnis 1 : 7,8) Soll auf der ,,1‘-Seite im Ersatzschaltbild 
eine Magnetisierungsimpedanz von %7,3, haben. Dadurch ist der Transformator im Ersatz- 
schaltbild zwar jetzt nicht mehr symmetrisch, es erscheint jedoch zweckmaBig, nicht nur das 
Ubersetzungsverhiltnis anstelle von w,:w, in der Form 1:w anzugeben, sondern auch dabei 
gleichzeitig die Magnetisierungsimpedanz grundsatzlich auf die ,,1‘-Seite zu verlegen. Der 
Sonderfall verschwindender Streuimpedanzen kann hiermit durch Weglassung der entsprechen- 
den Zweige (in Bild 5 2. und 233) behandelt werden. Der halb-ideale Transformator mit Magneti- 
slerungsimpedanz wird jetzt einfach durch die beiden Zweige 7 und 8 reprasentiert. Hierbei 
kann man, wie in Bild 5 geschehen, alle Querzweige, sofern sie nicht Parallelzweige sind (wie 1 
und 5), als Baumzweige einfithren. Wir beginnen wie in AGNI mit der Numerierung der Zweige 
bei den unabhangigen Zweigen. 


Der dem Netz zugeordnete Graph nach Bild 5b besteht zunachst noch aus zwei gegenseitig 
entkoppelten Teilen. Die Kopplung ergibt sich erst durch die nachfolgend eingefiihrten Uber- 
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setzungsverhaltnisse. Um dies in Form einer Matrizen-Transformation schreiben zu kénnen, 
soll jetzt neben dem Vektor der Zweigspannungen (u,.. . u9) ein transformierter Vektor von 
Spannungen (u,...,, u{, u49) eingefiihrt werden. Die Spannungen selbst sind die gleichen, es 
fehlt aber im Vektor der gestrichenen GréBen wu. Das entspricht einer Elimination dieser Span- 
nung aufgrund der Zwangsbedingung. 


Generator 


Bild 5. Beriicksichtigung eines Transformators in der einphasigen Darstellung eines Energictibertragungssystems entsprechend Bild 4 
durch Hinfiihrung von Zwangsbedingungen (Maschenmethode). a Originalschaltbild, b zugehériger orientierter Graph mit voll- 
standigem Baum und System unabhangiger Zweige (die Streuimpedanzen erscheinen hier als Zweige!). 


Wir erhalten dann 
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= = 

= 

o Sa Sa 
= 
~ 

= 

~ 


= 
SS 


Ss 
oO 

I 

& 
© 


im Matrizen 
Uns VVEU on (30) 


Fir die inversen Gleichungen der Stréme existiert ein Gleichungssystem mit der konjugiert- 
transponierten Matrix W# 
4 = 4 


Kah ee 


I 


I 
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in Matrizen 


iv = Wei. (31) 
Hier ist also 
+4 0 fo) ) ro) fo) 0 | 
o +1 O O O O 0) O 
0) o +41 (0) O O O 0) 
O (0) Oo -+1 (6) O O (e) 
W = 0 0 0 o +1 O 0 0) (32) 
O (e) (0) (0) Oo +1 O O 
O (0) ) 0 O 7,8) 0 0) 
O O O O O oO +1 O 
| o re) ) 0 0 ) Gretsi 


Der transformierte Strom 7, ist hier nicht identisch mit 7,, er ist vielmehr der Strom, der durch 
27, flieBt, also der auf die ,,1‘-Seite bezogene Magnetisierungsstrom. Fiihrt man jetzt eine 
Impedanzmatrix Z’ ein mit den von den Strémen 7, . . . 4, 74, 24) durchflossenen Impedanzen, so 
ergibt sich, wenn man in jeden solchen Zweig analog zu AGNI, Kap.2 auch noch eine EMK e’ 

einfiihrt, die Beziehung 


uv =Z'’' +e’. (33) 
Hierbei ist also 
ae) ) ) 
ONNEce. 
/ . . . . . 
i O 7 0 1 -2es) 00 0 (34) 
Or” Ove 0 ann O 
OF “Or were: oO 0 219 10 


Zwischen den ungestrichenen Zweispannungen und -str6men einerseits und den Maschen- 
spannungen und -strémen andererseits bestehen aufgrund der Topologie des in Bild 5b wieder- 
gegebenen Graphen noch die folgenden Beziehungen 


u® = C,u (35) 
und 
i=Cr (36) 


mit 


No 
+- 

(o) e 

-- 

i Ae} 


JO mA NWN HW 
| 
a 
| 
par 
je) (s) 
oo) KS) 
Oo So 
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Man beachte die Entkopplung der Zweige 1, 2; 6, 7 einerseits und der Zweige 3, 4, 5; 8, 9, 10 
anderseits, welche sich im unteren Teil der C-Matrix durch eine entsprechende Kastchenbildung 
mit Nullmatrizen auBerhalb der Hauptdiagonale ausdriickt (in Gl. (37) durch Strichelung an- 
gedeutet). 


Aus den Gln. (30), (31), (33), (35) und (36) ergibt sich die Beziehung 


w=CwW2ZW/ CP+C.We’'. (38) 
Auch hier ist der Vektor der Maschenspannungen gleich null 
ei OF (39) 


so daB man erhalt 
uo C, WZ W; Cio + C,W e = 0 
oder 
GW ZW Ci=—C.We'. (40) 
Die in Kap. 2.21 beschriebene Maschenmethode ist in dieser Methode enthalten. Die Zweig- 
impedanzmatrix mit Kopplungskoeffizienten ist namlich 


ZN ITN Se (41) 


Da8 diese Matrix mit derjenigen nach Gl. (9) nicht exakt iibereinstimmt, liegt natiirlich daran, daB 
dort die Streuimpedanzen in den entsprechenden Vierpol mit einbezogen sind. Ferner ist 


e= Wer". (42) 
Als Maschenimpedanzmatrix ergibt sich 
Za, Vin Gre. Goe, C; (43) 


womit man eine Ubereinstimmung mit Gl. (14) erhalt. 
Hier soll nun aber ein anderer Weg gegangen werden. Man kann namlich auch zusammen- 
fassen 
CG W => G ) bzw. WwW G => Cc* (44) 
und hat dann den Vorteil, daB Z’ als Diagonalmatrix (die sich im Digitalrechner einfacher spei- 
chern 1aBt) erhalten bleibt. Gl. (40) geht hiermit tiber in 


CrZeGris a Gres (45) 
JEON eee GG: (45a) 


Diese Gleichung ist natiirlich wieder dieselbe Gleichung wie (14). Beide Gleichungen sind ja 
aus (40) hervorgegangen, nur die Klammersetzung bei der Ausrechnung ist eine andere. Da- 
durch, daB man W* mit C zusammengefaBt hat, erhalt man eine modifizierte Indizenzmatrix, 
in welcher das Ubersetzungsverhiltnis bereits beriicksichtigt ist 


(=A ) O ) ©) 
o +1 O 6) 
O Ol ae O 
) Oo --1 0) 
C= O 0) On ak (46) 
—1 —1 ©) ) (0) 
ail U7 3) 10 0) 
O (e) Oo —1 O 
O O o +1 —1 
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Sie geht aus der Matrix C nach Gl. (37) dadurch hervor, daB man die 8. Zeile mit wis) durch- 
multipliziert und zur 7. Zeile hinzuzahlt und alsdann die 8. Zeile streicht. Man kann C auch 
unschwer aus Bild 5b direkt ablesen. DaB jetzt die beiden Teile des Netzes nicht mehr ent- 
koppelt sind, zeigt sich in dem unteren Teil der Matrix dadurch, daB eine Partitionierung und 
Kastchenbildung mit Nullmatrizen auBerhalb der Hauptdiagonale hier nicht mehr méglich ist. 
Der untere Teil von C enthalt jetzt zwei Parallelreihen von nichtverschwindenden Elementen. 
Dies entspricht einer kettenartigen Struktur des Netzes. Im Sonderfall des 1:1-Transformators, 
also mit w7,s) = 1, geht die Inzidenzmatrix C in diejenige einer einfachen Kettenschaltung 
ohne Transformator tiber. Die beiden Querzweige 7 und 8 sind dabei in einen Zweig zusammen- 
geschrumpft, der die Magnetisierungsimpedanz enthalt. 


Stellt man also die Maschengleichungen nach Gl. (45) auf, so hat man bei digitaler Berech- 
nung zwei Vorteile: 1. Z’ ist Diagonalmatrix, 2. C, die man iibrigens auch ohne Multiplikation auf- 
stellen kann, hat pro Transformator eine Zeile weniger als C. 


Eine Anwendung dieser Methode findet man bei H. W. HALE und J. B. Warp [78]. 


3.2 Zwangsbedingungen in der Knotenpunktsmethode 


Um Zwangsbedingungen in der Knotenpunktsmethode einfiithren zu konnen, muB man jetzt 
beachten, daB die Admittanzmatrix des halb-idealen Transformators mit Magnetisierungs- 
impedanz nicht existiert, dafiir aber die Admittanzmatrix des halb-idealen Transformators mit 
Streuimpedanz. Wir verlegen also die Streuimpedanz in die Primar- und Sekundarzweige des 
Transformators. Die Magnetisierungsimpedanz kann jetzt allerdings nicht mehr in der tiblichen 
Form beriicksichtigt werden, da die nachtragliche Hinzufiigung einer Querimpedanz zu einer 
anderen Ersatzschaltung fithren wiirde. Da man ohnehin in den meisten Fallen auf eine Beriick- 
sichtigung der Magnetisierungsimpedanz verzichtet, soll sie hier einfach weggelassen werden. 


Generator Last Transformator Leitung Last 


Bild 6. Beriicksichtigung eines halb-idealen Transformators in der cinphasigen Darstellung eines Energietibertragungssystems 
entsprechend Bild 4. durch LHinfithrung von Zwangsbedingungen (Knotenpunktsmethode oder Schnittmengenmethode). 
a Originalschaltbild, b zugehoriger orientierter Graph mit vollstandigem Baum und Systeme unabhingiger Zweige. 


Bild 6a zeigt die Ersatzschaltung eines Netzes nach Bild 4 (ohne Magnetisierungs-Impedanz in 
der Ersatzschaltung des Transformators). Der diesem Netz zugeordnete Graph ist in Bild 6b 
wiedergegeben. Die Unterscheidung von unabhangigen Zweigen und Baumzweigen ist hier noch 
ohne Bedeutung. Sie erhalt erst bei der Anwendung der Schnittmengenmethode einen Sinn. Wir 
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beginnen mit der Beriicksichtigung des Ubersetzungsverhaltnisses in den Stromgleichungen und 
erhalten 


ene 
ro 
= 
<) 


I 


in Matrizen 


wee (47) 


Auch hier existiert wieder fiir die inversen Gleichungen der Spannungen eine Beziehung mit der 
zu W* konjugiert-transponierten Matrix W, und zwar 


mM = U4 
U, = U,— Wes, Us 
“= Uy 
Ca Ug 
in Matrizen 
/ 
ue = W,u (48) 
mit 
4 ) Omreso 
O 1 Ome 0 
—1 
Vy ll Seal acs) el oe (40) 
O ) pete iO 
fe) e) O 1 


Wahrend das Gleichungssystem (47) — von der Zwangsbedingung zwischen Primar- und Se- 
kundarstrom abgesehen — nur Identifizierungen enthalt, wird in (48) eine neue Spannung und 
zwar die auf die Primarseite bezogene Streuspannung gema8 der Gleichung 


Yt 5 


Up = Wy — Wea, 's) Ma = be 22, 3) — WE, 8) #9 2, 8) = 12 22,8) + WR» WES) b2 2, 9 
= 15 (ee, 3) + w~? 4 3) = 22 2,8) (50) 
eingefiihrt. Man kann nun wieder analog zu AGNI, Kap. 3 in jeden Zweig eine EMK e’ einfiihren 
und erhalt die Beziehung 


i 2 ie 
oder 
f= Ze! (uv —e) = Yuu =e’) (51) 
mit 
(cor ro) ) on 5} 
fo) 20,3) 0 O 
Yo =| 0 0) Cis 1G (52) 
0 0) Cute ee 


oe 
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Zwischen den ungestrichenen Zweigstromen und -spannungen einerseits und den Knotenpunkts- 
strémen und -spannungen andererseits bestehen aufgrund der Topologie des in Bild 6b wieder- 
gegebenen Graphen die folgenden Beziehungen 


it=Ki (53) 
und 
u.= Kus (54) 
mit 
es 
| C1, 02) 43d So Cen ee 
ne 
Kia Shs +1 +1 +1 O° BY cele ae eee (55) 
2 |—1 —1 ) ) o —1 ) ) 
3 | ) o —1 —1 ) —1 ) 
4 ) 0.4 adie epee 0) Og 


Aus den Gln. (47), (48), (51), (53) und (54) erhalt man schlieBlich 


PK WEY? Wik UK ay ce: (56) 
und mit 
it=0 (57) 
ein Gleichungssystem fiir die Knotenpunktsspannungen ut 
K.W* YoW2Kp ul = KV ea cee (58) 


Die Matrix des Systems ist singular, denn es ist wieder jede Gleichung die negative Summe aller 
iibrigen Gleichungen. Damit ist eine Gleichung entbehrlich. Ferner ist auch hier wieder jede 
Spannung #, bis auf eine willkiirliche Konstante bestimmt. Wir kénnen also z. B. Spannung 
u~ = 0 setzen (Spannung des Bezugsknotenpunkts). Streichen wir auch noch die %-te (dem Be- 
zugsknotenpunkt zugeordnete) Gleichung, so entspricht dies bei K der Streichung der %-ten 
Zeile. Nennen wir wieder die so entstandene reduzierte Knotenpunktszweiginzidenzmatrix K, 
so erhalt man in der gleichen Weise wie in Kap. 3 von AGNI das System, das genau die erforder- 
lichen Gleichungen enthalt, nadmlich 


KW Ye WYK Ue = NY re: (59) 
oder 
Y? Ue = AC Wy oee. (59a) 
Diese Gleichung entspricht mit 
2 VV * YAWe (60) 
und 
e° = Wyre (61) 


Gl. (26a). Auch hier unterscheidet sich die Methode der Zwangsbedingungen von der Methode 
eingefiihrter Kopplungskoeffizienten nach Kap. 2.22 nur dadurch, da8 die Klammern beim 
Ausrechnen des Gleichungssystems (59) anders gesetzt werden. Um die Diagonalgestalt von Y’ 
zu erhalten, bilden wir eine modifizierte Inzidenzmatrix gemaB 


K=KW bw. K=KW. (62) 
Fiir unser Beispiel lautet diese Matrix 
+1 1—wes o +1 +1 +1 +1 
7 —1 —1 ) o —1 ) 0) 
K= . : 6 
O a Was) == oO 10) ==! O ( 3) 
O (o) +1 —1 (0) o —1 
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Die Matrix K geht aus K dadurch hervor, daB man die 3. Spalte mit — ®(e') multipliziert, 
zur 2. Spalte hinzuaddiert und alsdann die 3. Spalte streicht. Die Eigenschaft aller K-Matrizen, 
daB die Summe aller Elemente in jeder Spalte verschwindet, bleibt fiir die K-Matrizen erhalten. 
Hierdurch erhalt man eine Knotenpunkts-Zweig-Inzidenzmatrix, in welcher das Ubersetzungs- 
verhaltnis der beiden Zweige 2 und 3 zueinander eingearbeitet ist. Man kann natiirlich auch 
leicht K direkt aus Bild 6b ablesen. Setzt man W(2, 3) = 1, So verschmelzen die beiden Zweige 2 
und 3 zu einem einzigen Zweig, der die Gesamtstreuimpedanz enthilt. 


Die Methode der Zwangsbedingungen hat auch hier wieder zwei Vorteile 1. Y’ ist Diagonal- 


matrix, 2. K, die man auch direkt ohne Multiplikation aufstellen kann, hat pro Transformator 


eine Spalte weniger als K. Das Gleiche gilt natiirlich auch fiir K. 


Das aufzuldsende Gleichungssystem kann nunmehr geschrieben werden 


KY’ Kut =K*Y’ e’ (64) 
oder 
Y Ut=KY’e’ (64a) 
mit 
Y= KY’K: (65) 


3.3 Zwangsbedingungen in der Schnittmengenmethode 


Die Einfiihrung von Zwangsbedingungen in der Schnittmengenmethode bereitet jetzt keine 
Schwierigkeiten mehr. Fiir die Aufstellung der Gleichungen kann Bild 6 zugrundegelegt werden, 
wobei die gleiche Matrix Y’ verwendet wird und auch die gleiche Matrix W. Zur Aufstellung der 
Schnittmengen kann man den in Bild 6b eingezeichneten vollstandigen Baum zugrundelegen. 
Man wird bei der Berechnung dual zur Maschenmethode und analog zur Knotenpunktsmethode 


H W = H zusammenfassen. Das aufzulisende Gleichungssystem lautet dann 


Y? u2 =H*Y'e’ (66) 
wobei 
Y? =H*Y’H, (67) 
und 
H=HW. (68) 


3-4 Hinreichende Bedingungen fiir die Auflésbarkeit der Maschen-, Knoten- 
punkts- und Schnittmengengleichungen bei Einfiihrung von Zwangsbedingungen 


Dasich die Methoden, welche die Transformatoren durch Einfiihrung von Zwangsbedingungen 
beriicksichtigen, immer auf den Fall der gekoppelten Zweige nach Kap. 2 zuriickfiihren lassen, 
kann die Frage nach hinreichenden Bedingungen fiir die Auflésbarkeit der Gleichungen leicht 
beantwortet werden. 


Im Falle der Maschengleichungen bilden wir 
Z.=W.2We (41) 
und wenden hierauf Satz A von Kap. 2.2 an. 
Im Falle der Knotenpunkts- oder Schnittmengengleichungen bilden wir 
Y=WEyY w, (60) 


und wenden hierauf Satz B von Kap. 2.2 an. 
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Man kann aber auch noch einfachere Kriterien erhalten und zwar dadurch, daB man Kri- 
terien aufstellt, in welchen die Diagonalmatrix Z’ bzw. Y’ zugrundegelegt wird. Wenn man nun 
analog zu Kap. 5.1 von AGNI vorgeht, so hat man also auch hier wieder nachzuweisen, da unter 
gewissen Voraussetzungen 


Zi, = —(Z° + Z}*) (69) 


positiv-definit ist, und somit nichtverschwindende Eigenwerte 4, von Z° zur Folge hat. Hierdurch 
wird 


det Z°= 1,40. (70) 


t=1 


Die Positiv-Definitheit von Z?, folgt genau dann, wenn bei nicht sémtlich verschwindenden Ma- 
schenstrémen die gesamte im Netz verbrauchte Wirkleistung positiv ist. Hierzu bilden wir 


P = Reif? Z° 9 = — (j0* Zo jo + 12 Z* j°*) = if* Zio e* Ge Ze CI. (71) 
Nun ist aber wegen der Symmetrie von Zj 
De. fo) 
Zi Red poh | eee (72) 


— 


O 7%, —7,2—7 


t = Anzahl der Transformatoren 


Die gesamte im Netz verbrauchte Wirkleistung ist nun aber auch gleich der Summe aller in den 
einzelnen Zweigen verbrauchten Wirkleistungen. Sieht man von den Zweigen, welche zu Trans- 
formatoren gehoren, ab, so gehGren die restlichen Zweige entweder zu den unabhangigen Zweigen 
oder zu den Baumzweigen. Die Transformatoren kénnen einen unabhangigen Zweig im Primar- 
zweig (,,1°‘-Zweig), im Sekundarzweig (,,w(¢,»)‘‘-Zweig), im Primar- und Sekundarzweig haben 
oder beide Zweige sind Baumzweige. Dain den unabhangigen Zweigen die Maschenstréme flieBen 
und diese unabhangig gewahlt werden kénnen, kann die Summe aller in den einzelnen Zweigen 
verbrauchten Wirkleistungen folgendermaBen geschrieben werden 


z—T 
70 . of . cA o. J 4 / / VF 
P=i* CFR CP =i, Ri= * hie ae 2 Di(1) Tum l2,./? aE (2) irri ca |, : ty 2 (3) Tw I, 
= le le le 
(73) 

mit 

. z2-—-T = 

Vy a Da Cuy ty D 

ab 


Die drei Summationen haben hierbei folgende Bedeutung 


a) = Summation iiber alle unabhangigen Zweige w und der unabhangigen Primar- 
zweige w der Transformatoren, jedoch ausgenommen solche unabhangigen Pri- 
marzweige mw, deren Sekundarzweige unabhangige Zweige der Transformatoren 
sind. 


2) = Summation tiber alle unabhangigen Sekundarzweige yu der Transformatoren. 


2/3) = Summation tiber alle Baumzweige w und von Transformatoren nur tiber die Primar- 
zweige, ausgenommen solche Baumzweige yu von Transformatoren, deren Sekundir- 
oder Primarzweige unabhangige Zweige sind (reine Baumzweige!). 


Es kommt jetzt darauf an, die ersten beiden Summen positiv-definit in den Strémen 2, der un- 
abhangigen Zweige zu machen. Hierzu ist ausreichend, daB 


1. die 7,,,, bzw. 77; dieser Summen positiv sind, : 


aie 
2. die w=) ~ 0 (dies sei automatisch erfiillt!). 
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Die iibrigen 7,,,, das sind die Realteile der reinen Baumzweig-Impedanzen, seien nicht-negativ. 
In diesem Fall wird die Summe der beiden Summen 24) und Xie) positiv und nur dann gleich 
null, wenn alle Maschenstréme verschwinden, und ferner Xg nicht negativ. Damit ist die Positiv- 
Definitheit der Matrix Z?, gesichert und det Z° ¢ 0, und es gilt somit der Satz: 


Damit das System der Maschengleichungen 
Le Wie Coe Ge (= Ce. (45) 


eines Netzes mit halb-idealen Transformatoren mit Magnetisierungs-Impedanzen eindeutige 
Lésungen hat, ist hinreichend, daB in irgendeinem zugrundegelegten System unabhangiger 
Zweige des Netzes die Impedanzen und in den Transformatoren, deren Primar- oder/und 
Sekundarzweige unabhangige Zweige sind, die Magnetisierungsimpedanzen positive 


| Realteile besitzen und die iibrigen Impedanzen keine negativen Realteile. 


Hierzu dual gilt auch der folgende Satz: 


Damit das System der Knotenpunktsgleichungen 


Vora K YG us =K Ye! (64) 
und ebenso auch das System der Schnittmengengleichungen 
Ye Teh Veh he (66) 


eines Netzes mit halb-idealen Transformatoren mit Streuimpedanz eindeutige Lésungen hat, ist 
hinreichend, daf in irgendeinem zugrundegelegten System von Baumzweigen des Netzes 
die Admittanzen und in den Transformatoren, deren Primar- oder/und Sekundiarzweige 
Baumzweige sind, die Streuadmittanzen positive Realteile besitzen und die tibrigen Ad- 
mittanzen keine negativen Realteile. 


Damit ist die Problematik der Auflésbarkeit der Gleichungssysteme fiir die drei Methoden 
mit modifizierten Inzidenzmatrizen hinreichend geklart. 


Zusammenfassung 


Die drei grundlegenden Methoden fiir die Berechnung von Zweipolnetzen mit Inzidenzmatrizen 
kann so erweitert werden, daB auch Transformatoren und andere Vierpole eines Netzes beriick- 
sichtigt werden kénnen. Lat sich ein Transformator (oder auch ein anderer Vierpol) mit Hilfe 
von Zweipolen darstellen, so bendtigt man keine Erweiterung der Theorie. Fiir beliebige Trans- 
formatoren mit reellem Ubersetzungsverhialtnis (auch ideale und halbideale) existiert immer eine 
Ersatzschaltung, die ausschlieBlich aus Zweipolen zusammengesetzt ist, solern man auf die 
Eigenschaft der Potentialtrennung verzichten kann (dain diesem Fall die Potentialtrennung nicht 
beriicksichtigt werden kann, handelt es sich bei diesen Ersatzschaltungen eigentlich nur um 
Ersatzschaltungen des Spartransformators). Da in der einphasigen Darstellung bzw. bei Dar- 
stellung in symmetrischen Komponenten von Energietibertragungssystemen es sehr haufig 
vorkommt, da8 Transformatoren mit komplexem Ubersetzungsverhaltnis auftreten, ist dieser 
Weg nicht gangbar. Ein Transformator mit komplexem Ubersetzungsverhaltnis kann auf keinen 
Fall durch ein Zweipolnetz dargestellt werden, da solche Transformatoren eine nichtsymmetrische 
Impedanz- bzw. Admittanzmatrix haben, die Impedanz- bzw. Admittanzmatrix eines Zweipol- 
netzes hingegen grundsatzlich immer symmetrisch ist. Ein einfacher Zweiwicklungstransformator 
dieser Art muB also hier notwendigerweise als Vierpol eingefiihrt werden, sofern man nicht das 
Inzidenzmatrizenschema modifiziert. Um solche Transformatoren als Vierpole einzufiihren, 
braucht man nur die entsprechenden Kopplungselemente au8erhalb der Hauptdiagonale der 
Zweig-Impedanz- bzw. Admittanzmatrizen zu besetzen und hierauf die drei Methoden anzu- 
wenden. Selbstverstandlich lassen sich nach dieser Methode auch Mehrwicklungstransformatoren 
und auch andere 2u-Pole in das Netz einfiihren, sofern die entsprechenden Impedanz- bzw. 
Admittanzmatrizen existieren. Allerdings kenn man nach dieser Methode nur halb-ideale Trans- 


Cail a 9h BRS <3 y 
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formatoren einfiihren, da fiir ideale Transformatoren Z- NOC 


matrix existiert. Das Gleiche gilt auch fiir die dritte Méglichkeit lod 
Inzidenzmattizen, welche auf den Fall eingefithrter Vierpole matickeelntye Se 
die digitale Behandlung hat man jetzt den Vorteil die Zweigadmittanzmatrix, die auBerl 
Hauptdiagonale jetzt wieder nur noch mit Nullen besetzt ist, in einfacher Weise na 
Vektor speichern zu kénnen. Damit sind alle Fragen fiir die systematische Berech: 
Spannungen, Strémen und Leistungen bei vorgegebenen Spannungsquellen (dies ist der 
wichtigere Fall) gentigend erlautert. Der Fall vorgegebener Stromquellen bedeutet hier n 
triviale Erweiterung auf die hier nicht naher eingegangen werden soll. Die Ermittlun, 
gewissen aus Netzen abgeleiteten Matrizen (Impedanz-, Admittanz-, gemischte Matrizen, 
lustmatrizen) wird in einer weiteren Arbeit behandelt. ae 
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